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Chapter  1 


Executive  Summary 

1.1  Significant  Contributions 

Neutralization  of  bacterial  aerosol  releases  is  critical  in  countering  bioterrorism  As 
a  possible  bacterial  aerosol  neutralization  method  that  avoids  the  use  of  chemicals, 
we  investigate  the  mechanical  instabilities  of  the  bacterial  cell  envelope  in  air  as  the 
bacteria  pass  through  aerodynamic  shocks.  To  carry  out  this  fundamental  investiga¬ 
tion  a  novel  experimental  impactor  system  is  designed  and  built  to  simultaneously 
create  a  controlled  and  measured  shock  and  to  collect  the  bacteria  after  they  pass 
through  the  shock.  In  the  impactor  system,  the  aerosol  flows  through  a  converging 
nozzle,  perpendicular  to  a  collection  surface  that  has  an  orifice  through  which  the 
shocked  bacteria  enter  a  deceleration  tube.  Both  experimental  measurements  of  the 
pressure  in  the  impactor  system  at  multiple  points  and  computational  fluid  dynamics 


I 


simulations  are  used  to  characterize  quantitatively  the  shocks  created  in  the  impactor. 
Specifically,  the  developed  computational  model  describes  the  evolution  of  both  the 
gas  and  particle  velocity  and  temperature  in  the  impactor  system  to  determine  the 
forces  exerted  on  the  bacterial  aerosol  as  they  pass  through  the  shock.  The  results 
indicate  that  the  developed  computational  model  predictions  compare  well  with  the 
experimental  pressure  measurements. 

The  bacteria  experience  relative  deceleration  because  of  sharp  velocity  changes  in 
the  aerodynamic  shock  created  in  the  experimental  impactor  system.  Computational 
model  results  indicate  that  vegetative  E.  coli  cells  require  a  critical  acceleration  of 
3.0  x  108  m/s2  compared  to  3  9  —  16  x  109  m/s2  for  B.  atropheus  spores  to  break-up. 
Computed  accelerations  in  the  impactor  system  reach  5.9  x  109  m/s2  which  predict 
break-up  of  both  vegetative  cells  and  spores.  Experimental  findings  demonstrate 
that  aerosolized  E.  coli  cells  and  B.  atropheus  spores  that  pass  through  aerodynamic 
shocks  created  in  the  experimental  impactor  system  are  an  order  of  magnitude  less 
likely  to  retain  viability  than  those  that  pass  through  the  impactor  at  conditions 
which  do  not  lead  to  the  generation  of  an  aerodynamic  shock,  and  therefore,  do  not 
reach  the  critical  acceleration  required  for  break-up.  However,  B.  atropheus  spores 
that  do  not  break-up  have  a  higher  likelihood  to  remain  live  compared  to  E.  coli 

The  novel  impactor  system  is  introduced  in  Chapter  2  of  Appendix  A.  Specif¬ 
ically,  the  advantages  of  using  the  impactor  system  are  highlighted  in  relation  to 
studying  the  effects  of  aerodynamic  shocks:  simultaneous  collection  and  control  of 
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shock  strength.  The  gas  dynamics  of  the  impactor  system  is  studied  with  both  ex¬ 
perimental  measurements  and  computational  model  simulations.  The  Mach  number 
of  the  aerodynamic  shock  is  shown  to  be  a  function  of  the  operating  conditions  as 
well  as  the  geometry  of  the  impactor.  The  experimental  measurements  verify  the 
computational  gas  dynamics  model  which  is  used  in  Chapter  3  of  Appedix  A. 

Chapter  3  of  Appendix  A  explores  the  motion  of  the  bacteria  as  they  travel  through 
the  impactor  system:  specifically,  the  deceleration  forces  that  act  on  the  bacteria  as 
they  passe  through  the  aerodynamic  shock.  These  forces  are  computed  for  the  differ¬ 
ent  geometries  and  operating  conditions  of  the  impactor  system  and  are  compared  to 
the  Rayleigh-Taylor  break-up  model.  Knowing  the  critical  Rayleigh- Taylor  accelera¬ 
tions  needed  to  break-up  a  cell  and  the  decelerations  achieved  in  the  impactor  system, 
we  are  able  to  make  predictions  on  which  cells  survive. 

In  Chapter  4  of  Appendix  A,  the  predictions  obtained  from  the  previous  chapters 
arc  tested  by  passing  aerosolized  E.  coli  vegetative  cells  through  the  experimental 
impactor  system.  Furthermore,  a  stochastic  random  walk  model  is  used  to  determine 
the  collection  efficiency  at  different  experimental  operating  conditions.  The  model 
developed  in  Chapter  3  of  Appendix  A  is  confirmed  with  the  experimental  viability 
measurements  obtained  using  a  fluorescence  flow  cytometer. 

Chapter  5  of  Appendix  A  provides  more  experimental  data  to  confirm  the  com¬ 
putational  predictions  using  B.  atropheus  spores  which  is  a  more  relevant  test  aerosol 
for  defense  applications. 
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Chapter  6  of  Appendix  A  summarizes  the  main  conclusions  of  the  dissertation 


(Appendix  A). 

Pages  vi  through  xvii  of  Appendix  A  contains  a  list  of  figures  with  descriptions. 
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Appendix  A 


P.  Sislian  Dissertation 


The  appendix  contains  the  dissertation  of  Patrick  Sislian  which  was  submitted  to  the 
University  of  Calfornia,  Los  Angeles  on  June  2  of  2010.  The  dissertation  includes 
ther  required  information  needed  in  the  final  technical  report  including  a  description 
of  the  equipment,  experimental  procedures  and  computational  codes. 

Section  2.2  describes  the  novel  impactor  system  design  and  includes  a  drawing  of 
the  parts. 

Section  2.3.1  outlines  the  Pitot  experimental  pressure  measurements  used  to  mea¬ 
sure  the  Mach  number  in  the  impactor  system.  Sections  4.2  and  5.2  describe  the 
experimental  procedures  and  setup  used  to  assess  the  effect  of  an  aerodynamic  shock 
on  E.  coli  vegetative  cells  and  B ,  atropheus  spores,  respectively. 

Section  2.3.2  details  the  gas  dynamic  equations  and  the  code  framework  used  to 
solve  these  equations.  Section  3.2  describes  the  equations  of  motion  for  and  their 


VI 


computational  implementation.  Section  4.3.3  describes  the  stochastic  particle  track¬ 
ing  model  used  in  FLUENT  to  calculate  the  collection  efficiency  of  the  bacterial  cells. 
Sample  codes  are  included  in  Appendix  B. 
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2.1  Schematic  of  the  impactor  system.  (1)  Converging  nozzle  with  an  exit 
diameter  of  either  0.5  mm  or  1  mrn,  (2)  and  (3)  a  flat  surface  with  a 
0.5  mm  or  1  mm  hole  combined  with  a  screw  that  accommodates  the 
collection  substrate  make  up  the  deceleration  tube  (these  parts  can  be 
replaced  with  a  Pitot  tube  to  measure  the  stagnation  pressure),  (4) 
spacers  allow  variation  of  the  plate  to  nozzle  distance  from  0  mm  to 
2  ram,  (5)  support  holding  the  deceleration  tube  and  flat  surface  at 
fixed  distance  from  the  nozzle,  (6)  exit  chamber,  and  (7)  Pitot  tube 
which  can  replace  parts  2  and  3  for  stagnation  pressure  measurements.  10 

2.2  Experimental  Mach  number  before  the  shock  (M i)  calculated  using 
Eq.  2.4  compared  to  measurements  from  Powell  [42].  M\  is  plotted  as 
a  function  of  xjd .  The  operating  conditions  are  the  same  for  both  sets 
of  measurements  with  Po  =  3.72  atm  and  Pi  —  1  atm  (Pj/Po  =  0.27). 

The  geometries  are  different .  17 
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2.3  Two  sets  of  two  Pitot  pressure  measurements  in  the  developed  im- 


pactor  system  presented  as  a  function  of  x/d  for  different  upstream 
and  downstream  pressures;  the  ratio,  Pi/Pn,  remains  fixed.  (Set  1) 

Pi/Po  =  0.1  with  (P0(atm),  Pi(atm))  =  {(2.5,  0.25).  (5.0, 0.50)},  and 
(set  2)  P\/Pq  =  0.2  with  (P0(atm),  P\{atm))  =  {(1.0, 0.20),  (2.5, 0.50)}.  18 

2.4  Experimental  Mach  number  before  the  shock  (M\)  in  the  impactor  sys¬ 
tem  of  Fig.  2.1  calculated  using  Eq.  2.4  as  a  function  of  x/d.  Curves 
represent  different  operating  conditions,  represented  by  the  downstream 
to  upstream  pressure  ratio,  P\/ Pq.  [Pi/Po,  Po,  Pi]  =  (0)  0.05,8.5  atm, 

0.4  atm  (O)  0.05,5  atm, 0.25  atm  (0)  0.16,2.5  atm, 0.4  atm  (A)  0.27, 

3.72  atm,  1  atm  (x)  0.25,1  atm, 0.25  atm  (o)  [0.32,1  atm, 0.32  atm] 

(+)  0.50, 1  atm, 0.5  atm .  19 

2.5  Computational  domain  and  boundary  conditions .  24 

2.6  Contour  plots  of  the  Mach  number  at  different  x/d  for  P0  =  1  atm  and 
Pi  =  0.14  atm.  The  z-axis  is  the  distance  along  the  centerline  from 

the  nozzle  outlet  and  is  given  in  meters .  26 

2.7  Contour  plots  (top  plots)  and  centerline  fluid  properties  (bottom  plots) 
for  (a)  Mach  number,  (b)  static  temperature  and  (c)  static  pressure. 

The  z-axis  is  the  distance  along  the  centerline  from  the  nozzle  outlet 


and  is  given  in  meters.  The  impactor  is  simulated  at  x/d  —  1.5  with 
Po  =  1  atm  and  Pi  =  0.14  atm .  28 
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2.8  (a)  Pitot  pressure  versus  xjd  at  P0  =  3.72  atm  and  Pi  =  1  atin  us¬ 
ing  both  experimental  and  computational  data  of  the  experimental 
impactor  system.  It  can  be  seen  that  the  Pitot  pressures  are  in  agree¬ 
ment  for  very  short  distances  but  diverge  as  xjd  increases,  (b)  Mach 
number  versus  xjd  for  the  same  conditions.  (□)  Maximum  center- 
line  Mach  number  extracted  from  computational  gas  dynamic  field, 

(O)  Mach  number  calculated  using  Eq.  2.4  and  computational  Pitot 
pressure  data,  (A)  Mach  number  calculated  using  experimental  Pitot 
pressure  measurements .  29 

2.9  (a)  Computational  (solid)  vs.  experimental  (dotted)  maximum  center- 
,  line  Mach  number  (Mi)  at  different  Pi/P0:  (□)  P\/Pq  =  0.14  (solid) 
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3.2  (a,  b)  Centerline  particle  and  gas  velocity  as  a  function  of  the  distance 
from  the  nozzle,  (c,  d)  Mean  free  path  of  the  gas  molecules  as  a  function 
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3.3  (a,  b)  Centerline  particle  and  gas  velocity  as  a  function  of  distance  from 
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x/d  —  1.4,  P0  =  5  atm  and  Pi  =  0.70  atm .  41 

3.4  Bacterial  particle  acceleration  (relative  to  gas  flow)  as  a  function  of 
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3.5  Maximum  centerline  acceleration  of  the  particle  as  a  function  of  nor¬ 

malized  distance  between  the  nozzle  and  the  plate  {x/d)  for  four  dif¬ 
ferent  (Pq,  Pi)  pairs .  46 
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4.1  Schematic  of  the  experimental  system.  (1)  Meinhard  nebulizer  with 
a  concentric  nozzle  mixing  a  liquid  bacterial  suspension  feed  and  a 
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to  0.6  inin,  (7)  support  holding  the  deceleration  tube  and  fiat  surface 

at  fixed  distance  from  the  nozzle,  and  (8)  exit  chamber .  53 
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4.3  Computational  domain  of  the  impactor.  The  dashed  lines  indicate 
different  virtual  surfaces.  These  surfaces  are  used  to  calculate  the 
particle  distribution  at  various  cross-sections  of  the  impactor  shown  as 

S-l,  S-2,  S-3  and  S-4 .  59 

4.4  (a,  b)  Centerline  particle  and  gas  velocity  and  (c,d)  Bacterial  particle 
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from  the  nozzle  exit.  Impactor  operating  conditions:  (a,c)  x  —  0.11 

and  (b,d)  x  —  0-98 .  64 

4.5  (a,b)  Turbulent  kinetic  energy  (m2/s2)  contours  of  the  gas  phase  in  the 
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4.6  Representative  trajectories  of  5  particles  released  from  the  inlet  of  the 
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ciencies  as  measured  by  the  number  of  particles  entering  the  collection 
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4.7  Distribution  of  particles  along  different  surfaces  (s-1,  s-2,  s-3,  s-4)  la¬ 
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4.8  (left)  Side-scatter  versus  forward-scatter  plot  of  all  particles  in  the  flow 
cytometer.  The  G2  gate  indicates  the  particle  sizes  which  represent  in¬ 
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(FLH-2)  versus  green  fluorescence  (FLH-1)  of  gated  cells  in  the  left  side 
plots.  The  L  gate  represents  the  live  cells  and  the  D  gate  represents 
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4.9  Fraction  of  live  cells  in  different  collected  samples,  (a)  All  experiments 
start  with  a  measurement  of  the  fraction  of  live  cells  in  the  cell  culture 


(susp),  the  negative  control,  and  in  the  nebulized  cells  (neb).  Error 
bars  show  the  95%  confidence  interval  using  a  normal  distribution  for  5 
measurements  per  sample,  (b)  The  cells  that  pass  through  the  impactor 
system  and  then  collected  in  the  collection  tube  at  x  —  0.98,  x  —  0.50, 
and  x  =  0.11  with  P$  =  1  atm.  The  ethanol  (EtOH)  suspension  is  the 
positive  control  Error  bars  show  the  95%  confidence  interval  using  a 
normal  distribution  for  15  measurements  per  sample . 

4.10  Fraction  of  cells  that  remain  alive  after  the  effect  of  impactor  operation 
(/;)  and  after  exposure  to  ethanol.  Live  cells  in  these  samples  are 
compaied  to  live  cells  in  the  negative  control  by  accounting  for  losses 
before  the  impactor  and  for  the  collection  efficiency  of  the  impactor. 
Error  bars  show  the  95%  confidence  interval  using  a  normal  distribution 
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5.1  Schematic  of  the  experimental  system.  (1)  Meinhard  nebulizer  with 
a  concentric  nozzle  mixing  a  liquid  bacterial  suspension  feed  and  a 
dispersion  gas  (/V2),  (2)  Nebulization  chamber  to  collect  large  droplets, 

(3)  Converging  nozzle  with  an  exit  diameter  of  d  —  0.5  mm,  (4)  and 

(5)  a  flat  surface  with  a  0.5  mm  hole  combined  with  a  screw  that 
accommodates  the  collection  substrate  make  up  the  deceleration  tube, 

(6)  1.4  mm  spacer  fixes  the  spacing  x  between  the  impactor  and  nozzle 
to  0.6  min,  (7)  support  holding  the  deceleration  tube  and  flat  surface 
at  fixed  distance  from  the  nozzle,  and  (8)  exit  chamber.  Major  losses 

in  the  experimental  system  are  shown  in  italic .  84 
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5.3  Bacterial  particle  acceleration  (relative  to  gas  flow)  as  a  function  of 
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5.4  Representative  trajectories  of  5  particles  released  from  the  inlet  of  the 
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5.5  Computational  domain  of  the  impactor.  The  two  dashed  lines  indicate 
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collection  efficiency  is  calculated  by  dividing  the  number  of  particles 
passing  S-2  by  that  of  S-l . . .  99 

5.6  (left)  Side-scatter  versus  forward-scatter  plot  of  all  particles  in  the 
flow  cytometer.  The  S  gate  indicates  the  particle  sizes  which  represent 
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5.10  Fraction  of  cells  that  remain  alive  after  the  effect  of  inipactor  opera¬ 
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ABSTRACT  OF  THE  DISSERTATION 
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Neutralization  of  bacterial  aerosol  releases  is  critical  in  countering  bioterrorism. 
As  a  possible  bacterial  aerosol  neutralization  method  that  avoids  the  use  of  chemicals, 
we  investigate  the  mechanical  instabilities  of  the  bacterial  cell  envelope  in  air  as  the 
bacteria  pass  through  aerodynamic  shocks.  To  carry  out  this  fundamental  investiga¬ 
tion,  a  novel  experimental  impactor  system  is  designed  and  built  to  simultaneously 
create  a  controlled  and  measured  shock  and  to  collect  the  bacteria  after  they  pass 
through  the  shock.  In  the  impactor  system,  the  aerosol  flows  through  a  converging 
nozzle,  perpendicular  to  a  collection  surface  that  has  an  orifice  through  which  the 
shocked  bacteria  enter  a  deceleration  tube.  Both  experimental  measurements  of  the 
pressure  in  the  impactor  system  at  multiple  points  and  computational  fluid  dynamics 
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simulations  are  used  to  characterize  quantitatively  the  shocks  created  in  the  impactor. 
Specifically,  the  developed  computational  model  describes  the  evolution  of  both  the 
gas  and  particle  velocity  and  temperature  in  the  impactor  system  to  determine  the 
forces  exerted  on  the  bacterial  aerosol  as  they  pass  through  the  shock.  The  results 
indicate  that  the  developed  computational  model  predictions  compare  well  with  the 
experimental  pressure  measurements. 

The  bacteria  experience  relative  deceleration  because  of  sharp  velocity  changes  in 
the  aerodynamic  shock  created  in  the  experimental  impactor  system.  Computational 
model  results  indicate  that  vegetative  E.  coll  cells  require  a  critical  acceleration  of 
3.0  x  108  m/s2  compared  to  3.9  —  16  x  10°  m/s2  for  B.  atropheus  spores  to  break-up. 
Computed  accelerations  in  the  impactor  system  reach  5.9  x  109  m/s2  which  predict 
break-up  of  both  vegetative  cells  and  spores.  Experimental  findings  demonstrate 
that  aerosolized  E .  colt  cells  and  B.  atropheus  spores  that  pass  through  aerodynamic 
shocks  created  in  the  experimental  impactor  system  are  an  order  of  magnitude  less 
likely  to  retain  viability  than  those  that  pass  through  the  impactor  at  conditions 
which  do  not  lead  to  the  generation  of  an  aerodynamic  shock,  and  therefore,  do  not 
reach  the  critical  acceleration  required  for  break-up.  However,  B.  atropheus  spores 
that  do  not  break-up  have  a  higher  likelihood  to  remain  live  compared  to  E .  coli. 
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Chapter  1 


Introduction 

1.1  Background  on  bacterial  aerosol  neutralization 
by  aerodynamic  shocks 

Natural  bacterial  aerosol  is  present  both  in  indoor  and  outdoor  environments  at  back¬ 
ground  concentrations.  The  suspended  bacterial  particles  usually  range  in  size  from 
1  (single  cells)  to  10  fim  (multiple  cells  or  cells  associated  with  debris)  [58].  With 
respect  to  bioaerosol  concentrations,  the  comparison  of  concentration  data  from  dif¬ 
ferent  environments  is,  in  general,  very  difficult.  For  outdoor  environments,  the  re¬ 
ported  literature  values  of  bioaerosol  concentration  range  from  80  cfu/ml  in  rural 
areas  [5,  27,  47,  57]  to  850  cfu/ml  in  urban  areas  [5,  15,  39,  41],  where  cfu  refers  to 
colony  forming  units.  Indoor  environments  usually  contain  less  bioaerosol.  Different 
sources  contribute  to  the  background  in  different  environments  such  as  soil,  water, 
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plants  and  human  activity  (anthropogenic  aerosol).  In  addition  to  bioaerosols  present 
in  indoor  and  outdoor  environments  at  background  concentrations,  intentionally  pro¬ 
duced  anthropogenic  bacterial  aerosols  can  potentially  be  used  as  bioterrorism  agents, 
which  are  added  to  natural  background  concentrations. 

Characteristics  of  a  bacterial  cell,  such  a s  degree  of  hydration,  growth  medium 
(e.g.,  agar  vs.  broth),  Gram  stain  (+ve  vs.  -ve),  and  metabolic  state  (vegetative  vs. 
spore)  affect  its  survival  [54,  55].  Survival  rates  of  different  bacterial  cells  are  also 
affected  by  exposure  to  different  air  environments  [54,  55].  Temperature,  relative  hu¬ 
midity,  ultraviolet  light,  and  atmospheric  pollutants  are  some  factors  that  contribute 
to  cell  viability  loss.  These  environmental  factors  minimally  injure  spores  compared 
to  both  Gram-negative  and  Gram-positive  vegetative  cells.  Intentionally  disseminat¬ 
ing  spores  such  as  Bacillus  anthracis  (anthrax)  can  cause  much  greater  harm  than 
their  vegetative  forms  because  of  their  extended  survival  in  air.  Consequently,  it  is 
critical  to  study  the  effectiveness  of  neutralization  techniques  on  spores  for  defense 
applications. 

With  respect  to  sensing,  different  approaches  for  biological  aerosols  have  been 
studied  and  can  be  classified  into  four  categories:  (1)  Nucleic  Acid-based,  (2)  structure- 
based,  (3)  chemistry-based  and  (4)  function-based  [58].  Assuming  the  detection  of 
these  agents  is  possible,  current  methods  are  not  effective  in  neutralizing  the  aerosol 
cloud  at  its  source  and  mostly  include  containment  of  the  aerosol  release  in  order 
to  minimize  human  exposure.  These  methods  are  limited  to  indoor  environments 
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and  employ  techniques  such  as  high  efficiency  particulate  air  (HEPA)  filtration,  elec¬ 
trostatic  precipitation,  steam  condensation,  ultraviolet  (UV)  inactivation  as  well  as 
diverting  airflows  [58]. 

The  use  of  aerodynamic  shocks  holds  promise  as  an  alternative  way  to  neutralize 
the  bacterial  aerosol  at  the  source.  The  shock  can  be  applied  at  the  point  of  release 
both  in  indoor  and  outdoor  air  with  the  advantage  of  avoiding  chemicals.  Some 
experimental  data  exist  in  the  literature  [25,  40,  56]  on  the  use  of  shock  waves  to 
neutralize  bacteria;  however,  these  studies  are  limited  to  bacterial  liquid  suspensions 
or  powders  and  do  not  assess  the  effects  of  a  shock  in  the  aerosol  phase.  We  aim  at 
understanding  the  fate  of  the  bioaerosol  particles  when  exposed  to  shocks  of  different 
velocity,  temperature  and  pressure  gradients  in  air. 

Although  detailed  data  for  the  effect  of  shocks  on  bioaerosol  particles  are  lacking, 
aerodynamic  shocks  have  been  shown  to  cause  lysis  of  bacteria  suspended  in  the  liquid 
phase  [40,  56].  Teshima  et  al.  (1995)  [56]  studied  the  biomechanical  effects  on  E.  coll 
by  attaching  a  test  tube  with  the  bacterial  suspension  at  the  end  of  a  shock  tube. 
The  peak  pressure  of  the  shock  was  varied  from  3  MPa  to  14  MPa  with  100  pulses 
and  bacterial  lysis  was  quantified  by  measuring  leaked  phenylalanine  dehydrogenase 
activity  [56].  At  14  MPa,  15%  of  the  E.  coh  cells  were  destroyed  [56].  The  researchers 
mentioned  the  importance  of  the  pulse  width  but  failed  to  account  for  its  effect  [56]. 
In  earlier  studies,  Lundbeck  et  al.  (1963)  [40  used  a  shock  peak  pressure  of  ~  6  MPa 
with  variable  number  of  pulses  to  neutralize  bacteria  suspended  in  the  liquid  phase 
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The  shocks  were  produced  using  a  piston  fitted  to  a  stainless  steel  tube  filled  with 
suspension.  The  study  showed  a  loss  of  viability  of  E.  coli  vegetative  cells  but  did  not 
show  any  significant  effect  on  B.  subtilis  spores  [40].  The  authors  also  demonstrated 
that  an  increase  in  concentrations  results  in  an  increase  of  bacterial  viability.  Our 
approach  to  shock  bioaerosol  particles  in  an  impactor  system  avoids  concentration 
effect  because  individual  cells  pass  through  the  shock;  this  allows  understanding  the 
mechanism  of  single  bacterial  lysis  by  an  aerodynamic  shock. 

More  recently,  Horneck  ct  al.  (2001)  [25]  show-cd  that  only  0.01%  of  dried  B. 
subtilis  powder  (mounted  between  two  quartz  plates)  survived  a  shock  peak  of  32  GPa 
with  1  pulse.  However,  the  temperature  of  the  cells  in  the  sample  reached  250  °C  [25] 
and  it  was  not  possible  to  separate  the  effect  of  the  shock  from  that  of  temperature. 

The  aerodynamic  shock  represents  a  sharp  change  in  fluid  properties  (velocity, 
temperature,  pressure).  Due  to  its  inertia,  a  bacterial  particle  suspended  in  air  will 
not  be  able  to  immediately  adjust  to  a  sharp  drop  in  gas  velocity  and  will  therefore 
experience  a  relative  deceleration.  Rayleigh-Taylor  instabilities  in  the  bacterial  mem¬ 
brane  arise  if  the  bacterial  particle  is  accelerated  (or  decelerated)  in  a  perpendicular 
direction  to  the  bacteria-gas  interface  [7.  28].  We  have  computed  the  critical  accel¬ 
erations  and  Weber  numbers  needed  for  several  bacterial  aerosols  to  break-up  (see 
Chapter  3).  Higher  surface  tension  and  smaller  diameter  of  the  bacterium  increase 
the  critical  accelerations  needed  for  break-up  (see  Chapter  3).  An  impactor  system 
has  been  built  and  modeled  to  achieve  different  bacterial  accelerations  through  a 
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well-defined  aerodynamic  shock  (see  Chapter  2).  After  passing  through  the  shock, 
the  bacterial  aerosol  is  collected  at  low  velocity  (<  10  m/s)  to  avoid  break-up  at  the 
point  of  impaction  [53/  In  the  impactor  system  the  aerosol  flows  through  a  converg¬ 
ing  nozzle,  perpendicular  to  a  collection  surface.  The  collection  surface  has  an  orifice 
through  which  the  shocked  bacteria  enter  a  deceleration  tube  where  their  velocity 
is  reduced  to  avoid  neutralization  by  impaction.  The  work  focuses  on  the  develop¬ 
ment  of  a  computational  gas  (flow  and  temperature)  and  particle  motion  model  for 
the  impactor  system  capable  of  predicting  the  effect  of  the  shock  on  the  bacterial 
aerosols.  Experimental  pressure  measurements  provide  information  about  the  flow- 
field  in  the  impactor  and  are  compared  with  computational  model  predictions.  Sharp 
changes  in  velocity,  temperature  and  pressure  across  the  shock  induce  instabilities  on 
the  bacterial  membrane.  The  magnitude  of  the  deceleration  of  the  bacterial  particle 
is  calculated  and  is  compared  with  theoretical  values  needed  for  the  neutralization  of 
different  bacteria. 

Experimental  measurements  of  cell  viability  in  the  impactor  system  provide  arc 
compared  to  the  theoretical  predictions.  Both  vegetative  cells  (Chapter  4)  and  spores 
(Chapter  5)  are  aerosolized  and  collected  after  operating  the  impactor  at  different 
operating  conditions.  A  mass  balance  calculation  reveals  the  percentage  of  cells  that 
break-up  after  the  operation  of  the  impactor. 
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1.2  Dissertation  objectives  and  structure 


The  objective  of  this  dissertation  is  to  study  for  the  first  time  the  effect  of  aero¬ 
dynamic  shocks  on  bacterial  aerosols.  The  dissertation  provides  a  theoretical  and 
computational  framework  that  allows  for  the  prediction  of  operating  conditions  and 
geometry  of  the  impactor  system  that  lead  to  bacterial  aerosol  neutralization.  A 
detailed  understanding  of  the  mechanisms  of  loss  of  cell  viability  is  critical  for  both 
defense  applications  as  well  as  for  the  design  of  bacterial  aerosol  samplers.  For  defense 
applications,  the  results  allow  for  the  design  of  aerodynamic  shocks  that  are  able  to 
neutralize  specific  bacterial  aerosols.  For  the  design  of  collection  devices,  the  results 
indicate  the  operating  conditions  which  will  prevent  break-up  of  cells  if  the  goal  is  to 
collect  live  cells  for  later  analysis. 

The  structure  of  this  thesis  is  as  follows:  first,  the  novel  impactor  system  is  in¬ 
troduced  in  Chapter  2.  Specifically,  the  advantages  of  using  the  impactor  system  are 
highlighted  in  relation  to  studying  the  effects  of  aerodynamic  shocks:  simultaneous 
collection  and  control  of  shock  strength.  The  gas  dynamics  of  the  impactor  system  is 
studied  with  both  experimental  measurements  and  computational  model  simulations. 
The  Mach  number  of  the  aerodynamic  shock  is  shown  to  be  a  function  of  the  operating 
conditions  as  well  as  the  geometry  of  the  impactor.  The  experimental  measurements 
verify  the  computational  gas  dynamics  model  which  is  used  in  Chapter  3. 

Chapter  3  explores  the  motion  of  the  bacteria  as  they  travel  through  the  im¬ 
pactor  system:  specifically,  the  deceleration  forces  that  act  on  the  bacteria  as  they 
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passe  through  the  aerodynamic  shock.  These  forces  are  computed  for  the  different 
geometries  and  operating  conditions  of  the  impactor  system  and  are  compared  to  the 
Rayleigh-Taylor  break-up  model.  Knowing  the  critical  Rayleigh-Taylor  accelerations 
needed  to  break-up  a  cell  and  the  decelerations  achieved  in  the  impactor  system,  we 
are  able  to  make  predictions  on  which  cells  survive. 

In  Chapter  4,  the  predictions  obtained  from  the  previous  chapters  are  tested  by 
passing  aerosolized  E.  coh  vegetative  cells  through  the  experimental  impactor  system. 
Furthermore,  a  stochastic  random  walk  model  is  used  to  determine  the  collection 
efficiency  at  different  experimental  operating  conditions.  The  model  developed  in 
Chapter  3  is  confirmed  with  the  experimental  viability  measurements  obtained  using 
a  fluorescence  flow  cytometer. 

Chapter  5  provides  more  experimental  data  to  confirm  the  computational  pre¬ 
dictions  using  B .  atropheus  spores  which  is  a  more  relevant  test  aerosol  for  defense 
applications. 

Chapter  6  summarizes  the  main  conclusions  of  this  dissertation. 
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Chapter  2 

Impactor  System:  Design  and 
Aerodynamics 

2.1  Introduction 

Aerodynamic  shocks  are  created  in  an  impactor  system  described  in  detail  in  this 
chapter.  The  impactor  geometry  offers  various  advantages  over  other  systems  where 
aerodynamic  shocks  are  created.  First,  the  impactor  allows  for  simultaneously  passing 
the  bacterial  aerosol  through  an  aerodynamic  shock  and  achieving  their  collection. 
Second,  the  impactor  system  isolates  the  neutralization  effect  of  the  aerodynamic 
shock  from  that  of  impaction,  dehydration  and  temperature.  These  two  advantages 
allow'  for  detailed  experimental  analysis  of  various  bacterial  aerosols  (see  Chapters  4 
and  5).  A  computational  fluid  dynamics  model  is  developed  in  this  chapter  which 
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will  be  used  in  the  break-up  model  of  the  bacterial  particles  in  Chapter  3. 


2.2  Design  of  the  impactor  system 

The  impactor  system  consists  of  a  plate  designed  for  the  collection  of  the  bacterial 
aerosol  (see  Fig.  2.1;  parts  2  &  3)  which  is  placed  perpendicular  to  the  gas  flow 
emerging  from  a  converging  nozzle  (see  Fig.  2.1;  part  1).  Isentropic  flow  theory  of  an 
ideal  gas  predicts  a  critical  downstream  (Pi)  to  upstream  ( Pq )  pressure  ratio  below 
which  the  flow  at  the  exit  of  the  nozzle  is  sonic  [38,  48],  according  to  the  following 
equation: 

Xcrit  =  A/P0  =  [2/(7  +  l)r/(7-1),  (2.1) 

where  7  is  the  heat  capacity  ratio  (7  —  1.4;  Xcrxt  =  0-53  for  air).  The  generated 
impinging  flow  on  a  plate  from  a  converging  nozzle  operating  under  sonic  conditions 
results  in  the  creation  of  a  standoff  shock  whose  properties  can  be  changed  by  vary¬ 
ing  impactor  geometry  and  operating  conditions  [2,  13,  14,  29,  42].  In  traditional 
impactors  [21,  221,  the  aerosol  is  collected  by  impaction  on  a  flat  surface.  Stewart  et 
al.  (1995)  [53]  showed  that  bacteria  collected  by  impaction  are  injured  [53];  therefore, 
assessing  damage  induced  by  the  standoff  shock  is  impossible  using  current  impactor 
designs.  To  avoid  loss  of  viability  by  impaction,  the  flat  surface  of  the  developed 
impactor  system  has  an  opening  (see  Figure  2.1;  part  2)  that  acts  as  a  virtual  surface 
through  which  the  bacteria  enter  a  stagnant  gas  deceleration  tube  (see  Figure  2.1; 
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Figure  2.1:  Schematic  of  the  impactor  system.  (1)  Converging  nozzle  with  an  exit 
diameter  of  either  0.5  mm  or  1  mm,  (2)  and  (3)  a  flat  surface  with  a  0.5  mm  or 
1  mm  hole  combined  with  a  screw  that  accommodates  the  collection  substrate  make 
up  the  deceleration  tube  (these  parts  can  be  replaced  with  a  Pitot  tube  to  measure 
the  stagnation  pressure),  (4)  spacers  allow  variation  of  the  plate  to  nozzle  distance 
from  0  mm  to  2  mm,  (5)  support  holding  the  deceleration  tube  and  flat  surface  at 
fixed  distance  from  the  nozzle,  (6)  exit  chamber,  and  (7)  Pitot  tube  which  can  replace 
parts  2  and  3  for  stagnation  pressure  measurements. 

parts  2  and  3).  Biswas  and  Flagan  (1988)  have  designed  a  similar  tube  to  prevent 
particle  reentrainement  [4].  The  final  velocity  of  the  bacterium  is  related  to  the  length 
of  the  tube  {51)  and  is  calculated,  assuming  that  the  air  in  the  deceleration  tube  is 
stagnant,  using  a  force  balance  of  the  form: 


dvV  dvP  rp  , 

mv~J£  =  =  Fd  + 


(2.2) 
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the  drag  force  Fp  is  given  by  Stokes-Cuiniingharn  law  as  follows: 


Fd  = 


—  3tt  /j,dpvp 

~C, 


(2-3) 


where  mp  is  the  bacterium  mass,  vp  is  the  bacterium  velocity,  dp  is  the  bacterium 
diameter,  x  is  the  axial  coordinate,  g  is  the  gravitational  constant,  and  n  is  the 
dynamic  viscosity  of  the  gas.  Cc,  the  Cunningham  slip  correction  factor,  accounts 
for  non-continuum  effects  (for  1  /zm  particles  Cc  =  1.164  [18]).  The  length  of  the 
deceleration  tube  (SI  =  (CcmPvp\x) / (37r fidp))  is  obtained  by  integrating  Eq.  2.2  (ne¬ 
glecting  the  effect  of  gravity)  between  the  tube  inlet  (x'  =  x)  and  the  impaction 
surface  at  the  end  of  the  deceleration  tube  ( x '  —  61  +  x).  The  length  of  the  tube, 
61,  varies  between  1  mm  to  4  mm  in  the  impactor  design  of  Fig.  2.1,  which  is  suffi¬ 
cient  to  collect  1  fj, m  particles  entering  the  deceleration  tube  at  velocities  ( vp\T. )  up 
to  680  m/s.  The  bacterium  experiences  deceleration  both  as  it  passes  through  the 
shock  and  in  the  boundary  layer  of  the  plate;  therefore,  a  tube  between  2  mm  and 
4  mm  in  length  is  sufficient  to  collect  bacteria  at  negligible  velocities.  More  detailed 
computational  modeling  of  the  impactor  system  (see  Section  3.2)  confirms  that  the 
bacterium  velocity  at  the  entrance  of  the  deceleration  tube  ( vp\x )  is  much  smaller 
than  vp\x  =  680  m/s  which  is  assumed  in  the  design  of  the  deceleration  tube.  A 
detailed  solution  of  Eqs.  2.2  and  2.3,  without  the  assumption  of  stagnant  gas,  reveals 
low  impaction  velocities  (<  10  m/s)  at  the  point  of  collection  (see  Section  3.2). 

The  deceleration  tube  (see  Fig.  2.1;  parts  2  and  3  combined)  is  held  on  a  support 
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(see  Fig.  2.1;  part  5).  The  distance  between  the  opening  of  the  deceleration  tube  and 
the  nozzle  outlet  (x)  can  be  adjusted  with  spacers  of  different  thickness  (see  Figure  2.1; 
part  4).  The  exit  chamber  (see  Fig.  2.1;  part  6)  has  two  outlets  (see  Fig.  2.1;  part  7); 
one  outlet  connects  to  a  vacuum  pump  and  the  other  connects  to  a  pressure  gauge  for 
Pitot  measurements.  The  upstream  pressure  (Pq)  and  the  downstream  pressure  (P\) 
are  measured  through  small  orifices  drilled  normal  to  the  side  walls  of  the  converging 
nozzle  and  of  the  exit  chamber,  respectively  (see  Fig.  2.1;  part  7).  A  valve  attached 
to  a  purge  stream  downstream  of  the  impactor  controls  P1?  whereas  a  nitrogen  tank 
with  a  pressure  regulator  controls  Pq. 

Depending  on  the  method  of  analysis,  different  collection  substrates  can  be  placed 
on  the  top  surface  of  the  screw  (see  Fig.  2.1;  part  3)  to  collect  the  shocked  bacteria.  A 
silicon  wafer  substrate  is  used  for  scanning  electron  microscopy  (SEM)  and  a  copper 
grid  is  used  for  transmission  electron  microscopy  (TEM).  Both  analysis  techniques 
reveal  qualitative  structural  data  of  the  collected  samples.  Biochemical  analysis  can 
also  be  carried  out  by  placing  a  black  polycarbonate  filter  on  the  top  surface  of 
the  screw  (see  Fig.  2.1;  part  3),  which  is  stained  during  post  processing  using  a 
LIVE/DEAD®BacLightTW  (cat.  number  L7012;  Life  Technologies  Corporation)  to 
assess  the  viability  of  the  collected  bacteria,  both  quantitatively  and  qualitatively. 
Experimental  results  and  analysis  of  the  collected  bacteria  are  discussed  in  detail  in 
Chapters  4  and  5. 

Parts  2  h  3  can  be  replaced  with  a  Pitot  tube  that  has  the  same  diameter  open- 
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ing  as  the  deceleration  tube  (0.5  mm  or  1.0  mm).  The  Pitot  tube  extends  beyond 
the  support  and  exit  chamber  (see  Fig.  2.1)  and  is  connected  to  a  pressure  gauge 
that  measures  the  stagnation  pressure.  The  alternate  Pitot  setup,  similar  to  the  ones 
used  in  previous  studies  [16,  45],  allows  experimental  characterization  of  the  standoff 
shock  (see  Section  2.3.1  below).  The  developed  system  therefore  has  the  advantage 
of  passing  the  bacteria  through  a  controlled  and  measured  aerodynamic  shock  and 
simultaneously  collecting  them  on  various  substrates  for  further  analysis.  The  simul¬ 
taneous  collection  eliminates  the  need  to  sample  the  bacterial  aerosol  compared  to 
shock  tube  geometries,  thereby  reducing  sampling  errors.  Both  the  operating  pres¬ 
sures  (Pq  and  Pi)  and  the  geometry  of  the  impactor  (nozzle  diameter,  d)  and  plate 
to  nozzle  distance,  x)  affect  the  properties  of  the  aerodynamic  shock  [13,  14].  Their 
effect  will  be  quantitatively  studied  below. 

2.3  Experimental  and  computational  studies  of  the 
impactor  system 

The  velocity,  pressure,  and  temperature  of  the  gas  in  the  impactor  system  for  dif¬ 
ferent  geometric  and  operating  conditions  are  important  parameters  in  determining 
the  forces  that  act  on  a  bacterium  as  it  passes  through  the  aerodynamic  shock.  Ex¬ 
perimental  measurements  are  compared  with  computational  fluid  dynamics  (CFD) 
model  predictions.  CFD  calculations  provide  velocity,  pressure  and  temperature  data 
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at  every  point  in  the  impactor  system  which  are  not  accessible  in  the  experimental 
impactor  system.  The  CFD  profiles  are  also  used  to  predict  the  particle  velocity  and 
temperature  in  the  impactor  system  (see  Chapter  3).  Therefore,  comparison  of  the 
experimental  and  computational  data  provides  realistic  predictions  of  the  eventual 
fate  of  different  bacterial  particles  as  they  pass  through  the  shock  and  are  collected 
in  the  system. 

2.3.1  Experimental  results:  Pitot  pressure  measurements 

One  of  the  directly  measured  variables  in  the  experimental  impactor  system  of  Fig.  2.1 
(see  Section  2.2)  is  the  Pitot  stagnation  pressure  (Pp)  after  the  shock.  While  a 
detailed  CFD  model  has  been  developed  for  the  experimental  impactor  system  (see 
Section  2.3.2  below),  in  order  to  be  able  to  compare  our  experimental  measurements 
with  results  available  in  the  literature,  we  first  develop  an  impactor  model  based  on 
appropriate  simplifying  assumptions  to  estimate  the  Mach  number  before  the  shock 
(Mi)  using  the  Pitot  pressure  measurement  ( Pp )  along  with  that  of  upstream  (Po) 
and  downstream  (Pi)  pressures. 

Specifically,  the  gas  is  considered  to  be  thermally  perfect  (i.e.,  the  enthalpy  ( h ) 
is  only  a  function  of  temperature;  h  =  h(T ))  and  caloncally  perfect  (i.e.,  the  heat 
capacity  (cp)  at  constant  pressure  is  not  a  function  of  the  temperature;  Cp  ^  cp(T)). 
Along  the  center  streamline,  the  flow  before  and  after  the  shock  is  assumed  to  be 
one  dimensional,  steady- state  and  isentropic  (i.e.,  adiabatic  and  frictionless).  These 
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conditions  lead  to  poi  ~  Po  and  po2  ~  Pp  where  poi  and  po2  are  the  stagnation 
(or  total)  pressures  before  and  after  the  shock,  respectively.  Both  Pq  and  Pp  are 
experimentally  measured  in  our  system  (see  Section  2.2).  The  flow  through  a  shock 
is  non-isentropic  (poi  ^  P02;  i  e.,  not  in  equilibrium).  Assuming  that  the  shock  is 
infinitely  thin  and  normal  to  the  center  streamline,  the  following  equation  can  be 
used  to  relate  the  stagnation  pressures  to  the  Mach  number  before  the  shock  [38,  48]: 


Pp  ~  P02 
Po  Poi 


1  + 


27 

7  +  1 


(Mi2  -  1) 


-1/(7— 1) 


(7  +  l)Mi2 


7/(7~l) 


(7-  l)A/,2  +  2 


(2.4) 


where  7  denotes  the  heat  capacity  ratio  and  M 1  is  the  Mach  number  before  the  shock. 
The  Mach  number  after  the  shock  (M2)  can  also  be  calculated  using  Mi  under  the 
above  flow  assumptions  using  the  following  formula  [38,  48]: 


M22 


1  +  l^Mi2 
7 Mt2  -  ^ 


(2.5) 


Although  the  experimental  pressure  measurements  do  not  provide  a  detailed  picture 
of  the  flow-field  in  the  impactor,  one  can  use  the  Mach  number  before  the  shock  as  a 
measure  of  the  velocity  change  (M\  —  M2)  over  the  thickness  of  the  shock  which  can 
be  estimated  to  be  on  the  order  of  one  mean  free  path  (A)  of  the  gas,  which  is  given 
by  [9]: 

A  =  2.15pT011/7^oi,  (2.6) 
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where  T0i  and  P0 1  are  the  static  temperature  and  pressure  before  the  shock,  re¬ 
spectively.  The  Mach  number  combined  with  the  shoek  thiekness  is  important  in 
determining  the  forces  exerted  on  the  bacterial  particles  (see  Sections  3.2  and  3.3). 

The  flow-field  assumptions  used  in  the  derivation  of  Eqs.  2.4  and  2.5  give  rise 
to  errors  in  determining  the  aetual  shock  Mach  number.  Detailed  CFD  simulations 
(see  Section  2.3.2)  can  be  used  to  eompute  a  better  estimate  of  the  Maeh  number 
and  will  therefore  be  used  in  the  development  of  a  bacterial  particle  break-up  model 
(see  Section  3.3).  Nonetheless,  the  experimental  Pitot  measurements  constitute  a 
basis  for  comparison  with  both  experimental  data  found  in  the  literature  and  with 
the  CFD  simulations.  The  experimental  impaetor  systems  found  in  the  literature  arc 
geometrically  different  from  the  impaetor  system  of  Figure  2.1  [16,  42,  45].  Although, 
impaetor  system  geometry  is  an  important  factor  that  affects  the  flow  characteristics 
of  the  impinging  jet,  a  comparison  is  still  useful  if  the  objective  is  to  determine 
whether  the  two  sets  of  measurements  (Pitot  pressure)  or  calculations  (Mach  number) 
fall  within  the  same  range.  Therefore,  even  when  the  operating  conditions  are  the 
same,  we  do  not  expect  the  measurements  and  calculations  to  agree  perfectly. 

Specifically,  Powell  calculated  an  oseillating  Maeh  number  as  a  function  of  the 
distance  between  the  nozzle  and  the  Pitot  tube,  x/d  [42].  Measurements  with  the 
current  system  (see  Fig.  2.1)  were  earried  out  to  mimie  Powell’s  operating  conditions 
( Pq  =  3.72  atm  and  Pi  =  1  atm)  (Fig.  2.2).  However,  the  geometries  differ  as  follows: 
(1)  d  =  7.9  mm  in  Powell’s  setup  [42]  compared  with  d  =  0.5  mm  in  the  impaetor 
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Figure  2.2:  Experimental  Mach  number  before  the  shock  (Mi)  calculated  using  Eq.  2.4 
compared  to  measurements  from  Powell  [42],  M\  is  plotted  as  a  function  of  x/d.  The 
operating  conditions  arc  the  same  for  both  sets  of  measurements  with  P0  =  3.72  atm 
and  P\  =  l  atm  (Pi/ Pq  =  0.27).  The  geometries  are  different. 

system  of  Fig.  2.1,  and  (2)  a  Pitot  tube  of  unspecified  geometry  without  the  flat  plate 
is  used  in  Powell's  case  [42],  thereby  avoiding  impinging  flow.  Due  to  this  fact,  in  the 
experimental  x/d  range  shown  in  Figure  2.2,  Powell’s  data  have  two  maxima  and  one 
minimum  whereas  in  the  experimental  measurements  of  the  impactor  system  such 
behavior  is  not  observed.  Both  measurements  though  predict  Mach  numbers  that  are 
in  the  same  range,  i.e.  Mi  is  between  1  and  2.75  (see  Figure  2.2). 

In  the  experimental  impactor  system  of  Fig.  2.1,  more  measurements  were  done 
with  geometries  of  varying  plate  to  nozzle  distance  (0.2  <  x/d  <  2.4;  limited  only 
by  the  spacer  thicknesses)  at  different  upstream  (1.0  atm<  P0  <  8.9  atm)  and  down¬ 
stream  pressures  (0.12  atm<  P\  <  0.9  atm).  We  compared  measurements  of  Pp/Pq  at 
constant  P\/Pq  for  different  upstream  and  dowmstream  pressures.  In  Figure  2.3,  the 
Pp/Po  measurements  are  identical  for  equal  P\/Pq  with  different  upstream  pressures. 
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These  results  indicate  that  the  main  variable  determining  the  ratio,  Pv/ Pq,  is  the  up¬ 
stream  to  downstream  pressure  ratio  instead  of  the  absolute  values  of  the  pressures. 
CFD  results  further  confirm  these  observations  (see  Section  2.3.2  below). 


Figure  2.3:  Two  sets  of  two  Pitot  pressure  measurements  in  the  developed  impactor 
system  presented  as  a  function  of  x/d  for  different  upstream  and  downstream  pres¬ 
sures;  the  ratio,  P\/Pq ,  remains  fixed.  (Set  1)  P\/Pq  —  0.1  with  (P0(atm),  P\{atm)) 
=  {(2.5, 0.25),  (5.0, 0.50)},  and  (set  2)  Pi/P0  =  0.2  with  (PQ(at.m),  Pi(citm))  = 
{(1.0,0.20),  (2.5,0.50)}. 

Fig.  2.4  shows  the  variation  of  the  Mach  number  (Mi)  with  x/d  for  different 
Pi/Pq.  The  Mach  number  increases  with  decreasing  P\/Pq  and  this  trend  is  more 
pronounced  at  higher  x/d  ( x/d  >  1)  which  is  observed  in  the  computational  results 
shown  in  Section  2.3.2  below  (see  Fig.  2.9-a).  Fig.  2.4  reenforces  the  findings  of  Fig.  2.3 
in  that  measurements  made  with  the  same  P\/Pq  but  different  P0  and  P\  yield  similar 
curves.  Specifically,  the  calculated  M\  for  P\f  Pq  —  0.05  is  identical  for  two  upstream 
pressures  P$  =  8.5  atm  and  5  atm.  The  highest  Mach  number,  calculated  using 
Equation  2.4,  with  the  current  setup  was  M\  —  3.93  ±  0.01  at  Pq  —  8.93  ±  0.01  atm. 
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Pi  =  0.41  ±0.01  atm  (P\/P0  =  0.05)  and  x/d  =  2.4.  The  errors  in  the  measurements 
are  computed  using  the  standard  deviation  of  3  independent  Pitot  measurements  and 
are  on  the  order  of  0.01  atm.  These  measurements,  along  with  the  CFD  predictions, 
are  essential  in  understanding  the  key  parameters  that  influence  bacterial  aerosol 
neutralization. 


Figure  2.4:  Experimental  Mach  number  before  the  shock  (Mi)  in  the  impactor  sys¬ 
tem  of  Fig.  2.1  calculated  using  Eq.  2.4  as  a  function  of  x/d.  Curves  represent 
different  operating  conditions,  represented  by  the  downstream  to  upstream  pressure 
ratio,  Pi/Pq.  [P\/Po,  Pq,  Pi]  =  (□)  0.05,8.5  atm,  0.4  atm  (Q)  0.05,5  atm, 0.25  atm 
(0)  0.16,2.5  atm, 0.4  atm  (A)  0.27,  3.72  atm.  1  atm  (x)  0.25,1  atm. 0.25  atm  (o) 
[0.32, 1  atm, 0.32  atm]  (+)  0.50, 1  atm, 0.5  atm 


2.3.2  Computational  results:  gas  flow  and  thermal  fields 


For  the  two-phase  aerosol  flow  (air-bacteria)  present  in  the  impactor  system,  the 
equations  for  gas  and  particle  dynamics  are  coupled.  At  low  aerosol  concentrations, 
however,  the  particle  dynamics  can  be  assumed  to  have  no  effect  on  the  gas  flow 
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field  [18],  and  the  equations  of  particle  velocity  and  temperature  can  be  decoupled 
from  those  of  the  gas  phase.  The  gas  dynamics  can  therefore  be  solved  independently 
and  the  results  are  then  used  to  determine  the  particle  trajectories  (see  Section  3.2 
below).  An  Eulerian  approach  is  used  to  solve  the  gas  flow  field.  Because  the  flow 
is  turbulent,  a  full  solution  to  the  Navier-Stokes  equations  requires  computationally 
intensive  Direct  Numerical  Simulation  (DNS).  In  this  work,  the  Reynolds-averaged 
Navier-Stokes  (RANS)  equations  are  solved  so  that  the  small-scale  turbulence  fluctu¬ 
ations  do  not  have  to  be  simulated  directly.  In  Reynolds  averaging,  variables  such  as 
velocity,  pressure  and  energy  are  decomposed  into  a  mean  and  fluctuating  components 
as  follows: 

0  =  0  +  0',  with  0  =  —  I  <f)dt ,  and  0  =0,  (2.7) 

A f  Jto 

where  0  is  a  scalar  property  of  the  fluid  (e.g.,  velocity  magnitude,  pressure,  tempera¬ 
ture,  energy  ...),  and  0  is  the  tiine-averaged  component  of  0  with  the  corresponding 
fluctuating  component  0'.  If  we  substitute  the  above  form  for  the  flow  variables  into 
the  instantaneous  Navier-Stokes  equations  and  take  a  time  average,  we  obtain  the 
equations  that  describe  the  gas  flow  field.  For  simplicity  and  to  provide  an  idea  of 
the  terms  included  in  these  equations,  we  provide  below  the  general  form  of  these 
equations  in  Cartesian  tensor  form  (after  dropping  the  overbars): 

|  +  (2.8) 
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where  p  is  the  density,  v  is  the  velocity,  p  is  the  static  pressure,  6  is  the  Kroneeker 
delta,  (i  is  the  viscosity,  and  t,  j,  and  l  are  the  indices.  To  close  the  form  in  Eq.  2.9, 
the  Reynolds  stress  term  —pv[v':  is  calculated  using  the  Boussinesq  hypothesis  [24]: 


-pv\v2  =  /i, 


dv\  dv., 
dxj  dxi 


2  /  dv[' , 

-3{pk  +  IMfrl)St,' 


(2.10) 


where  pt  is  the  turbulent  viscosity  and  k  is  the  turbulent  kinetic  energy.  Different 
turbulence  models  were  employed  and  compared  giving  similar  results  (see  Remark  1). 
Neglecting  the  effect  of  gravity,  the  turbulence  model  used  is  the  standard  k  —  e  model 
which  is  given  by 
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Cic  =  1.44,  C2c  =  1.92,  CM  =  0.09,  a*  -  1.0,  a,  =  1.3, 


where  e  is  the  rate  of  dissipation,  and  Gk  is  the  generation  of  turbulent  kinetic  energy 


due  to  mean  velocity  gradients.  <Jk  and  ae  are  the  turbulent  Prandtl  numbers.  Y\i  is 


the  dilatation  dissipation  term,  which  represents  the  contribution  of  the  fluctuating 


dilatation  in  compressible  turbulence  to  the  overall  dissipation  rate.  Ym  is  given  by 


Ym  =  2  peM2 


(2.13) 


where  Mt  is  the  turbulent  Mach  number,  defined  as 


Aft  =  y/k/tf- 


(2.14) 


where  a  (=  y/'yRT)  is  the  speed  of  sound  and  7  (= cp/cy )  is  the  ratio  of  specific  heat. 
The  final  conservation  law  solved  is  the  energy  equation  for  the  gas  which  is  given, 
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in  Cartesian  tensor  form,  as  follows: 
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where  E  is  the  internal  energy,  a  is  the  thermal  diffusivity,  cPg  is  the  heat  capacity 
at  constant  pressure,  and  T  is  the  static  temperature. 

A  numerical  solution  program  was  written  within  FLUENT®,  a  CFD  software,  to 
solve  the  above  equations  for  the  computational  domain  shown  in  Fig.  2.5.  Because 
the  impactor  is  axisymmetric,  a  cross  sectional  grid  that  covers  half  of  the  impactor 
system  is  used.  A  different  mesh  is  created  for  different  distances  between  the  nozzle 
and  the  impactor  wall  ( x/d ).  The  distance  shown  in  Fig.  2.5  is  x/d  =  4.0  with  a  nozzle 
diameter  of  0.5  mm.  The  grid  is  more  dense  (0.01  mm  without  adaption;  0.0025  mm 
with  adaption)  in  the  area  between  the  nozzle  outlet  and  the  deceleration  tube  inlet 
in  order  to  capture  the  details  of  the  shock.  Both  the  inlet  and  outlet  boundaries 
are  defined  by  stagnation  pressures  fixed  at  the  pressures  shown  in  Table  2.1  below 
and  temperatures  fixed  at  300  K  with  the  gas  considered  to  be  ideal.  The  walls  are 
defined  by  zero  heat  flux  (adiabatic)  and  no-slip  boundary  conditions.  As  for  the 
centerline,  wre  use  axisymmetric  boundary  conditions.  Full  multi-grid  (FMG)  initial¬ 
ization  is  used  with  the  first-order  upwind  scheme  to  obtain  a  convergent  solution  to 
the  momentum  equation.  Next,  the  energy  and  turbulence  equations  are  solved  with 
a  first-order  scheme  followed  by  a  second-order  upwind  scheme  to  obtain  convergent 
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solutions  with  residuals  of  continuity,  velocity,  k  and  e  of  less  than  10“4  and  energy 
residual  of  less  than  10“6.  Subsequently,  three  grid  adaptation  steps  are  applied  to 
refine  the  mesh  in  areas  with  high  pressure  gradients,  specifically  in  the  area  between 
the  nozzle  exit  and  the  deceleration  tube  inlet.  Li  et  al.  [34,  37,  35,  36]  used  a  simi¬ 
lar  approach  to  solve  for  a  high  velocity  oxygen-fuel  (HVOF)  thermal  spray  process 
involving  a  converging-diverging  nozzle.  A  parametric  study  with  different  upstream 


and  downstream  pressures  (see  Table  2.1)  is  conducted  to  compare  with  the  exper¬ 
imentally  calculated  Mach  numbers  (see  Figs.  2.8  and  2.9)  as  well  as  to  understand 
the  effects  of  shocks  with  different  properties  on  the  bacterial  aerosol.  Cases  1  and  4 
in  Table  2.1  have  different  upstream  pressures  as  boundary  conditions;  however,  the 
ratio  P\/Pq  is  constant.  Case  5  is  for  comparison  with  experimental  conditions  that 
are  similar  to  Powell’s  [42]. 

Theoretical  and  experimental  observations  indicate  the  presence  of  instabilities  in 
the  flow  structure  of  the  impinging  jets  in  the  presence  of  an  impaction  plate  [2,  29,  42]. 
Understanding  of  the  properties  of  an  under-expanded  free  (absence  of  impaction 
plate)  jet  is  essential  for  the  study  of  instabilities  for  an  impinging  jet  [42] .  In  a  free  jet, 
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Table  2.1  Inlet  and  outlet  pressure  conditions  for  parametric  analysis  of  the  impactor 
system  with  variable  nozzle  to  plate  distance  ( x/d ).  (x/d)c  refers  to  the  value  of  x/d 
where  instabilities  in  the  flow  arise,  xjd  values  are  represented  as:  (min  :  increment 
:  max). _ 


Case 

x/d 

Po  (atm) 

P\  (atm) 

Pi /a 

{x/d)c 

1  (base) 

(0.2 

0.1  : 

4.0) 

1.00 

0.14 

0.14 

1.5 

2 

(0.2 

0.1  : 

2.4) 

1.00 

0.25 

0.25 

0.9 

3 

(0.2 

0.1 

2.4) 

1.00 

0.32 

0.32 

0.9 

4 

(0.2 

0.1  : 

2.4) 

5.00 

0.70 

0.14 

1.6 

5  (for  comparison) 

(0.2 

0.1  . 

2.4) 

3.72 

1.00 

0.27 

- 

as  the  flow  emerges  from  the  nozzle,  it  expands  and  accelerates  to  supersonic  velocities 
creating  a  shock  after  which  the  gas  velocity  is  subsonic  [42Jt  The  subsonic  flow 
expands  further  to  regain  sonic  velocity  thereby  matching  the  nozzle  conditions  [42]. 
The  aerodynamic  cell  length  is  defined  as  the  length  of  the  region  between  the  nozzle 
exit  and  the  point  where  the  free  jet  flow  regains  a  near  nozzle  state.  Instabilities 
in  the  flow  structure  (e.g.,  recirculation  bubbles)  arise  when  an  impingement  plate  is 
placed  at  a  distance  equal  to  an  aerodynamic  cell  length  from  the  nozzle  [42j .  Jurcick 
et  al.  (1989)  29]  computationally  observed  these  instabilities  by  solving  the  time- 
dependent  Navier-Stokes  equations  for  P\/Po  =  0.16  and  x/d  >  1.7.  Recirculation 
bubbles  in  the  impingement  zone,  which  lies  at  the  center  of  the  flow-field  between 
the  shock  and  the  plate,  are  an  indication  of  instabilities  [2,  29].  Further  study  of  the 
structure  of  the  instabilities  is  outside  the  scope  of  this  work;  however,  calculating  the 
critical  x/d  [{x/d)c}  at  which  instabilities  arise,  using  the  developed  computational 
fluid  dynamics  model,  is  important  in  analyzing  the  particle  motion  in  the  impactor 
system  (see  Section  3.2  below).  Further  analysis  of  the  flow  structure  for  Pq  “  1  atm 
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Figure  2.6:  Contour  plots  of  the  Mach  number  at  different  x/d  for  Pq  =  1  atm  and 
P\  —  0.14  atm.  The  x-axis  is  the  distance  along  the  centerline  from  the  nozzle  outlet 
and  is  given  in  meters. 

and  P\  =  0.14  atm  (P\/Po  =  0.14)  reveals  the  appearance  of  recirculation  bubbles  at 
x/d  >  1.5,  which  is  consistent  with  the  range  x/d  >  1.7  reported  by  Jurcick  et  al. 
(1989)  [29].  The  recirculation  bubble  velocities  increase  as  x/d  increases.  As  P\/Pq 
decreases  recirculation  bubbles  form  at  lower  x/d  with  lower  velocities.  Furthermore, 
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as  can  be  seen  in  Fig.  2.6,  the  shock  structure  evolves  from  shocks  forming  away 
from  the  centerline,  to  stand-off  shocks,  to  diamond  shock.  At  xjd  <  0.5,  the  shocks 
form  away  from  the  centerline  along  the  impactor  walls  because  the  nozzle  area  is 
not  the  critical  flow  area.  As  xjd  increases,  a  stand-off  shock  forms  perpendicular  to 
the  centerline  with  increasing  Mach  numbers.  For  xjd  >  2.0,  a  shock  diamond  forms 
with  a  subsequent  second  shock.  For  the  different  operating  conditions,  our  study 
will  focus  on  0.5  <  xjd  <  ( x/d)c  where  the  shock  is  stable  and  along  the  path  that 
the  particle  follows  to  settle  in  the  deceleration  tube.  Table  2.1  refers  to  the  (x/d)c 
values  where  instabilities  arise  for  all  the  operating  conditions  studied. 

The  centerline  properties  of  the  fluid  are  typical  of  the  flow-field  that  a  particle 
experiences  before  settling  on  the  collection  plate  for  analysis.  In  Section  3.2  the 
particle  is  released  at  the  centerline.  Therefore,  the  study  will  focus  on  the  centerline 
properties  of  the  complex  flow  field  that  develops  in  the  experimental  impactor  of 
Fig.  2.1.  Fig.  2.7  shows  both  the  contour  and  centerline  temperature,  pressure  and 
Mach  numbers  for  xjd  —  1.5,  P0  =  1  atm  and  Px  =  0.14  atm.  As  the  gas  exits  the 
nozzle  at  x  =  0.062  m,  it  expands  reaching  supersonic  velocities  before  the  shock. 
Correspondingly  the  static  temperature  and  pressure  decrease  as  the  gas  expands.  In 
Figure  2.7,  the  shock  at  0.0625  m,  is  observed  as  a  sharp  change  in  the  fluid  properties, 
which  is  responsible  for  the  forces  experienced  by  the  bacteria  (see  Sections  3.2  and 
3.3  below).  In  the  contour  plots  of  Mach,  temperature  and  pressure,  a  standoff  shock 
perpendicular  to  the  direction  of  the  flow  can  be  observed.  The  Mach  contour  pattern 
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Figure  2.7:  Contour  plots  (top  plots)  and  centerline  fluid  properties  (bottom  plots) 
for  (a)  Mach  number,  (b)  static  temperature  and  (c)  static  pressure.  The  x-axis  is 
the  distance  along  the  centerline  from  the  nozzle  outlet  and  is  given  in  meters.  The 
impactor  is  simulated  at  x/d  =  1.5  with  Pq  =  1  atm  and  Pi  =  0.14  atm. 


is  similar  to  the  ones  observed  in  experimental  shadowgraphs  and  wave  diagrams  in 
the  literature  [2,  29]. 

In  the  CFD  model,  the  Pitot  pressure  is  the  centerline  static  pressure  at  the  end  of 
the  deceleration  tube  (see  Remark  3).  The  CFD  calculations  are  made  using  idealized 
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assumptions,  discussed  above,  such  as  perfect  gas,  axisymmetric  flow,  and  frictionless 
walls.  Whereas,  in  the  experimental  Pitot  measurements  these  assumptions  are  not 
made  and  the  pressure  is  measured  directly.  Therefore,  we  expect  differences  in  the 
measured  and  computed  values  of  the  Pitot  pressure.  In  Fig.  2.8-a,  the  computational 
results  yield  higher  Pitot  pressures.  However,  the  extracted  maximum  centerline 
Mach  numbers  represent  a  more  accurate  picture  of  the  gas  dynamics;  specifically, 
at  x/d  <  0.5.  Equation  2.4  is  for  supersonic  Pitot  pressure  measurements  and  does 
not  account  for  subsonic  conditions.  At  x/d  <  0.5  the  flow  along  the  centerline  is 


x/d  x/d 

(a)  (b) 


Figure  2.8:  (a)  Pitot  pressure  versus  x/d  at  Po  =  3.72  atm  and  P\  =  l  atm  using 
both  experimental  and  computational  data  of  the  experimental  impactor  system  It 
can  be  seen  that  the  Pitot  pressures  are  in  agreement  for  very  short  distances  but 
diverge  as  x/d  increases,  (b)  Mach  number  versus  x/d  for  the  same  conditions.  (□) 
Maximum  centerline  Mach  number  extracted  from  computational  gas  dynamic  field, 
(O)  Mach  number  calculated  using  Eq.  2.4  and  computational  Pitot  pressure  data, 
(A)  Mach  number  calculated  using  experimental  Pitot  pressure  measurements. 


subsonic  as  discussed  above  (see  Fig.  2.6).  Therefore,  Eq.  2.4,  which  yields  Mach 
numbers  >  1  at  x/d  <  0.5,  is  not  accurate.  Eq.  2.4  also  does  not  take  into  account 
the  formation  of  two  shocks  (see  Fig.  2.6)  after  certain  x/d  giving  rise  to  additional 
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errors  at  higher  x/d. 


Even  though  a  perfect  match  is  not  expected  as  discussed  above,  the  trends  in  the 
computational  maximum  centerline  Mach  numbers  are  similar  to  the  experimental 
results  calculated  using  Eq.  2.4  as  shown  in  Fig.  2.9-a.  The  Mach  number  increases 
with  decreasing  P\/Pq  and  the  effect  is  more  pronounced  after  x/d  of  about  1  For 
the  same  upstream  to  downstream  pressure  ratio  with  different  absolute  values  for 
the  pressures,  the  maximum  centerline  Mach  number  is  identical  (see  Fig.  2.9-b). 
This  fact,  is  also  confirmed  for  Mach  numbers  calculated  from  the  experimental  Pitot 
measurements  of  Fig.  2.3  as  well  as  the  Mach  numbers  shown  in  Fig.  2.4. 


(a) 


(b) 


Figure  2.9:  (a)  Computational  (solid)  vs.  experimental  (dotted)  maximum  centerline 
Mach  number  (Mi)  at  different  P\/Pq'-  (□)  P\/Pq  =  0.14  (solid)  and  0.15  (dotted), 
(O)  0.25.  (0)  0.32.  All  the  upstream  pressures,  Pq,  are  equal  to  1.  (b)  Computa¬ 
tional  Mi  vs.  x/d  at  constant  P\/Pq  =  0.14  for  different  upstream  and  downstream 
pressures. 


Although  temperature  effects  will  not  be  considered  in  the  discussion  of  the  forces 
experienced  by  the  particle  (see  Section  3.3  below),  we  acknowledge  that  this  effect 
exists  and  will  be  the  focus  of  further  study.  The  minimum  temperature  on  the  cen- 
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terline  decreases  as  P\/Pq  increases  (see  Fig.  2.10).  For  constant  P0/Pi,  the  minimum 
temperature  is  identical,  mimicking  the  Mach  trend  of  Fig.  2.9. 


Figure  2.10:  Minimum  centerline  temperature  of  the  gas  as  a  function  of  the  nor¬ 
malized  distance  between  the  nozzle  and  the  plate  {x/d)  for  the  4  different  (Pq,Pi) 
pairs. 

In  summary,  the  computational  and  experimental  gas  flow-field  analysis  yield  sim¬ 
ilar  trends  as  discussed  above.  The  absolute  values  differ  because  of  the  assumptions 
in  calculating  the  Mach  number  employed  in  the  developed  models.  The  more  de¬ 
tailed  CFD  model  will  be  used  in  determining  the  forces  exerted  on  the  bacterial 
particles  as  they  travel  through  the  experimental  impactor  system  of  Fig.  2.1. 

Remark  1  The  centerline  stagnation  pressures  at  the  end  of  the  deceleration  tube 
(Pitot  pressures)  were  compared  using  Spalart-Allmaras,  standard  k  —  e,  RNG  k  —  e, 
and  Realizable  k  —  e  models  and  the  results  were  found  to  have  a  maximum  deviation 
of  4%.  The  operating  conditions  for  the  comparison  were  x/d  =  0.6,  Pq  =  5.48  atm 
and  P\  =  0.08  atm.  Based  on  these  results,  the  standard  k  —  e  model  was  chosen  in 
all  the  simulations  discussed  in  the  study. 
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Remark  2  As  noted  m  Section  2.3.1,  the  shock  thickness  is  on  the  order  of  the  mean 
free  path  (X)  of  the  gas.  In  the  simulation  results  shown  in  Fig.  2.7,  the  centerline 
velocity  drop  occurs  over  a  range  of  about  42  fim  from  a  maximum  Mach  number  of 
2.52  to  less  than  1.  The  maximum  mean  free  path  along  the  centerline  is  0.51  fun 
which  is  much  less  than  the  distance  over  which  the  Mach  number  decreases.  CFD 
comparisons  are  based  on  flow  equations  based  on  the  assumption  of  continuity  and 
are  therefore  not  able  to  capture  non-continuum  shock  properties  such  as  the  shock 
thickness.  However,  this  has  a  minimal  effect  on  computing  the  forces  acting  on  the 
bacterium  as  can  be  seen  in  Section  3.2  below. 

Remark  3  All  the  reported  computational  Pitot  pressures  are  extracted  from  static 
pressure  data  along  the  centerline  and  at  the  end  of  the  deceleration  tube.  A  new 
mesh  was  also  created  to  capture  the  length  of  the  Pitot  tube  (7  cm)  of  Fig.  2.1-7. 
Calculations  using  the  new  Mesh  at  Pq  =  3.72  atm,  P\  =  1  atm  and  x/d  =  1.6  reveal 
a  Pitot  pressure  equal  to  that  at  the  end  of  the  deceleration  tube  (Pp  =  3.00  atm). 
The  Pitot  pressures  are  also  equal  by  comparing  the  centerline  pressure  to  the  average 
pressure  at  the  bottom  of  the  deceleration  tube.  Therefore,  reporting  the  centerline 
static  pressure  at  the  end  of  the  deceleration  tube  is  sufficient. 
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2.4  Conclusions 


An  impactor  system  was  developed  capable  of  simultaneously  passing  the  bacterial 
aerosols  through  a  controlled  shock  and  collecting  them  at  low  velocities.  The  Mach 
number  of  the  aerodynamic  shock  can  be  varied  by  changing  the  operating  conditions 
(P\,Pq  and  P\ / Pq)  and  the  geometry  (i /d.)  of  the  impactor  system.  A  maximum  Mach 
number  of  3  is  possible  with  the  current  design.  Both  experimental  measurements  and 
detailed  computational  fluid  dynamics  models  were  used  to  analyze  the  Mach  number 
of  the  shock  with  close  agreement.  The  computed  Mach  numbers  using  experimental 
Pitot  measurements  were  in  good  agreement  with  Powell’s  experimental  work  [42], 
given  the  geometrical  differences  between  the  two  systems. 
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Chapter  3 


Bacterial  Motion  and  Break-up 

3.1  Introduction 

Understanding  the  forces  acting  on  the  bacterial  membrane  as  a  vegetative  cell  or 
spore  passes  through  the  aerodynamic  shock  is  critical  and  determines  whether  the 
bacterium  will  remain  viable.  Therefore,  a  model  is  needed  to  assess  whether  the 
force  the  bacterium  experiences  as  it  goes  through  the  impactor  is  sufficient  to  cause 
loss  of  viability.  In  Chapter  2  we  developed  a  computational  model  to  understand  the 
gas  dynamics  and  this  model,  will  be  used  in  this  chapter  to  compute  the  deceleration 
experienced  by  the  bacterium  as  it  passes  through  the  shock. 

In  developing  a  model  to  understand  the  effect  of  the  shock,  we  assume  that 
break-up  of  the  cell  is  needed  to  achieve  cell  neutralization;  although  this  is  a  more 
stringent  requirement  than  loss  of  viability.  The  critical  force  needed  to  break-up  a 
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cell  is  calculated  using  the  Rayleigh-Taylor  instability  model  which  is  compared  to 


values  from  the  bacterial  motion  model  to  provide  predictions. 

3.2  Bacterial  motion  in  the  impactor  system 

The  bacterial  motion  in  the  impactor  system  is  computed  in  Lagrangian  coordinates 
utilizing  the  gas  flow  and  temperature  fields  calculated  using  the  CFD  model  of 
Section  2.3.2.  Due  to  the  very  low  bacterium  loading,  the  coupling  between  the  gas 
phase  and  the  bacterial  (particulate)  phase  can  be  considered  to  be  one-way,  i.e. , 
the  gas  flow  field  solution  is  obtained  by  neglecting  the  presence  of  bacteria  and  is 
subsequently  used  unmodified  in  the  calculation  of  the  bacterium  trajectories.  In 
the  bacterium  model  development,  it  is  also  reasonable  to  assume  that  coagulation  is 
negligible  and  that  particle  size  does  not  change  during  flight. 

The  particle  motion  in  the  impactor  system  is  governed  by  Newton’s  second  law 
and  is  given  by 


(3.1) 


where  mp  is  the  particle  mass,  vp  is  the  particle  velocity,  t  is  the  time,  and  Ft  are  the 


different  forces  that  can  act  on  the  particle  including  the  drag  force,  thermophoretic 
force,  gravitational  force,  the  basset  history  term,  etc.  For  typical  impactor  operat¬ 
ing  conditions,  the  major  force  acting  on  a  bacterial  particle  is  the  drag  force,  and 
the  other  forces  can  be  neglected.  As  a  result,  the  force  balance  of  Equation  3.1 
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applied  to  bacterial  particles  traveling  through  the  developed  impactor  system  takes 


the  following  form: 


(3.2a) 


dx 


(3.2b) 


where  vp  and  vg  are  the  velocity  of  the  bacterial  particle  and  gas,  respectively,  pg 
is  the  density  of  the  gas,  Av  is  the  projected  area  of  the  bacterium  on  the  plane 
perpendicular  to  the  flow  direction,  and  Co  is  the  drag  coefficient.  In  the  supersonic 
Hows  occurring  in  the  impactor  system  compressibility  and  non-continuum  effects 
have  to  be  considered  in  the  evaluation  of  the  drag  coefficient  Co-  Specifically,  Co 
is  a  function  of  both  particle  Reynolds  (Rep  =  pgdp \vg  —  vp\ /  pg)  and  particle  Mach 
( Mav  —  \vp  —  vg\/s\  s  is  the  speed  of  sound)  numbers.  For  the  simple  case  where  the 
ratio  of  particle  to  gas  absolute  temperature,  Tp/Tg ,  is  approximately  unity,  data  for 
Co  —  /(i?ep,  Map)  is  given  in  Clift,  Grace  and  Weber  [9].  At  each  grid  point  the  Rep 
and  Map  are  computed  and  a  corresponding  value  for  Co  is  evaluated. 

The  bacterial  particle  is  assumed  to  be  spherical  at  all  times.  The  density  (pp) 
of  bacterium  is  1  g/ml  [59].  The  dimension  of  bacterial  particles  varies  depending 
on  the  species  and  the  state  of  the  individual  bacterium.  B.  subtilis ,  for  example, 
has  a  width  of  0.7  —  0.8  pm  and  a  length  of  2  —  3  pm  [33].  Most  bacteria  diameters 
are  in  the  micrometer  range,  and  we  will  therefore  assume  that  they  are  spheres  of 
dp  —  1  pm  for  simplicity.  The  effect  of  particle  size  is  discussed  in  Section  3.3  below. 
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It  is  also  reasonable  to  assume  that  the  bacterial  particles  are  heated  with  negli¬ 
gible  internal  resistance  and  that  the  internal  bacterial  temperature  gradient  can  be 
ignored  due  to  the  small  Biot  number  (ratio  of  heat  transfer  coefficient  due  to  the 
boundary  layer  versus  the  internal  heat  conductance;  see  Remark  4  for  a  detailed 
analysis  of  this  issue.).  The  mass  of  the  particle,  mp,  is  constant.  Therefore,  we 
can  model  the  evolution  of  the  bacterial  temperature  accounting  for  the  heat  transfer 
between  the  particle  and  the  gas  with  the  following  first-order  ordinary  differential 
equation: 

=  KPV  -  T,),  (3,3) 

where  Cpp  is  the  particle  heat  capacity  at  constant  pressure  and  Afp  is  the  surface  area 
of  the  particle.  The  heat  transfer  coefficient  ( h )  is  computed  using  the  Ranz  and 
Marshall  correlation  [43,  44] 


h  =  ^2  +  0.6Rei/2Pr1/3],  (3.4) 

dp 

where  kg  is  the  thermal  conductivity  of  the  gas  and  Pr  is  the  Prandtl  number  (Pr  = 
CpUg/kg).  The  thermodynamic  properties  of  the  bacterium  are  extracted  from  Datta 
[12]- 

A  fifth  order  Runge-Kutta  method  is  used  in  numerically  integrating  Equations  3.2a, 
3.2b  and  3.3  using  the  gas  velocity  and  temperature  fields  obtained  from  the  compu¬ 
tational  model  of  Section  2.3.2.  The  particle  velocity  and  temperature  are  calculated 
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every  10-7  m  along  the  centerline  of  the  impactor.  These  solutions  are  calculated 
for  all  conditions  shown  in  Table  2.1,  including  flow-fields  where  instabilities  are  seen 
(see  Remark  5). 

Compared  to  the  maximum  gas  Mach  (Mi)  shown  in  Figure  2.9-a,  the  maximum 
particle  Mach  number  ( Mav  =  max  \vp  —  vg\ /s;  s  is  the  speed  of  sound)  does  not 
exhibit  a  clear  trend  versus  x/d  (see  Figure  3.1-a).  Mav  is  approximately  equal  for  all 
operating  conditions  up  to  x/d  «  0.8.  There  is  a  peak  that  occurs  for  Pq  =  1.00  atm 
and  P\  —  0.25  atm  at  x/d  w  1.1.  This  might  be  because  of  the  effect  of  the  instabilities 
(see  Remark  5  for  a  discussion  on  this  issue).  At  higher  x/d,  the  pressure  pairs 
(Po,  P\)  =  (1.00  atm, 0.14  atm)  and  (5.00  atm, 0.70  atm)  have  higher  particle  Mach 
numbers  similar  to  the  gas  Mach  number;  however,  these  are  at  x/d  values  that  exhibit 
instabilities.  The  minimum  centerline  particle  temperature  variation  versus  x/d  of 
Figure  3.1-b  is  similar  to  the  temperature  variation  of  the  gas  phase  in  Figure  2.9-b. 
For  equal  ratio  of  P\/Pq  the  centerline  particle  temperature  drop  is  higher  for  higher 
Po.  Therefore,  although  the  gas  flow  and  temperature  are  constant  at  constant  Pi/Po, 
the  different  absolute  values  of  Po  and  P\  affect  the  bacterial  motion  and  temperature 
in  the  impactor  system.  Figures  3.2  and  3  3  show  the  particle  and  gas  velocity  profiles 
as  well  as  the  mean  free  path  as  a  function  of  the  distance  from  the  nozzle  ( x )  for 
different  operating  conditions.  In  Figure  3.2,  x/d  and  P0  are  held  constant  at  1  and 
1  atm,  respectively,  while  Pi  assumes  two  values  0.14  atm  and  0.25  atm.  The  mean 
free  path  is  expected  to  be  lower  for  Pi  =  0.25  atm  and  can  be  seen  by  comparing 
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Figure  3.1  (a)  Maximum  centerline  Mach  number  and  (b)  minimum  centerline  tem¬ 
perature  of  the  particle,  as  a  function  of  the  normalized  distance  between  the  nozzle 
and  the  plate  ( x/d )  for  four  different  (Po,  P\)  pairs. 


Figures  3.2-c  and  3.2-d.  The  instabilities  in  the  flow  can  be  seen  in  Figure  3.2-b  where 
the  centerline  gas  velocity  assumes  a  negative  velocity  because  of  the  recirculation 
bubble.  The  particle  stops  at  the  point  where  the  gas  velocity  is  zero  (see  Remark  5 
for  more  details).  In  Figure  3.3,  x/d  and  P\ / Po  are  held  constant  at  1.4  and  0.14, 
respectively,  while  Po  assumes  two  values  1  atm  and  5  atm.  Although  the  maximum 
centerline  Mach  numbers  were  the  same  for  both  cases  (see  Figure  2.9-a),  the  gas 
velocity  drop  and  mean  free  path  are  different.  The  drop  from  maximum  gas  velocity 
to  100  m/s  occurs  over  a  length  of  150  fim  and  240  fum  for  Po  =  1  atm  and  Po  =  5  atm, 
respectively.  This  will  affect  the  forces  experienced  by  the  bacterial  particles  as  they 
go  through  the  impactor  system  (see  Section  3.3  below). 


Remark  4  The  computational  analysis  of  the  impactor  revealed  a  maximum  Biot 
number  of  0.31  (x/d  =  0.5,  P0  =  5  atm  and  Pi  =  0.7  atm.)  using  the  following 
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Figure  3.2:  (a,  b)  Centerline  particle  and  gas  velocity  as  a  function  of  the  distance 
from  the  nozzle,  (c,  d)  Mean  free  path  of  the  gas  molecules  as  a  function  of  distance 
along  the  centerline  with  Om  marking  the  exit  of  the  nozzle.  Conditions:  (a,  c) 
x/d  =  1.0,  Pq  -  1  atm  and  P\  —  0.14  atm,  (b,  d)  x/d  —  1.0,  Pq  =  1  atm  and 

Pi  =  0.25  atm. 


formula  [26] 


(3.5) 


where  kv  is  thermal  conductivity  of  the  bacterium  and  is  given  by  Datta  to  be  0.44  W/m  K 
[12].  Although  the  maximum  Biot  in  the  parametric  study  is  greater  than  0.1 ,  it  is 
still  well  less  than  1,  and  thus ,  equation  3.3  is  a  good  approximation  since  the  internal 
resistance  to  heat  transfer  is  well  less  than  the  external  one. 
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Figure  3.3:  (a,  b)  Centerline  particle  and  gas  velocity  as  a  function  of  distance  from 
the  nozzle,  (c,  d)  Moan  free  path  of  the  gas  molecules  as  a  function  of  distance  along 
the  centerline  with  0  m  marking  the  exit  of  the  nozzle.  Conditions:  (a,  c)  x/d  =  1.4, 
P0  =  1  atm  and  Pi  —  0.14  atm,  (b,  d)  x/d  =  1.4,  P0  =  5  atm  and  Pi  =  0.70  atm. 


Remark  5  In  the  unstable  flow  conditions  ( i.e presence  of  recirculation  bubbles  in 
the  flow-field)  discussed  m  Section  2.3.2,  the  use  of  Equations  3.2a  and  3.2b  does  not 
capture  the  particle  behavior  in  the  experimental  impactor  system.  Specifically ,  in  this 
case ,  the  computations  reveal  that  the  particles  stop  before  reaching  the  collection  plate 
and  in  some  cases  before  entering  the  deceleration  tube.  This  type  of  behavior  arises 
because  of  the  negative  velocities  of  the  gas  on  the  centerline  caused  by  the  recirculation 
bubble.  Furthermore,  the  particle  motion  Equations  3.2a  and  3.2b  do  not  account  for 
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forces  such  as  gravity,  Brownian  motion ,  and  turbulent  dispersion  which  may  need 
to  be  included  to  improve  the  accuracy  of  the  particle  motion  trajectories  calculations 
at  unstable  flow  conditions.  The  effects  of  Brownian  motion  and  turbulent  dispersion 
on  the  particle  trajectories  and  accelerations  will  be  the  subject  of  future  work  To 
study  these  forces  more  detailed  simulations  are  necessary  where  a  cloud  of  particles 
is  released  from  the  centerline  and  average  properties  are  calculated . 


3.3  Bacterial  envelope  instability 


Instabilities  in  the  bacterial  membrane  arise  if  the  bacterial  particle  is  accelerated 
(or  decelerated)  in  a  perpendicular  direction  to  the  bacteria-gas  interface  (direction 
parallel  to  the  centerline  rr-direction)  and  are  referred  to  as  Rayleigh-Taylor  insta¬ 
bilities  [7,  28  .  The  relative  acceleration  creates  waves  in  the  surface  of  the  bacteria 
with  characteristic  wavelengths.  These  waves  can  grow  uncontrollably  if  not  for  the 
stabilizing  effect  of  surface  tension.  However,  once  a  critical  wavelength  is  reached, 
surface  tension  cannot  compensate  for  the  growth  of  waves  on  the  bacterial  surface 
and  the  bacterium  breaks  up.  This  phenomenon  has  been  studied  extensively  by 
Chandrasekar  [7]  and  Joseph  [28].  Specifically,  a  bacterium  is  predicted  to  break  up 
if 

dry  >  Ac  =  27T.  (3.6) 

V  Pp^C 
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where  Ac  is  the  critical  wavelength  and  is  directly  related  to  the  critical  bacterial 
acceleration  ac.  In  Equation  3.6,  a  is  the  surface  tension,  pp  is  the  particle  density 
and  dp  is  the  particle  diameter.  Other  conditions,  which  have  to  be  satisfied  for  the 
critical  acceleration  to  cause  a  break-up,  are  expressed  in  terms  of  two  dimensionless 
numbers  [28].  Specifically, 


We  >  Wec  =  12(1  +  1.077 Oh1-6) 


(3.7a) 


where 

\\rp  _  Pg^vr  ~  %-)2^p 

a 


(3.7b) 


and 


Oh  = 


Pp 

(pP<V)1/2 


(3.7c) 


where  We  is  the  Weber  number,  which  is  the  ratio  between  the  inertial  force  exerted 
on  the  bacterial  particle  and  the  particle  surface  tension  force.  Oh  refers  to  the  Ohne- 
sorge  number,  which  is  the  ratio  between  the  viscous  forces  and  the  surface  tension 
force.  vPx  and  v9x  are  the  particle  and  gas  velocity  in  the  x-direction,  respectively.  We 
can  therefore  infer  bacterial  membrane  break-up  by  computing  the  relative  velocity 
(vPx  ~~  v9x)  of  the  bacterium,  as  well  as  the  maximum  acceleration  achieved.  Since 
the  particle  grid  and  the  gas  grid  do  not  coincide,  the  gas  velocity  is  linearly  interpo¬ 
lated  to  obtain  data  at  the  specific  positions  where  the  particle  velocity  is  calculated. 
Specifically,  the  relative  velocity  is  then  numerically  differentiated  with  respect  to 
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time  to  give  the  maximum  acceleration  as  follows: 


= 


u9xx 


u9*i+ i 


ti-U 


+  1 


(tj  tx)  +  ^gxx  >  ^  G  [ti, 


(3.8a) 


max  —  max 

j 

where  i  and  j  are  the  indices  of  gas  and  particle  grid  points,  respectively,  and  tx  and 
1 7  are  the  times  at  the  grid  point  i  and  j,  respectively.  Only  grid  points  between  the 
nozzle  exit  and  the  deceleration  tube  entrance  are  considered  in  the  calculation  of  the 
maximum  acceleration  in  the  x-direction. 

Table  3.1:  Biological  cell  diameter  and  surface  tension  [61]  (not  limited  to  bacteria) 


Cell 

dp  (/im) 

a  (N/m) 

E.  coli 

1 

7.5  x  10"3 

D .  carota 

60 

4.5  x  101 

M.  hungatei 

0.44 

3.5  -  5 

C.  eugametos 

16 

3.8  x  101 

B  Emersomi 

20 

3.2  x  101 

(VJ.J+1  “  VP*j  +  l)  ~  ~  VP*j  ) 


tj  +1  tj 


(3.8b) 


Table  3.2:  Critical  shock  properties  needed  to  induce  bacterial  break-up  for  different 
biological  cells.  The  maximum  Weber  number  is  calculated  for  1  particles  (see 
Remark  6) 


Cell 

ac  (  m/s2) 

Oh 

Wec 

^mai 

E.  coli 

3.0  x  10s- 

3.7  x  10-1 

1.4  x  101 

2.9  x  101 

D.  carota 

4.9  x  108 

6.1  x  10-4 

1.2  x  101 

3.0  x  10'1 

M.  hungatei 

0.7  -  1.0  x  1012 

2.1  -  2.5  x  lO”2 

1.2  x  101 

2.0  -  2.8  x  10'2 

C.  eugametos 

5.9  x  109 

1.3  x  10'3 

1.2  x  101 

9.4  x  10-2 

B.  Emersonii 

3.2  x  1010 

1.3  x  10“3 

1.2  x  101 

1.4  x  10'1 
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Figure  3.4  shows  the  acceleration  of  the  bacterial  particles  relative  to  the  gas 
stream.  The  bacterial  envelope  becomes  unstable  at  a  critical  acceleration  (ac  = 
(4n2a)/(ppdp2))  which  depends  on  the  surface  tension,  density  and  diameter  of  the 
bacterium.  The  values  of  ac  for  different  bacteria  are  shown  in  Table  3.1  along  with 


the  Weber  and  Ohnesorg  numbers.  If  the  maximum  calculated  Weber  number  in  the 
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Figure  3.4:  Bacterial  particle  acceleration  (relative  to  gas  flow)  as  a  function  of  dis¬ 
tance  from  the  nozzle  outlet,  (a)  The  operating  conditions  are  constant  for  all  three 
curves  at  P0  =  1  atm  and  Pi  =  0.14  atm.  The  xjd  varies  for  each  case  from  0.5  to 
1.4.  (b)  The  x/d  is  held  constant  at  1.4  as  well  as  Pi/Pq  at  0.14.  The  acceleration  is 
higher  for  the  case  where  Pq  =  5  atm. 


impactor  system  ( Wemax ),  exceeds  the  critical  Weber  number  Wec ,  then  conditions 
are  satisfied  to  achieve  bacterial  break-up  at  the  critical  accelerations  (ac).  Spores 
(i.e.,  properties  similar  to  M.  hungatei )  have  a  higher  surface  tension  (lower  Wemax 
and  higher  ac)  than  vegetative  cells  (i.e.,  E.  coli)  and  hence,  they  are  more  difficult 
to  damage  mechanically.  Table  3.1  indicates  that  all  of  the  cells  except  E.  coll ,  does 
not  satisfy  the  critical  Weber  number  requirement;  basically,  the  cells  do  not  achieve 
the  energies  needed  for  break-up.  The  computational  results  shown  in  Figure  3.5 
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Figure  3.5:  Maximum  centerline  acceleration  of  the  particle  as  a  function  of  nor¬ 
malized  distance  between  the  nozzle  and  the  plate  (x/d)  for  four  different  (P0,  P\) 
pairs. 

indicate  that  the  maximum  accelerations  achieved  in  the  developed  experimental 
impactor  system  are  on  the  order  of  109  —  1010  m/s2  including  data  for  a  wide  range 
of  x/d  values.  Comparing  this  data  with  the  critical  accelerations  in  Table  3  1,  we  can 
predict  the  possibility  of  break-up  of  certain  bacterial  particles  in  the  experimental 
impactor  system.  For  example,  E.  coli ,  a  vegetative  gram-negative  bacterium,  requires 
an  acceleration  of  3.0  x  10s  m/s2  and  we,  therefore,  predict  that  the  experimental 
impactor  is  capable  of  neutralizing  this  bioaerosol.  On  the  other  hand  M.  hugatei, 
an  archaea,  (see  Chapter  5  for  B .  atropheus  data)  requires  an  acceleration  of  10  x 
1012  m/s2  which  can  not  be  achieved  by  the  developed  experimental  impactor  system 

Remark  6  Particle  size  affects  the  relative  velocity  between  the  bacterial  particle  and 
the  gas,  which  leads  to  a  change  m  the  Weber  number.  In  this  work,  the  Weber  number 
of  a  0.5  fim  bacterial  particle  was  compared  to  that  of  a  1.0  ptm  bacterial  particle  at 
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Figure  3.6:  Maximum  centerline  acceleration  of  the  particle  as  a  function  of  the 
particle  size  in  impactor  system  operated  at  P0  —  1  atm  and  Pi  =  0.14  atm  for  two 
values  of  x/d. 

Pq  =  1  atm,  Pi  =  0.14  atm  and  x/d  —  1.4.  The  difference  in  value  was  ^  5%, 
which  is  not  significant  considering  that  the  maximum  Weber  number  of  M.  hungatei 
(dp  =  0.44  pirn)  is  much  smaller  than  the  critical  one  (see  Table  3.1).  The  amax 
values  also  change  with  different  particle  sizes  as  shown  in  Figure  3.6 .  However, 
these  changes  are  not  significant  since  they  are  of  the  same  order  of  magnitude. 

3.4  Conclusions 

The  computational  fluid  dynamics  model  of  the  experimental  impactor  system  of 
Figure  2.1  developed  in  Chapter  2  was  used  to  predict  the  evolution  of  both  the 
gas  and  particle  velocity  and  temperature  profiles  in  the  impactor  system.  For  this 
chapter  these  profiles  were  used  to  compute  the  critical  acceleration  needed  to  break¬ 
up  different  bacterial  aerosols  using  the  Rayleigh-Taylor  instability  equations.  Our 
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models  predict  that  the  bacterial  accelerations  achieved  in  the  impactor  system  of 
Figure  2.1  are  on  the  order  of  109  —  1010  m/s2  for  various  conditions  of  x/d,  Pq  and 
Pi.  Certain  bacterial  aerosols  which  have  critical  accelerations  less  than  the  ones 
computed  in  the  impactor  system  are  predicted  to  break-up. 
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Chapter  4 

Effect  of  Aerodynamic  Shock  on  E. 
coli 

4.1  Introduction 

The  computational  model  developed  in  Chapter  3  predicts  the  break-up  of  E  colt 
vegetative  cells  in  the  impactor  system.  In  order  to  validate  our  computational  model, 
E.  colt  cells  are  aerosolized  and  passed  through  the  experimental  impactor  system. 
Different  operating  conditions  are  used  in  order  to  validate  both  break-up  and  no 
break-up  predictions  in  the  system.  Measurements  of  viability  are  carried  out  on  dye 
labeled  cells  using  fluorescence  flow  cytometry  [49] . 
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4.2  Experimental  materials  and  methods 


4.2.1  Experimental  setup 

The  impactor  system  discussed  in  Chapter  2  was  used  in  the  experimental  setup  and 
is  shown  in  Fig.  4.1.  The  bacterial  suspension  (see  Section  4.2.2)  is  aerosolized  using 
a  capillary  nebulizer  (TR-30-A1,  Meinhard  Glass  Products)  at  a  gas  (A^)  critical 
flow  rate  of  0.2  L/min  (see  Fig.  4.1;  parts  1  and  2).  The  impactor  (see  Fig.  4.1; 
parts  3  thru  8)  operates  at  a  critical  flow'  rate  of  2.3  L/min  at  shock  conditions.  As 
shown  in  Fig.  4.1,  a  high  efficiency  particulate  air  (HEPA)  filter  is  installed  between 
the  nebulizer  (part  2)  and  the  impactor  system  (part  3)  to  allow  particle-free  air 
to  compensate  for  the  higher  impactor  flow  rate,  thereby  maintaining  the  upstream 
impactor  pressure  at  1  atm. 

Nebulization 

The  phosphate  buffered  saline  (PBS)  (BP2438-4,  Fisher  Scientific)  suspension  con¬ 
taining  the  E.  coli  cells  (see  Section  4.2.2)  is  fed  via  a  syringe  pump  (KDS410,  Kd 
Scientific)  at  a  liquid  flow  rate  Q  =  0.05  mL/min  for  20  min.  The  liquid  is  fed  at  a 
specified  bacterial  concentration  (see  Section  4.2.2)  determined  by  the  optical  density 
measurement  at  600nm  (OD6 oo)-  Nitrogen  is  supplied  to  the  nebulizer  at  a  pressure 
Pn  via  a  pressurized  tank.  When  the  liquid  is  mixed  with  the  air  in  the  annular  re¬ 
gion  of  the  nebulizer  tip,  droplets  containing  bacterial  cells  are  formed.  The  bacterial 
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aerosol  concentration  and  droplet  size  can  be  adjusted  by  changing:  Q,  OD6 oo  and 
Pn,  Increasing  Q  and  decreasing  Pn  increases  the  average  droplet  diameter  (and  vol¬ 
ume).  Increasing  ODsoo  increases  the  number  of  cells  per  volume  of  liquid  nebulized. 
Collectively,  increasing  Q  and  OD60o  and  decreasing  Pn,  increases  the  number  of  cells 
per  droplet.  Q  and  OL)$ 0o  also  increase  the  number  of  cells  input  into  the  system  in  a 
given  time.  In  all  the  experiments,  the  nebulizer  was  operated  at  Q  =  0.05  mL/min, 
ODqoq  =  0.1,  and  Pn  =  2  atm.  The  cell  concentration  at  OD^oo  =  0.1,  as  measured 
using  fluorescence  flow-cytometry  (see  Section  4.2.3),  was  2.66±0.28  x  107  cells/mL. 
The  expected  average  nebulized  droplet  size  is  10  /xm,  which  translates  to  about  1 
cell  in  every  1000  droplets  assuming  that  no  aggregates  are  formed  in  the  suspension. 

In  previous  liquid  suspension  shock  studies,  it  has  been  shown  that  an  increase 
in  cell  concentration  results  in  an  increase  in  bacterial  viability  [40].  The  controlled 
bacterial  nebulization  process  produces  individual  cells  suspended  in  air  that  can  pass 
isolated  through  the  shock,  thereby  avoiding  the  concentration  effect  seen  in  liquid 
suspension  studies. 

Impactor  system 

The  impactor  system  consists  of  a  plate  designed  for  the  collection  of  the  bacterial 
aerosol  (see  Fig.  4.1;  parts  4  and  5)  which  is  placed  perpendicular  to  the  gas  flow 
emerging  from  a  converging  nozzle  (see  Fig.  4  1;  part  3).  Isentropic  flow  theory  of  an 
ideal  gas  predicts  a  critical  downstream  (Pi)  to  upstream  (P0)  pressure  ratio  below 
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which  the  flow  at  the  exit  of  the  nozzle  is  sonic  [38,  48],  according  to  the  following 
equation: 

Xcrit  =  (Pl/P0)crit  =  [2/(7  +  (4.1) 

where  7  is  the  heat  capacity  ratio  (7  —  1.4;  Xcrit  ~  0.53  for  N2).  The  generated 
impinging  flow  on  a  plate  from  a  converging  nozzle  operating  under  sonic  conditions 
results  in  the  creation  of  a  standoff  shock  whose  properties  can  be  changed  by  varying 
impactor  geometry  (distance  between  nozzle  and  plate  x  and  diameter  of  nozzle  d) 
and  operating  conditions  [2,  13,  14,  29,  42,  51].  In  traditional  impactors  [21,  22],  the 
aerosol  is  collected  by  impaction  on  a  flat  surface.  Stewart  et  al.  (1995)  showed  that 
bacteria  collected  by  impaction  are  injured  [53];  therefore,  assessing  damage  induced 
by  the  standoff  shock  is  impossible  using  current  impactor  designs.  To  avoid  loss 
of  viability  by  impaction,  the  flat  surface  of  the  developed  impactor  system  has  a 
central  opening  (see  Fig.  4.1;  part  4)  that  acts  as  a  virtual  surface  through  which  the 
bacteria  enter  a  stagnant  gas  deceleration  tube  (see  Fig.  4.1;  parts  4  and  5).  Biswas 
and  Flagan  (1988)  have  designed  a  similar  tube  to  prevent  particle  reentrainement  [4]. 
Computational  results  (see  Chapters  2  and  3)  indicate  that  the  bacteria  are  collected 
at  low  impact  velocities  (<  10  m/s)  in  the  deceleration  tube  and  consequently  will  not 
be  injured  due  to  impaction.  The  distance  between  the  opening  of  the  deceleration 
tube  and  the  nozzle  outlet  (x)  can  be  adjusted  with  spacers  of  different  thickness  (see 
Fig.  4.1;  part  6).  In  the  current  experimental  study,  the  spacer  is  1.4  mm  thick.  A 
valve  attached  to  a  purge  stream  downstream  of  the  impactor,  controls  whereas 
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Figure  4.1:  Schematic  of  the  experimental  system.  (1)  Meinhard  nebulizer  with  a 
concentric  nozzle  mixing  a  liquid  bacterial  suspension  feed  and  a  dispersion  gas  ( N2 ), 
(2)  Nebulization  chamber  to  collect  large  droplets,  (3)  Converging  nozzle  with  an  exit 
diameter  of  d  —  0.5  mm,  (4)  and  (5)  a  flat  surface  with  a  0.5  mm  hole  combined  with 
a  screw  that  accommodates  the  collection  substrate  make  up  the  deceleration  tube, 

(6)  1.4  mm  spacer  fixes  the  spacing  x  between  the  impactor  and  nozzle  to  0.6  mm, 

(7)  support  holding  the  deceleration  tube  and  flat  surface  at  fixed  distance  from  the 
nozzle,  and  (8)  exit  chamber. 
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the  upstream  is  left  at  ambient  pressure  (1  atin). 


The  deceleration  tube  is  filled  with  600  ju L  of  PBS,  which  keeps  the  liquid  level 
to  about  4  mm  away  from  the  entrance  of  the  deceleration  tube.  The  PBS  provides  a 
hydration  environment  for  the  collected  cells  and  ensures  that  the  collected  cells  do  not 
lose  viability  because  of  drying  (this  has  been  experimentally  verified).  The  impactor 
system  with  x/d  =  1.2  was  operated  for  20  min  at  three  pressure  ratios:  x  =  0.11, 0.50 
and  0.98.  As  a  control  to  test  any  low-pressure  effects,  600  /jl L  of  bacterial  suspension 
(sec  Section  4  2.2)  at  O/9600  =  0.01  was  placed  in  the  deceleration  tube  and  the 
impactor  was  run  with  particle- free  air  for  20  min  at  the  three  operating  pressures. 
We  experimentally  found  that  cells  retain  their  viability  under  low  pressure  conditions 
when  compared  with  the  starting  bacterial  suspension. 

4.2.2  Bacterial  cell  culture  preparation 

E.  coh  K-12  wild-type  strain  (MG1655)  was  used  as  a  standard  bacterial  aerosol  to 
study  the  effect  of  aerodynamic  shocks  on  vegetative  bacterial  cells  in  the  impactor 
system.  All  cell  cultures  were  initiated  from  a  stock  E.  coh  suspension  (80%  glycerol 
solution)  stored  at  —20  °C.  5  mL  of  Luria  Bertani  (LB)  broth  (Fisher)  was  inoculated 
with  3  /XL  of  stock  suspension  in  a  test  tube.  The  culture  was  incubated  in  a  shaker 
(Excella  E24,  New  Brunswick  Scientific)  at  37  °C  at  200  rpm  for  approximately  16  hrs. 
Subsequently,  2  mL  of  culture  was  centrifuged  at  4400  rpm  for  5  min  (Centrifuge  5702, 
Eppendorf)  and  washed  (3X)  with  4  mL  PBS.  After  washing,  the  optical  density  of 
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the  suspension  was  measured  at  a  wavelength  of  600  nm  (OZ)600)  using  a  visible  light 
spectrophotometer  (Genesis  20,  Thermo  Scientific).  The  cell  suspension  was  diluted 
to  ODq oo  =  0.1  to  be  used  in  the  experiments  described  in  Section  4.2.1.  Another 
dilution  at  ODb0 0  —  0.01  was  prepared  for  suspension  flow  cytometry  measurements, 
but  not  nebulized. 

Positive  and  negative  controls  were  prepared  and  carried  out  before  the  start 
of  every  experiment.  The  negative  control  is  the  cell  suspension  in  PBS  diluted  to 
OD(y oo  =  0.01  as  described  above.  The  negative  control  serves  as  the  starting  point  for 
viability  measurements.  It  is  the  measurement  with  most  viable  cells.  The  positive 
control  was  prepared  by  centrifuging  2  mL  of  the  cell  culture  at  4400  rpm  and  washing 
with  4  mL  70%  ethanol.  The  ethanol-washed  cells  were  incubated  for  15  min  at  room 
temperature.  The  suspension  was  then  centrifuged  at  4400  rpm  and  washed  (2x) 
with  4  mL  PBS.  The  positive  control  was  then  diluted  to  OD60o  =  0.01.  The  effect 
of  aerodynamic  shock  was  then  compared  to  a  standard  ethanol  disinfection  method 
(positive  control)  used  in  liquid  suspension. 

4.2.3  Bacterial  viability  analysis 

Bacterial  viability  is  determined  using  LIVE/DEAD®£?acLightTAf  (cat.  number 
L7012;  Life  Technologies  Corporation)  3,  32].  The  assay  contains  two  components: 
(A)  3.34  mM  STYO  9,  a  green-fluorescent  nucleic  acid  stain,  in  dimethyl  sulfoxide 
(DMSO)  and  (B)  20  mM  propidium  iodide  (PI),  a  red-fluorescent  nucleic  stain  in 
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DMSO.  STYO  9  penetrates  both  intact  and  damaged  membranes,  while  PI  oidy  pen¬ 
etrates  damaged  membranes.  A  stock  dye  solution  is  prepared  at  the  start  of  every 
experiment  by  mixing  5  /xL  of  each  of  components  A  and  B  with  10  /xL  PBS.  All 
controls  and  experimental  samples  mentioned  in  Sections  4.2.1  and  4.2.2  were  stained 
with  the  stock  dye  solution  to  a  final  concentration  of  8.35  /xM  and  50  /xM  of  STYO 
9  and  PI,  respectively. 

The  stained  samples  were  analyzed  using  a  fluorescence  flow  cytometer  (FACScan™, 
BD  Biosciences).  The  cytometer  was  operated  at  12  /xL/min  for  20  s  with  a  488  run 
excitation  laser.  Red  fluorescence  was  measured  with  a  high  pass  filter  at  630  run 
(FLH-2)  and  green  fluorescence  was  measured  at  a  bandpass  of  520  nm  (FLH-1). 
The  forward  scatter  (FSC)  and  side  scatter  (SSC)  were  also  measured  for  the  parti¬ 
cles  passing  the  laser.  Photomultiplier  tubes  were  used  to  amplify  the  optical  signal 
for  FLH-1,  FLH-2  and  SSC  at  850,  800  and  375  volts,  respectively.  A  photodiode  was 
used  to  detect  the  FSC  signal. 

4.2.4  Mass  balance  analysis 

A  0.2  fi m  black  polycarbonate  filter  (Milhpore)  was  placed  in  an  in-line  25  mm  filter 
holder  inserted  after  the  nebulization  chamber  and  before  the  HEPA  filter  of  the 
experimental  setup  shown  in  Fig.  4.1.  After  running  the  experiment  for  20  min, 
the  filter  was  placed  in  a  Millipore  filter  holder  and  1  mL  of  STYO  9  stock  was 
passed  through  the  filter  staining  all  cells.  After  staining,  the  filter  was  placed  on 
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a  microscope  slide  and  5  fiL  of  mounting  oil  (cat.  number  L7012;  Life  Technologies 
Corporation)  was  added  to  the  surface  of  the  filter,  and  covered  with  a  glass  slip. 
Images  of  the  filter  were  taken  with  a  fluorescence  microscope  at  an  objective  lens 
magnification  of  40 x  (Zeiss).  Each  image  represents  a  filter  area  of  3.8  x  104  /jm2. 
Images  of  10  randomly  selected  areas  on  the  filter  were  taken  to  get  an  average  cell 
count  per  image  area,  which  was  then  multiplied  by  the  total  area  of  the  filter  to 
get  the  total  number  of  cells  collected  on  the  in-line  filter.  A  total  of  7  independent 


(a)  (b) 


Figure  4.2:  Fluorescence  microscopy  images  taken  at  40X  objective  lens  with  an  area 
of  3.8  x  104  /mi2.  The  cells  are  stained  with  STYO  9  dye  which  penetrates  both  live 
and  dead  cells.  The  images  are  taken  for  counting  purposes  and  not  for  determination 
of  cell  viability,  (a)  image  of  E.  coli  from  cell  suspension  in  PBS  at  OD60 o  =  0.01 
on  a  polycarbonate  filter,  (b)  image  of  nebulized  E.  coli  collected  on  an  inline  filter 
between  the  nebulizer  and  the  entrance  of  the  impactor  system  operated  at  a  flowrate 
of  2.3  L/min  for  20  min. 


experiments  with  different  starting  cell  cultures  were  carried  out  The  cells  collected 
on  the  filter  were  assumed  to  be  equal  to  the  cells  entering  the  impactor  system. 
Fig.  4.2-b  shows  a  sample  image  of  bacterial  aerosol  collected  on  a  filter  before  the 
impactor  entrance  and  stained  with  STYO  9.  Fig.  4.2-a  is  an  image  from  OD600  — 
0.01  PBS  E.  coli  suspension  and  is  shown  for  comparison.  We  can  see  that  both 


57 


Fig.  4.2-a  and  b  have  a  single  layer  of  cells  and  can  be  therefore  easily  be  counted. 


4.3  Computational  results 

4.3.1  Gas  flow  field  and  particle  dynamics 

For  the  two-phase  aerosol  flow  (air-bacteria)  present  in  the  impactor  system,  the  equa¬ 
tions  for  gas  and  particle  dynamics  are  one-way  coupled  at  low  aerosol  concentrations 
18] ;  the  existence  of  aerosol  has  little  effect  on  the  gas  while  the  aerosol  flow  is  de¬ 
termined  by  the  gas  velocity.  Therefore,  the  gas  dynamics  are  solved  independently 
using  the  Reynolds  time-averaged  Navier-Stokes  equations  shown  in  Cartesian  tensor 
form: 


(4.2) 


and 


(4.3) 


where  p  is  the  density,  v  is  the  velocity,  p  is  the  static  pressure,  5  is  the  Kronecker 
delta,  p.  is  the  viscosity,  and  i ,  j ,  and  l  are  the  cartesian  indices.  To  close  the  form  in 
Eq.  4.3,  the  Reynolds  stress  term  —pv\v^  is  calculated  using  the  Boussinesq  hypothesis 
[24] .  A  standard  k  —  e  turbulence  model  is  used  in  the  gas  flow  equations. 
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Figure  4.3:  Computational  domain  of  the  impactor.  The  dashed  lines  indicate  differ¬ 
ent  virtual  surfaces.  These  surfaces  are  used  to  calculate  the  particle  distribution  at 
various  cross-sections  of  the  impactor  shown  as  S-l,  S-2,  S-3  and  S-4. 


The  computational  model  details  used  to  solve  both  gas  and  particle  dynamics 
(flow  and  temperature)  for  the  computational  domain  shown  in  Figure  4.3  arc  dis¬ 
cussed  in  Chapter  2  and  3  [51].  In  this  Chapter,  temperature  effects  are  not  analyzed. 
After  obtaining  the  gas  flow  field,  the  particle  trajectory  is  solved  using  Newton’s  sec¬ 
ond  law: 

i CDpgAp{vg  -  vp)\vg  -  Vp | ,  (4.4a) 


P’ 


(4.4b) 


where  vp  and  vg  are  the  velocity  of  the  bacterial  particle  and  gas,  respectively,  mp 
is  the  mass  of  the  particle,  pg  is  the  density  of  the  gas,  Ap  is  the  projected  area 
of  the  bacterium  on  the  plane  perpendicular  to  the  flow  direction,  and  Co  is  the 
drag  coefficient.  The  thermophoretic,  gravitational,  and  Basset  history  term  arc 
neglected  in  Eq.  4.4a.  The  drag  force  is  assumed  to  be  the  major  force  acting 
on  the  bacterial  particle  which  is  typical  for  impactor  systems.  The  bacterium  is 
simulated  as  a  1  gm  sphere  with  a  density  pv  of  1  g/mL  which  is  representative 
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of  values  reported  in  literature  [33,  59].  All  calculations  are  made  to  mimic  the 
experimental  conditions  in  the  study  with  the  x/d  =  1.2  and  both  shock  (\  =  P\/Pq  = 
0.11,0.50)  and  non-shock  (x  —  0.98)  conditions.  The  detailed  conditions  are  shown 
in  Section  4.3.3.  The  computational  results  provide  predictions  for  bacterial  envelope 
instability  (see  Section  4.3.2)  and  the  collection  efficiency  of  the  impactor  system  at 
different  operating  conditions  (see  Section  4.3.3). 


4.3.2  Bacterial  envelope  instability 

The  acceleration  (or  deceleration)  of  a  bacterium  relative  to  the  gas  phase  creates 
instabilities  in  the  bacterial  envelope  and  break-up  of  the  particle  [7,  28.  51].  Specif¬ 
ically,  waves  are  created  on  the  bacterial  surface  perpendicular  to  the  direction  of 
acceleration  (x-direction)  which  are  stabilized  by  the  bacterial  surface  tension.  Un¬ 
der  the  assumption  of  viscous  potential  flow  around  a  spherical  cell,  an  estimate  of 
the  critical  acceleration  needed  to  overcome  the  stabilizing  effect  of  surface  tension 
can  be  computed  using  the  derivation  in  [28]  (see  also  Chapter  10  of  [7])  as  follows: 


ac  —  47t2 


Ppdp 


(4.5) 


where  a  is  the  surface  tension,  pp  is  the  particle  density  and  dp  is  the  particle  diameter. 
Other  sufficient  conditions,  which  have  to  be  satisfied  for  the  critical  acceleration  to 
cause  a  break-up,  are  expressed  in  terms  of  two  dimensionless  numbers  [28].  Specifi- 
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cally, 


We  >  Wec  =  12(1  +  1.077Oh16) 


(4.6a) 


where 

_  Pg(VPx  ~  %  )2(^P 
a 


(4.6b) 


and 


Oh  = 


Pp 

(, PvdvG )1/2 


(4.6c) 


where  We  is  the  Weber  number,  which  is  the  ratio  between  the  inertial  force  exerted 
on  the  bacterial  particle  and  the  particle  surface  tension  force.  Oh  refers  to  the  Ohnc- 
sorge  number,  which  is  the  ratio  between  the  viscous  forces  and  the  surface  tension 
force.  vPt  and  v9i  are  the  particle  and  gas  velocity  in  the  z-direction  (i.e.  direction 
of  flow  along  centerline),  respectively.  We  can  therefore  infer  bacterial  membrane 
break-up  by  computing  the  relative  velocity  (vPi  —  v9i)  of  the  bacterium,  its  well  as 
the  maximum  acceleration  achieved  (see  Chapter  3).  The  criteria  that  both  decel¬ 
eration  and  Weber  number  have  to  be  above  the  critical  values  are  sufficient  for  the 
bacterial  cell  to  break-up  as  they  pass  through  the  shock.  However,  if  only  the  We¬ 
ber  number  is  above  its  critical  value  the  cells  may  experience  membrane  vibrations, 
boundary  layer  stripping,  or  even  break-up  [28].  Furthermore,  the  model  does  not 
account  for  loss  of  viability  without  break-up  which  is  a  more  stringent  requirement 
than  other  mechanisms. 

Bacterial  cells  with  a  lower  surface  tension  and  larger  diameters  require  lower 
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Table  4.1:  Biological  cell  properties  [61]  (not  limited  to  bacteria)  and  critical  shock 


properties  needed  to  induce  bacterial  break-up. 


Cell 

dp  (/an) 

a  (N/m) 

ac  (m/s2) 

Oh 

Wec 

E.  coli 

1 

7.5  x  10~3 

3.0  x  108 

3.7  x  10"1 

14 

D.  carota 

60 

4.5  x  101 

4.9  x  108 

6.1  x  10"4 

12 

M.  hungatei 

0.44 

3.5  —  5 

0.7-  1.0  x  1012 

2.1  -  2.5  x  10“2 

12 

C.  eugametos 

16 

3.8  x  101 

5.9  x  109 

1.3  x  10-3 

12 

B.  Emersomi 

20 

3.2  x  101 

3.2  x  1010 

1.3  x  10~3 

12 

accelerations  to  break-up.  It  is  important  to  note  that  the  surface  tension  values 
reported  in  Table  4.1  are  derived  from  the  turgor  pressure  [61].  However,  the  surface 
tension  can  vary  with  deformation  generated  by  fluid  flow  over  the  cell  membrane. 
Such  a  variation  of  the  cell  surface  tension  under  fluid  flow  may  limit  the  accuracy 
of  the  critical  acceleration  predictions  of  Table  4.1.  In  such  a  case,  more  accurate 
models  of  the  variation  of  cell  surface  tension  under  stress  generated  by  the  fluid  flow 
(similar  to  the  works  of  [30,  31,  52,  60])  should  be  utilized  and  combined  with  the  fluid 
flow  and  particle  motion  models  developed  here  to  make  accurate  critical  acceleration 
predictions.  We  also  note  that  previous  work  [31]  has  demonstrated  for  the  first  time 
the  use  of  inertial  forces  (expressed  in  terms  of  the  Weber  number  and  acceleration)  in 
achieving  cell  break-up  in  the  context  of  a  high-pressure  homogenization  system.  As 
shown  in  Table  4.1,  vegetative  cells  represented  by  E.  coli  require  critical  accelerations 
(ac)  of  3.0  x  108  m/s2.  On  the  other  hand,  spores  represented  by  M.  hungatei,  an 
archea,  require  accelerations  of  about  four  orders  of  magnitude  more  than  vegetative 
cells.  In  our  previous  parametric  study  of  the  impactor  system  [51]  where  Pi,  Po 
and  x/d  were  varied,  it  was  shown  that  we  can  achieve  accelerations  on  the  order  of 
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109  —  1010  in/s2  in  our  current  system  operating  under  shock  conditions.  Furthermore, 
it  is  possible  to  achieve  accelerations  below  10b  m/s2  by  operating  the  impactor  at 
non-shock  conditions. 

A  shock  is  represented  as  a  sharp  discontinuity  in  fluid  properties  (velocity,  tem¬ 
perature,  pressure)  and  can  be  seen  in  Figure  4.4-a.  Due  to  its  inertia,  a  bacterial 
particle  is  not  able  to  follow  the  sharp  drop  in  gas  velocity  and  will  therefore  expe¬ 
rience  a  relative  deceleration  of  4.95  x  109  m/s2  at  x  ~  0.11  and  2.02  x  108  m/s2  at 
X  =  0.50  (both  above  the  critical  acceleration  of  E.  coli).  However,  in  the  non-shock 
condition  with  x  =  0.98  (see  Figure  4.4-b)  the  gas  properties  do  not  experience  sharp 
changes.  In  the  incompressible  flow  the  gas  decelerates  because  of  the  larger  area 
due  to  the  expansion  after  exiting  the  nozzle.  The  bacterial  particle  is,  therefore, 
able  to  adjust  to  the  gas  velocity  and  experience  a  much  lower  relative  deceleration 
of  3.89  x  106  m/s2.  The  deceleration  of  the  bacterial  particle  in  the  non-shock  case 
is  not  large  enough  (<  3.0  x  108  m/s2,  see  Table  4.1)  to  cause  the  break-up  of  the 
cell  envelope.  Figures  4.4-c  and  d  show  the  acceleration  of  the  bacterial  particle  as  a 
function  of  the  distance  from  the  exit  of  the  nozzle.  In  the  shock  case  with  x  =  0.11, 
the  peak  acceleration  lasts  for  a  short  distance  and  is  3  orders  of  magnitude  larger 
than  the  non-shock  case.  We  therefore  operated  the  impactor  system  at  both  the 
shock  (x  =  0.11,0.50)  and  non-shock  (x  =  0.98)  conditions  to  obtain,  in  addition  to 
the  computational  data,  experimental  bacterial  viability  data  (see  Section  4.4). 

Remark  7  At  x  ~  0.98,  the  acceleration  of  the  bacterial  particle  shown  in  Fig- 
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Figure  4.4:  (a,  b)  Centerline  particle  and  gas  velocity  and  (c  d)  Bacterial  particle 
acceleration  (relative  to  gas  flow)  as  a  function  of  the  distance  from  the  nozzle.  The 
entrance  of  the  collection  tube  is  at  a  distance  6  x  10-4  m  from  the  nozzle  exit. 
Impactor  operating  conditions:  (a,c)  X  —  0.11  and  (b,d)  X  =  0.98. 

ure  4-4‘d  experiences  a  sharp  change  at  a  position  greater  than  0.55  mm  from  the 
exit  of  the  nozzle.  The  change  is  because  of  the  stagnant  gas  at  the  entrance  of  the 
deceleration  tube ,  especially  since  Brownian  motion  is  neglected  in  these  calculations. 
The  maximum  acceleration  is  calculated  at  a  position  approximately  0.4  mm  from  the 
nozzle . 
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Remark  8  The  particle  accelerations  were  calculated  by  tracking  a  bacterium  released 
on  the  centerline  of  the  impactor,  where  the  aerodynamic  shock  has  the  highest  Mach 
number.  Therefore ,  the  calculated  bacterial  accelerations  represent  a  maximum  in  a 
distribution  if  the  particle  accelerations  were  to  be  calculated  at  all  points.  These 
off- centerline  accelerations  can  represent  values  lower  than  the  critical  value  of  3.0  x 
10s  m/s 2  required  to  break-up  the  E.  coh  cells.  In  the  experiments,  the  E.  coh  cells 
will  experience  the  distribution  of  accelerations.  Therefore,  we  do  not  expect  all  cells 
to  lose  viability ,  especially  since  they  pass  the  shock  at  different  points  as  explained  in 
Section  f.3.3. 
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Figure  4.5:  (a,b)  Turbulent  kinetic  energy  (m2/s2)  contours  of  the  gas  phase  in  the 
converging  region  of  the  nozzle.  A  region  of  high  turbulence  can  be  seen  in  the 
converging  part.  (c,d)  Mach  contours  at  the  exit  of  the  nozzle.  Conditions:  (a,c) 
shock  with  x  =  0.11  and  (b,d)  non-shock  with  x  =  0.98.  Note  that  the  scales  are 
different  for  the  different  conditions. 
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4.3.3  Bacterial  collection  efficiency  in  the  impactor 

The  collection  efficiency  (77)  is  defined  as  the  number  of  particles  that  pass  through 
the  collection  tube  entrance  versus  the  number  of  particles  that  enter  the  impactor 
inlet.  The  bacterial  collection  efficiency  (77)  for  all  three  impactor  operating  condi¬ 
tions  (x  =  0.11,0.50,0.98)  are  evaluated  to  determine  the  overall  bacterial  aerosol 
neutralization  rate  in  Section  4.4.  For  a  highly  turbulent  flow  (see  Fig.  4.5-a,  b)  and 
low-density,  small-size  aerosols,  particle  dispersion  due  to  gas  velocity  fluctuations 
could  be  significant  and  should  be  accounted  for  in  the  calculation  of  collection  effi¬ 
ciency.  The  gas  dynamics,  described  in  Section  4.3.1,  is  based  on  a  turbulent  model 
and  only  provides  a  Reynolds-averaged  velocity  while  the  velocity  fluctuation  terms 
are  not  explicitly  described.  The  random  walk  model  is  employed  in  this  work  to 
study  the  effect  of  velocity  fluctuation  on  bacterial  aerosol  with  reasonable  compu¬ 
tational  load.  Stochastic  tracking  of  the  particles  is  used  in  order  to  account  for  the 
turbulence  dispersion  effects  of  the  gas  flow.  In  Section  4.3.2,  the  average  values  of 
the  gas  velocity  were  used  to  evaluate  the  maximum  relative  acceleration  of  the  par¬ 
ticles  on  the  centerline.  In  the  stochastic  tracking  model,  however,  the  gas  velocity 
is  represented  as  an  average  plus  a  fluctuating  term  when  calculating  the  particle 
velocity,  vp  in  Equation  4.4a.  The  fluctuating  velocity  is  kept  constant  for  time  in¬ 
tervals  corresponding  to  the  lifetime  of  the  eddies.  The  following  equations  describe 
the  stochastic  behavior: 

Vg~v  +  vf  (4.7a) 
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(4.7b) 


v  =  c\4/2  =  cvW3 


ti  =  0.15  k/e 

(4.7c) 

t(  =  -tilogr 

(4.7d) 

where  ^  is  the  instantaneous  gas  velocity,  F  is  the  average  velocity,  and  t/  is  the 
velocity  fluctuation  term,  xt  is  a  function  of  a  normally  distributed  random  number, 
and  i/'2  is  the  local  mean  square  of  the  fluctuations,  v '2  is  calculated  based  on  the 
assumption  of  isotropy  and  is  given  by  vn  =  2/c/3,  with  k  being  the  turbulence 
kinetic  energy  in  the  k  —  e  model.  Finally,  ti  is  the  lagrangian  time  interval,  e  is 
the  turbulence  dissipation  rate,  r€  is  the  lifetime  of  the  eddies,  and  r  is  a  uniformly 
distributed  random  number  between  0  and  1.  Fig.  4.6  shows  the  path  of  5  particles 
released  from  the  inlet  as  they  exit  the  nozzle  of  the  impactor  system  for  both  x  —  0.11 
(Fig.4.6-a)  and  \  =  0.98  (Fig.4.6-b).  In  calculating  r;,  10,000  particles  are  released 
at  the  inlet  of  the  impactor.  The  particles  arc  then  counted  as  they  pass  the  different 
virtual  surfaces  shown  in  Fig.  4.3  as  $  —  1,  s  —  2,  .s  —  3  (nozzle  exit)  and  s  —  4 
(deceleration  tube  entrance)  (see  Remark  9).  All  the  boundaries  shown  in  Fig.  4.3, 
except  pressure  inlet  and  pressure  outlet,  are  assumed  to  be  reflection  boundaries. 
The  collection  efficiency  is  defined  as  follows: 
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where  nx  is  the  number  of  particles  entering  the  impactor  (see  also  Remark  9)  and 
n'c  is  the  number  of  particles  that  enter  the  collection  tube.  nfc  is  calculated  under 
the  assumption  of  no  break-up.  77  is  used  to  determine  the  total  number  of  live  cells 
that  would  be  collected  if  the  shock  were  to  have  no  effect  (see  Section  4.4.2).  For 
different  operating  conditions,  the  collection  efficiency  (77)  is  reported  in  Table  4.2. 
The  maximum  acceleration  (amax)  and  Weber  number  {Wemax),  which  enable  the 
prediction  of  the  experimental  break-up,  are  also  reported  in  Table  4.2  The  critical 
acceleration  conditions  that  satisfy  the  requirements  for  break-up  as  discussed  in 
Section  4.3.2  are  met  for  x  —  0.11  and  0.50;  however,  the  conditions  for  the  Weber 
number  are  not.  Increasing  the  upstream  pressure  to  5  atm  will  increase  the  We  to 
15  while  keeping  the  maximum  acceleration  similar  as  shown  in  Table  4.2.  However, 
as  reported  in  Section  4.4.2.  cells  will  still  break-up  under  non-critical  Weber  number 
conditions.  The  Reynolds  number  at  the  exit  of  the  nozzle  and  the  maximum  Mach 
number  of  the  gas  are  also  reported  to  give  insight  into  the  flow  characteristics  of  the 
impactor.  The  effect  of  these  different  conditions  is  described  in  detail  in  Section  4.4. 

Table  4.2:  Critical  parameters  for  different  operating  conditions  with  x/d  =  1.2 


Po  (atm) 

x{—  P1/P0) 

Omai  (m/s2) 

H  ^max 

M max 

Re 

V 

1.0 

0.98 

1.67  x  10* 

1.7  x  10“2 

0.13 

9.9  x  102 

0.35 

1.0 

0.50 

2.02  x  10s 

3.5  x  10° 

0.96 

5.8  x  103 

0.72 

1.0 

0.11 

4.96  x  109 

3.1  x  10° 

2.36 

5.4  x  103 

0.51 

5.0 

0.11 

5.98  x  109 

1.5  x  101 

2.32 

3.2  x  104 

0.38 

Remark  9  Despite  the  reflective  boundary  conditions  employed  in  the  simulations, 
not  all  10,000  particles  released  from  the  inlet  of  the  impactor  are  accounted  for  at 
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Figure  4  6:  Representative  trajectories  of  5  particles  released  from  the  inlet  of  the 
impactor  for  both  (a)  x  —  0.11  and  (b)  x  =  0-98.  The  collection  efficiencies  as 
measured  by  the  number  of  particles  entering  the  collection  tube  are  0.51  and  0.36 
for  x  =  0.11  and  x  =  0.98,  respectively. 


the  exit  of  the  nozzle  (Fig.  4-3;  s  —  3).  For  x  =  0.11,  0.50  and  0.98  only  7,  763,  9840 
and  9500  particles  are  counted ,  respectively .  For  \  —  0.11  and  0.50,  105  integration 
steps  were  used  to  calculate  the  particle  track.  For  \  —  0.98,  2  x  10rj  and  4  x  105  inte¬ 
gration  steps  were  used:  in  the  first  case  3840  particles  were  counted  and  m  the  second 
9500.  Both  integration  steps  resulted  in  a  collection  efficiency  of  0.35.  Therefore,  in 
the  calculation  of  collection  efficiency  the  number  of  particles  exiting  the  nozzle  are 
assumed  to  be  the  total  number  of  particles  entering  the  impactor. 


Remark  10  Fig.  ^.7  shows  the  particle  distribution  at  different  surfaces  (s  —  l,  5  —  2, 
5  —  3,  and  s  —  4)  of  the  impactor  system.  The  particle  distribution  at  the  entrance  of 
the  collection  tube  is  shown  in  Fig.  4-7  9  and  h.  where  the  dotted  line  indicates  the 
edge  of  the  entrance.  In  the  x  —  0.98  case ,  the  peak  occurs  away  from  the  center  and 
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this  can  be  expected  for  x/d  <  4  according  to  Burwash  et  al.  [6].  Similar  distributions 
have  also  been  observed  by  other  researchers  [1,  \6\.  The  particle  distribution  in  the 
shock  condition  has  a  peak  at  the  center  of  the  collection  tube  entrance. 
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Figure  4.7:  Distribution  of  particles  along  different  surfaces  (s-1,  s-2,  s-3,  s-4)  labeled 
in  Fig.  4.3.  The  z-axis  is  the  distance  (r)  in  meters  from  the  centerline  and  the  y-axis 
is  the  percentage  of  total  particles  released  from  the  inlet  of  the  impactor.  The  dashed 
and  dotted  lines  indicate  the  wall  of  the  impactor  at  s-1,  s-2  and  s-3  respectively.  At 
s-4  the  dotted  line  indicates  the  entrance  of  the  collection  tube,  (left  plots)  Particles 
released  under  shock  conditions  and  (right  plots)  particles  released  under  non-shock 
conditions. 
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4.4  Cell  viability  results  and  discussion 


Quantitative  experimental  measurements  of  cell  viability  were  obtained  via  fluo¬ 
rescence  flow  cytometry  (see  Section  4.4.1)  and  fluorescence  microscopy  (see  Sec¬ 
tion  4.2.4).  Experimental  data  are  supplemented  with  computational  data  on  the 
collection  efficiency  of  the  impactor  (see  Section  4.3.3)  for  the  different  operating 
conditions.  The  goal  is  to  assess  the  fraction  of  cells  that  remain  alive  after  the  effect 
of  impactor  operation  (//): 


nr 


(4.9) 


where  nc  is  the  experimentally  measured  number  of  live  cells  including  the  break-up 
model  and  n'c  is  the  number  of  live  cells  that  enter  the  collection  tube  under  the 
assumption  of  no  break-up  (i.e.,  all  cells  are  alive  -  see  Eq.  4.11).  ft  for  the  different 
operating  conditions  will  be  compared  to  assess  the  effect  of  the  aerodynamic  shock 
on  the  bacterial  aerosol.  The  following  sections  will  detail  the  data  used  to  calculate 

ft- 


4.4.1  Fluorescence  cytometry 

The  experimental  fluorescence  cytometry  protocol  described  in  Section  4.2.3  is  used 
to  obtain  the  data  shown  in  Fig.  4.8.  The  data  shown  in  Fig.  4.8  represents  one  set  of 
experiments.  A  total  of  15  sets  of  experiments  were  conducted,  each  from  a  different 
starting  cell  culture,  to  obtain  statistically  significant  results  in  the  calculation  of 
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The  left  dot-plots  in  Fig.  4.8  represent  the  side-scatter  (y-axis)  and  forward-scatter 
(z-axis)  of  a  particle  passing  the  laser  beam  in  the  flow  cytometer.  Each  dot  represents 
a  particle’s  side  and  forward  scatter  intensities.  The  side  scatter  is  a  measure  of  the 
granularity  and  complexity  while  the  forward  scatter  is  related  to  the  surface  area 
(i.e.  size)  of  the  particle.  The  population  of  bacterial  particles,  enclosed  in  the  gate 
denoted  by  G2/G3  in  Fig.  4.8,  is  determined  by  running  a  PBS  only  solution  to  isolate 
the  background  enclosed  in  the  gate  denoted  by  Gl.  The  right  dot-plots  shows  the 
red  fluorescence  (FL2-H)  versus  green  fluorescence  (FL1-H)  intensity  for  the  gated 
(i.e.  inside  the  polygon  G2/G3)  cell  population.  The  analysis  of  the  fluorescence 
plots  provides  information  on  the  viability  of  the  cell  population.  In  Fig.  4.8-d  the 
EtOH  cell  population  can  clearly  be  seen  to  lie  in  the  top  left  half  denoted  by  D.  Live 
cells  can  be  counted  in  the  bottom  right  half  denoted  by  L.  Similar  patterns  have 
also  been  observed  in  the  literature  [3].  Furthermore,  Berncy  et  al.  [3]  indicate  the 
presence  of  an  intermediate  state  where  the  outer  membrane  in  gram-negative  cells  is 
compromised  The  intermediate  state  population  can  be  seen  as  the  population  with 
the  highest  intensity  green  fluorescence  in  the  L  gate.  For  the  purposes  of  the  current 
study,  the  intermediate  state  cell  population  is  grouped  with  the  live  population  in 
order  to  make  the  requirements  for  neutralization  more  stringent. 

The  fraction  of  live  cells  (cells  in  bottom  right  half  divided  by  total  gated  cells  or 
L/(L+D))  in  all  analyzed  samples  is  shown  in  Fig.  4.9.  The  error  bars  are  the  95% 
confidence  interval  limits  on  the  data  using  a  normal  fit.  The  data  does  not  take  into 
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(a)  (b) 


(g)  (h) 


Figure  4.8:  (left)  Side-scatter  versus  forward-scatter  plot  of  all  particles  in  the  flow 
cytometer.  The  G2  gate  indicates  the  particle  sizes  which  represent  intact  cells. 
G1  indicates  the  debris  population,  (right)  Red  fluorescence  (FLH-2)  versus  green 
fluorescence  (FLH-1)  of  gated  cells  in  the  left  side  plots.  The  L  gate  represents  the 
live  cells  and  the  D  gate  represents  the  dead  cells.  (a,b)  Cell  suspension  in  PBS, 
(c,d)  cell  suspension  in  EtOH,  (e,f)  cells  collected  after  passing  impactor  operated  at 
X  =  0.98.  (g,h)  x  =  0.50  and  (i,j)  x  =  0.11. 
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account  the  mass  balance  oil  the  system  to  indicate  the  overall  neutralization  rate. 
The  data  as  shown  in  Fig.  4.9  ignores  the  break-up  of  cells  since  only  intact  cells 
are  analyzed.  Therefore,  only  a  mass  balance  on  the  overall  system  can  lead  to  the 
determination  of  the  overall  cell  neutralization  rate  (1  —  fi)  or  f[.  However,  Fig.  4.9 
provides  raw  data  on  the  collected  samples. 


(a)  (b) 

Figure  4.9:  Fraction  of  live  cells  in  different  collected  samples,  (a)  All  experiments 
start  with  a  measurement  of  the  fraction  of  live  cells  in  the  cell  culture  (susp),  the 
negative  control,  and  in  the  nebulized  cells  (neb).  Error  bars  show  the  95%  confidence 
interval  using  a  normal  distribution  for  5  measurements  per  sample,  (b)  The  cells 
that  pass  through  the  impactor  system  and  then  collected  in  the  collection  tube  at 
X  =  0.98,  x  —  0.50,  and  x  =  0.11  with  P0  =  1  atm.  The  ethanol  (EtOH)  suspension 
is  the  positive  control.  Error  bars  show  the  95%  confidence  interval  using  a  normal 
distribution  for  15  measurements  per  sample. 


The  percentage  of  live  cells  in  the  nebulized  sample  can  differ  from  that  in  the 
suspension  because  of  the  operating  conditions  of  the  nebulizer.  The  data  in  Fig.  4.9- 
a,  represents  5  experiments  with  different  starting  cultures.  The  fraction  of  live  cells 
in  the  suspension  (0.910  ±  0.014)  is  equal  to  that  of  the  nebulization  (0.938  ±  0.005) 
with  p  >  95%  using  a  two  tailed  t-test  with  unequal  variance.  Therefore,  the  effect 
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of  nebulization  can  be  neglected  in  the  calculation  of  the  neutralization  rate. 

In  Fig.  4.9-b,  15  experiments  with  different  starting  cultures  were  conducted  for 
each  of  the  operating  conditions  of  the  impactor.  Both  shock  conditions  =  O  il  and 
0.50)  have  statistically  equal  fraction  of  live  cells  while  they  are  lower  statistically  from 
the  non-shock  condition  —  0  98).  The  cells  exposed  to  ethanol  (positive  control) 
have  the  lowest  fraction  of  live  cells. 

4.4.2  Mass  balance 

To  evaluate  the  effect  of  particle  acceleration  in  an  aerodynamic  shock  a  mass  balance 
on  the  entire  experimental  system  has  to  be  done  to  account  for  all  losses.  Bacte¬ 
rial  aerosol  is  lost  to  the  walls  of  the  system  between  the  exit  of  the  nebulizer  and 
the  entrance  of  the  impactor.  Determining  the  losses  is  important  in  the  overall 
neutralization  rate  and  is  calculated  as  follows: 


(f)= —  (4.10) 

Tin 

where  nn  is  the  total  number  of  cells  that  are  nebulized  and  nx  is  the  total  number 
of  cells  that  enter  the  impactor.  Experimental  measurements  (see  Section  4.2.4) 
determine  the  actual  number  of  particles  that  enter  the  impactor  after  losses  in  the 
nebulization  process  and  in  the  tubes  connecting  the  nebulization  chamber  to  the 
impactor  inlet.  The  measured  0  is  0.088  ±  0.029  for  all  operating  conditions  of  the 


linpactor  system.  Losses  to  the  walls  of  the  nozzle  after  the  impactor  entrance  are 


neglected.  Nonetheless,  not  all  the  particles  that  enter  the  impactor  are  collected  for 
analysis  as  detailed  in  Section  4.4.1.  By  assuming  no  cell  break-up,  we  can  calculate 
the  number  of  live  cells  that  enter  the  collection  tube  (nfc)  as: 


n'c  =  nsl(pr) 


(4.11) 


where  n.,j  is  the  number  of  live  cells  inputed  to  the  system  by  nebulizing  the  sus¬ 
pension  and  is  obtained  experimentally  by  counting  the  cells  in  the  lower  right  half 
(L)  of  Fig.  4.8-b.  The  collection  efficiency  of  the  impactor  for  the  different  operating 
conditions  is  listed  in  Table  4.2.  Eq.  4.9  then  becomes: 


h  = 


nc 


(4.12) 


where  nc  is  obtained  by  counting  the  number  of  cells  in  the  lower  right  half  (L)  of 
Fig.  4.8-f  or  equivalent  data  for  the  different  impactor  operating  conditions  (Fig.  4.8-h, 

j)- 


4.4.3  Effect  of  the  aerodynamic  shock 

By  combining  fluorescence  cytometry,  microscopy  and  computational  collection  effi¬ 
ciencies  in  Eq.  4.10,  we  can  obtain  the  values  for  /;  for  the  different  operating  condi¬ 
tions  in  the  impactor  system.  The  data  is  shown  in  Fig.  4.10  with  the  95%  confidence 
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interval  on  the  normal  fit  of  15  independent  experiments  for  each  condition  and  us¬ 
ing  error  propagation.  In  the  experimental  impactor  system  operating  at  x  =  0.98, 
the  bacterial  aerosol  does  not  experience  the  critical  acceleration  discussed  in  Sec¬ 
tion  4.3.2:  therefore,  the  /;  is  an  order  of  magnitude  higher  than  the  conditions  which 
lead  to  acceleration  greater  than  the  critical  (>  3.00  x  108  m/s2).  The  acceleration  at 
X  =  0.50  is  2.02  x  10s  m/s2  which  is  approximately  equal  (considering  experimental 
variation  in  pressure  and  assumptions  in  calculating  the  computational  accelerations 
mention  in  Section  4.3)  to  the  critical  acceleration.  In  this  case,  a  reduction  in  cell 
viability  can  be  seen  that  is  statistically  the  same  as  that  observed  for  x  =  0.11  with 
an  acceleration  of  4.96  x  109  m/s2. 

It  is  important  to  note  that  the  Weber  numbers  in  Fig.  4.10  for  x  =  0.11,  0.50 
and  0.98  carried  at  Pq  —  1  atm  are  in  the  range  [0.0017, 3.5]  (see  Table  4.2)  which 
is  lower  than  the  computational  model-predicted  critical  Weber  number  of  14  needed 
for  bacterial  break-up.  However,  Fig.  4.10  clearly  demonstrates  that  this  estimate  is 
conservative  and  an  order  of  magnitude  reduction  in  cell  viability  can  be  attained  at 
a  lower  Weber  number  than  14  if  the  critical  acceleration  is  achieved.  To  investigate 
this  point  further,  wre  carried  out  experiments  (see  Remark  11)  at  a  Weber  number 
of  15  and  an  acceleration  of  5.98  x  109  m/s2  (both  in  the  range  needed  for  bacterial 
break-up)  with  x  =  0.11  at  Pq  —  5  atm  and  found  no  statistical  difference  in  fraction 
of  cells  remaining  live  (/)  =  0.047±0.016  at  P0  =  5  atm  compared  to  /;  =  0.035±0.012 
at  Pq  =  1  atm).  Therefore,  the  sharp  decelerations  that  the  cells  experience  after  they 
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Figure  4.10:  Fraction  of  cells  that  remain  alive  after  the  effect  of  impactor  operation 
( fi )  and  after  exposure  to  ethanol.  Live  cells  in  these  samples  are  compared  to  live 
cells  in  the  negative  control  by  accounting  for  losses  before  the  impactor  and  for  the 
collection  efficiency  of  the  impactor.  Error  bars  show  the  95%  confidence  interval 
using  a  normal  distribution  for  15  measurements  per  sample. 

pass  through  the  shock  is  the  main  reason  for  which  they  lose  viability. 

Remark  11  Five  experiments  were  carried  out  at  a  higher  upstream  pressure ,  Pq  = 
5  atm.  (We  =  15),  to  investigate  the  influence  of  Weber  number  on  cell  viability. 
In  order  to  run  the  impactor  system  of  Fig.  fl  at  an  upstream  pressure  of  Pq  = 
5  atm ,  a  high  pressure  chamber  was  designed  to  encase  the  nebulizer ,  and  the  HEPA 
filter  was  removed  to  avoid  exposure  to  atmospheric  pressure  before  the  impactor. 
Furthermore ,  the  nebulizer  was  operated  at  Pn  =  10  atm  to  maintain  the  same  Pn/Po 
ratio  as  the  atmospheric  pressure  experiments.  The  bends  in  the  tubing  connecting  the 
nebulizer  to  the  impactor  were  also  modified;  this  modification  might  slightly  change 
the  losses  as  measured  by  0,  however ,  in  our  calculations  we  assumed  the  losses  as 
calculated  in  Section  4-4-2  stayed  constant  for  all  experiments.  All  other  conditions 
were  maintained  as  described  in  Section  2.  The  fraction  of  live  cells  in  the  collection 
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tube  in  these  5  experiments  was  found  to  be  0.332  ±  0.215  and  the  corresponding  ft 
was  computed  to  be  fi  —  0.047  ±  0.016. 

4.5  Conclusions 

In  this  chapter,  the  computational  model  presented  in  Chapters  2  and  3  was  val¬ 
idated  by  experimental  viability  analysis  on  vegetative  gram-negative  E.  coli  cells. 
The  E.  coli  cells  were  used  as  a  model  bacterial  aerosol  with  a  computational-model 
predicted  critical  acceleration  of  3.0  x  108  m/s2  to  achieve  bacterial  break-up.  The 
experimental  impactor  system  was  operated  under  three  conditions  to  provide  max¬ 
imum  accelerations  of  4.96  x  109  m/s2,  2.02  x  108  m/s2  and  1.67  x  106  m/s2;  values 
above,  equal,  and  below  the  predicted  critical  acceleration  for  E.  coli ,  respectively. 
The  experiments  show  that  for  accelerations  higher  than  the  critical  value  the  shock 
has  a  considerable  effect,  fi  =  0.055  ±  0.020  for  x  =  0.50  and  fi  =  0.035  ±  0.012  for 
X  =  0.11.  Under  the  non-shock  condition,  ft  =  0.451  ±  0.150  which  is  an  order  of 
magnitude  higher  survival  rate  than  that  found  at  higher  accelerations.  These  results 
confirm  the  computational  model  predictions  of  Chapters  2  and  3.  The  computa¬ 
tional  model  can  therefore  be  used  to  predict  the  effect  of  an  aerodynamic  shock  on 
bacterial  cells  given  their  surface  tension  and  size. 
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Chapter  5 


Effect  of  Aerodynamic  Shock  on  B. 
atropheus  spores 

5.1  Introduction 

In  Chapter  4,  Gram-negative  E .  coll  vegetative  cells  were  shown  to  break-up  at  decel¬ 
erations  >  3.0 x  108  m/s2.  The  impactor  was  operated  to  reach  computed  decelerations 
of  5.0  x  109  m/s2  which  predict  break-up  of  the  cells.  The  experimentally  measured 
fraction  of  cells  remaining  live  was  0.035  ±  0.012,  confirming  the  computational  pre¬ 
dictions. 

This  Chapter  focuses  on  the  effect  of  aerodynamic  shocks  on  spores  which  are  more 
relevant  as  bioterrorism  threats.  In  laboratory  settings,  Bacillus  subtilis  var  niger , 
reclassified  as  Bacillus  atropheus  T91,  is  suggested  as  a  surrogate  for  B .  anthracis  that 
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emulates  the  spore’s  physical  characteristics  but  avoids  its  pathogenicity  [17].  The  two 
spores  have  similar  surface  tension  values  [8].  Both  computational  and  experimental 
methods  are  used  to  study  the  effect  of  aerodynamic  shocks  on  D.  atropheus  spores 
and  to  obtain  the  fraction  of  cells  that  remain  alive  after  the  effect  of  impactor 
operation  (ft). 

5.2  Experimental  materials  and  methods 

5.2.1  Experimental  setup 

The  experimental  system  discussed  in  Chapter  4  [49]  is  used  to  determine  the  effect 
of  the  aerodynamic  shock  on  B.  atropheus  spores.  The  details  of  the  system  (see 
Fig.  5.1)  are  described  below.  The  vegetative  E.  coli  suspension  was  substituted  with 
B.  atropheus  spore  suspension  (see  Section  5.2,2)  in  the  nebulizer  (see  Section  5.2.1) 
[49].  The  resulting  aerosol  travels  to  the  impactor  which  was  operated  at  different 
downstream  to  upstream  pressure  ratios  (see  Section  5.2.1)  to  capture  both  break-up 
and  no  break-up  conditions.  Quantitative  analysis  of  the  survival  rate  was  carried 
out  using  fluorescence  flow  cytometry  (see  Section  5.2.3). 

Nebulization 

The  B.  atropheus  spores  suspension  (see  Section  5.2.2)  was  fed  via  a  syringe  pump 
(KDS410,  Kd  Scientific)  at  a  liquid  flow  rate  Q  =  0.05  mL/min  for  20  min  to  a 
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Bacterial 

Suspension 


Figure  5.1:  Schematic  of  the  experimental  system.  (1)  Meinhard  nebulizer  with  a 
concentric  nozzle  mixing  a  liquid  bacterial  suspension  feed  and  a  dispersion  gas  (7V2) 
(2)  Nebulization  chamber  to  collect  large  droplets,  (3)  Converging  nozzle  with  an  exit 
diameter  of  d  —  0.5  mm,  (4)  and  (5)  a  flat  surface  with  a  0.5  mm  hole  combined  with 
a  screw  that  accommodates  the  collection  substrate  make  up  the  deceleration  tube, 

(6)  1.4  mm  spacer  fixes  the  spacing  x  between  the  impactor  and  nozzle  to  0.6  mm, 

(7)  support  holding  the  deceleration  tube  and  flat  surface  at  fixed  distance  from  the 
nozzle,  and  (8)  exit  chamber.  Major  losses  in  the  experimental  system  are  shown  in 
italic. 
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capillary  nebulizer  (TR-30-A1,  Meinhard  Glass  Products).  Nitrogen  (/V2)  gas  was 
supplied  at  a  critical  flow  rate  of  0.2  L/min  through  the  annular  region.  The  liquid  was 
fed  at  a  specified  bacterial  concentration  (see  Section  5.2.2)  determined  by  the  optical 
density  measurement  at  600  nm  (ODb oo).  Nitrogen  was  supplied  to  the  nebulizer  at  a 
pressure  Pn  via  a  pressurized  tank.  In  all  the  experiments,  the  nebulizer  was  operated 
at  Q  =  0.05  mL/min,  ODb 00  =  0.5,  and  Pn  =  2  atm.  The  spore  concentration  at 
OD 6oo  =  0.5,  as  measured  using  fluorescence  flow-cytometry  (see  Section  5.2.3),  was 
5.02  ±  1.24  x  105  spores/inL.  The  average  expected  nebulized  droplet  size  is  10  ^in, 
which  translates  to  approximately  1  cell  in  every  1000  droplets  assuming  no  aggregates 
are  formed  in  the  droplets. 

Impactor  system 

A  detailed  description  of  the  impactor  system  was  given  in  Chapters  2  and  4  [51,  49]. 
Briefly,  the  impactor  system  consists  of  a  plate  designed  for  the  collection  of  the 
bacterial  aerosol  (see  Fig.  5.1;  parts  4  and  5)  which  is  placed  perpendicular  to  the  gas 
flow  emerging  from  a  converging  nozzle  (see  Fig.  5.1;  part  3).  Iscntropic  flow  theory  of 
an  ideal  gas  predicts  a  critical  downstream  (Pt)  to  upstream  (P0)  pressure  ratio  (Xcnt) 
below  which  the  flow  at  the  exit  of  the  nozzle  is  sonic  [38,  48].  For  Nitrogen  (A2),  Xcrit 
is  0.53  [51,  49].  The  generated  impinging  flow  on  a  plate  from  a  converging  nozzle 
operating  under  sonic  conditions  results  in  the  creation  of  a  standoff  shock  whose 
properties  can  be  changed  by  varying  impactor  geometry  (distance  between  nozzle 
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and  plate  x  and  diameter  of  nozzle  d)  and  operating  conditions  [2,  13,  14,  29,  42,  51]. 
The  deceleration  tube  (see  Fig.  5.1;  parts  4  &  5)  is  filled  with  600  /xL  of  phosphate 
buffer  saline  (PBS)  solution  (BP2438-4,  Fisher  Scientific),  which  keeps  the  liquid  level 
to  about  4  mm  away  from  the  entrance  of  the  deceleration  tube.  The  plate  to  nozzle 
distance  is  set  at  x/d  =  1.2.  The  impactor  system  is  operated  at  an  upstream  pressure 
of  P0  =  1  atm  with  three  different  pressure  ratios  x  —  P\/Po  =  0.98,  0.50  and  0.11  for 
20  min.  After  each  experimental  run  the  volume  of  PBS  remaining  in  the  deceleration 
tube  was  measured  and  used  in  mass  balance  calculations  to  determine  the  fraction 
of  cells  surviving  the  impactor  operating  conditions  (see  Section  5.4),  Fig.  5.1  shows 
the  losses,  0  and  77,  of  bacterial  particles  in  different  parts  of  the  experimental  system 
which  are  discussed  in  Sections  5.2.1  and  5.3.2. 

Particle  losses  in  experimental  setup 

Losses  of  spores  to  the  walls  of  the  system  before  entering  the  impactor  (1  —  0)  were 
accounted  for  using  calculations  developed  Chapter  4  [49].  Specifically,  vegetative  E. 
coli  cells  were  collected  on  an  in-line  filter  placed  before  the  impactor  system  [49]. 
Fluorescence  microscopy  images  were  taken  to  count  the  fluorescent  bacteria  on  the 
stained  filter.  These  images  were  then  analyzed  to  calculate  the  fraction  of  cells 
entering  the  impactor  (0),  which  is  defined  as 
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where  nn  is  the  total  number  of  cells  nebulized  and  nt  is  the  total  number  of  cells 
entering  the  impactor.  The  parameter  0  was  calculated  to  be  0.088  ±  0.029  from 
the  in-line  filter  counts  [49].  Although  the  calculation  was  made  using  experimental 
data  on  E.  coli  vegetative  bacteria  (whose  diameter  is  ~  1  fim),  the  same  value  of  0 
was  used  for  B.  atropheus  spores  (whose  diameter  is  ss  0.5  ^m).  Two  processes  are 
responsible  for  losses  in  the  system  before  the  impactor:  (1)  losses  of  large  nebulized 
liquid  droplets  on  the  nebulization  chamber  (see  Fig.  5.1;  part  2)  walls,  and  (2)  losses 
of  bacterial  aerosol  in  the  tubes  leading  to  the  impactor.  Losses  in  the  tubes  can 
either  be  diffusional  or  inertial  (impaction).  Diffusional  losses  (1  —  nout/nlTl )  were 
calculated  for  circular  tubes  using  the  following  expression  [23]: 

flout /nin  =  1-5.5  [—)  +3.77^— J  (5.2) 


where  riout  and  Hjn  are  the  number  of  particles  exiting  and  entering  the  tube,  re¬ 
spectively,  D  is  the  diffusion  coefficient,  L  is  the  length  of  the  tube,  and  Qg  is  the 
volumetric  flow  rate  of  gas.  For  both  1  pm  and  0.5  pm  particles  the  ratio,  nm[/ntn, 
is  calculated  to  be  >  99%.  Therefore,  diffusional  losses  (1  —  nout/nm)  are  negligible. 
To  calculate  the  impaction  losses  the  following  relation  was  used  [23]: 


flout/  ff  in  —  1  2$ 


/  pAU'  \ 

V18  HgDj 


(5.3) 
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where  pv  is  the  particle  density,  dv  is  the  particle  diameter,  Ug  is  the  average  gas 
velocity  in  the  tubes  leading  to  the  impactor,  pg  is  the  gas  viscosity,  Ds  is  the  tube 
diameter  and  9  is  the  bend  angle.  The  term  in  parentheses  in  Eq.  5.3  is  also  referred 
to  as  the  Stokes  number.  The  ratio,  nout/nin,  was  also  >  99%  for  both  particle 
sizes.  Therefore,  most  of  the  losses  before  the  impactor  (1  -  <t>)  are  attributed  to  the 
identical  nebulization  process  (see  Section  5.2.1))  for  ~  1  /im  E.  coll  vegetative  cells 
and  ^  0.5  pm  B.  atropheus  spores  under  identical  conditions  with  0  =  0.088  ±0.029. 

5.2.2  Spore  preparation 

Bacillus  atropheus  (strain  ATCC  9372)  was  used  as  the  standard  bacterial  aerosol  to 
study  the  aerodynamic  shock  effect  on  bacterial  spores  in  the  impactor  system.  In 
our  previous  study,  E .  coli  was  used  as  a  test  aerosol  for  vegetative  bacterial  cells  [49]. 
The  methods  described  in  this  section  were  adopted  with  minor  modifications  to  the 
standard  procedures  for  culturing  and  washing  of  Bacillus  spores  [20].  All  vegetative 
cultures  were  initiated  from  a  stock  B.  atropheus  suspension  (80%  glycerol  solution) 
stored  at  —20  °C.  5  mL  of  Nutrient  broth  (Fisher)  was  inoculated  with  10  pL  of  stock 
suspension  in  a  test  tube.  The  vegetative  culture  was  incubated  in  a  shaker  (Excella 
E24,  New  Brunswick  Scientific)  at  32  cC  and  200  rpm  for  approximately  22  —  25  hrs. 
To  start  the  spore  culture,  1  mL  of  the  vegetative  B.  atropheus  was  inoculated  in 
15  mL  of  2xSG  media  that  was  prepared  as  described  in  standard  procedures  [20]. 
The  spore  culture  was  incubated  for  12  days  at  32  °C  and  200  rpm,  allowing  time  for 
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the  cells  to  exhaust  nutrients  in  the  medium,  thereby  inducing  sporulation. 

After  the  incubation  period  the  spore  culture  was  washed  over  4  days.  On  day 
1,  the  culture  was  centrifuged  at  4400  rpm  for  15  min  (Centrifuge  5702,  Eppendorf). 
Three  layers  were  observed  with  the  krwer  layer  containing  the  free  spores  and  the 
top  2  layers  containing  cell  debris  [20].  The  suspension  was  washed  3  times  with  cold 
autoclaved  deionized  (DI)  water,  keeping  all  three  layers  and  then  shaking  in  an  ice 
bath  overnight.  Similar  procedures  were  repeated  for  days  2  and  3  with  the  exception 
of  resuspending  and  decanting  the  top  two  layers  instead  of  keeping  them.  Most  of 
the  cell  debris  layers  were  removed  after  the  third  day.  On  day  4  the  suspension  was 
centrifuged  at  4400  rpm  and  wrashed  3  times  using  cold  PBS.  The  optical  density  of 
the  spore  suspension  wras  measured  at  a  wavelength  of  600nm  (ODqqq)  using  a  visible 
light  spectrophotometer  (Genesis  20,  Thermo  Scientific)  before  storage  at  4  °C. 

To  check  the  degree  of  sporulation,  starting  suspensions  w'ere  treated  with  lysozyme, 
which  breaks  dowrn  the  peptidoglycan  outer  membrane  in  Gram-positive  vegetative 
cells  such  as  B.  atropheus.  Since  spores  do  not  have  an  exposed  peptidoglycan  layer, 

4 

lysozyme  is  expected  to  have  no  effect.  The  suspension  was  diluted  to  OD$ oo  —  0.05 
and  lysozyme  (L6876,  Sigma  Aldrich)  was  added  for  a  solution  concentration  of 
2.8  x  10-1  g/mL  and  incubated  for  30min  at  32  cC  and  200  rpm.  The  lysozyme 
treatment  was  also  used  on  vegetative  B.  atropheus.  The  fraction  of  spores  was  de¬ 
termined  by  comparing  viability  analysis  (see  Section  5.4.1)  of  the  spore  suspension 
to  the  vegetative  cells. 
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At  the  start  of  each  experiment,  a  positive  and  a  negative  control  were  prepared 
by  diluting  the  spore  suspension  to  OD600  =  0.05.  The  positive  control  was  incu¬ 
bated  in  a  water  bath  at  100  °C  for  over  2  hours  to  neutralize  the  spores  11  The 
negative  control  was  left  untreated  and  is  therefore  expected  to  have  the  highest  con¬ 
centration  of  viable  spores.  The  samples  collected  after  the  impactor  operation  is 
compared  to  both  the  positive  and  negative  control  to  assess  the  effect  of  the  shock 
(see  Section  5.4.2). 

5.2.3  Bacterial  viability  analysis 

Bacterial  viability  is  determined  using  LIVE/DEAD(g)^acLight™  (cat.  number 
L7012;  Life  Technologies  Corporation)  [3,  32].  The  assay  contains  two  components: 
(A)  3.34  mM  STYO  9,  a  green-fluorescent  nucleic  acid  stain,  in  dimethyl  sulfoxide 
(DMSO)  and  (B)  20  mM  propidium  iodide  (PI),  a  red-fluorescent  nucleic  acid  stain, 
in  DMSO.  STYO  9  penetrates  both  intact  and  damaged  membranes,  while  PI  only 
penetrates  damaged  membranes.  A  stock  dye  solution  is  prepared  at  the  start  of  ev¬ 
ery  experiment  by  mixing  5  /iL  of  each  of  components  A  and  B  with  10  /iL  deionized 
water  (DI).  All  controls  and  experimental  samples  mentioned  in  Sections  5.2.1  and 
5.2.2  were  stained  with  the  stock  dye  solution  to  a  final  concentration  of  8.35  /iM  and 
50  (jM  of  STYO  9  and  PI,  respectively. 

The  stained  samples  were  analyzed  using  a  fluorescence  flow  cytometer  (FACScan™, 
BD  Biosciences).  The  cytometer  was  operated  at  12  /iL/min  for  20s  with  a  488  nm 
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excitation  laser.  Red  fluorescence  was  measured  with  a  high  pass  filter  at  630  inn 
(FL2-H)  and  green  fluorescence  was  measured  with  a  bandpass  filter  of  520  nm  (FL1- 
H).  The  forward  scatter  (FSC)  and  side  scatter  (SSC)  were  also  measured  for  the 
particles  passing  the  laser.  Photomultiplier  tubes  were  used  to  amplify  the  optical 
signal  for  FL1-H,  FL2-H  and  SSC  at  730,  648  and  462  volts,  respectively.  A  photo¬ 
diode  was  used  to  detect  the  FSC  signal. 

5.3  Computational  results 

5.3.1  Gas  flow  field 

At  low  aerosol  concentrations  the  equations  for  gas  and  particle  dynamics  are  one¬ 
way  coupled  [18].  First,  the  gas  dynamics  equations  are  solved  and  the  gas  velocity  is 
used  to  calculate  the  particle  dynamics.  A  FLUENT®  computational  fluid  dynamics 
(CFD)  code  was  developed  to  solve  both  gas  and  particle  dynamics  in  the  experimen¬ 
tal  impactor  system  of  Fig.  5.1  51].  Details  of  the  model  were  described  in  Chapter  2. 
The  Reynolds  time-averaged  Navier-Stokes  equations  in  Cartesian  tensor  form  are: 

‘5-4> 
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(5.5) 


where  p  is  the  density,  v  is  the  velocity,  p  is  the  static  pressure,  6  is  the  Kronecker 
delta,  /i  is  the  viscosity,  and  i,j,  and  l  are  the  Cartesian  indices.  To  achieve  closure  of 


Eq.  5.5.  the  Reynolds  stress  term  —pv[vf  is  calculated  using  the  Boussinesq  hypothesis 
[24].  A  standard  k  —  e  turbulence  model  is  used  in  the  gas  flow  equations.  A  detailed 
computational  parametric  study  indicated  that  the  maximum  centerline  Mach  number 
(Mmax  increases  with  increasing  x/d  and  decreasing  P\/Pq)  (see  Chapter  2)  r51 
Furthermore,  at  higher  x/d  (>  1.4)  the  aerodynamic  shock  is  unstable  [51]  x/d  =  1.2 
had  the  highest  Mmax  with  no  flow  instabilities  (see  Chapter  2)  [51]  Operating 
the  irnpactor  at  three  different  pressure  ratios  \  =  R98  0.50,  and  0.11  captures 
gas  flow  fields  with  no  aerodynamic  shock,  at  critical  nozzle  flow  (Xcrit)^  and  strong 
aerodynamic  shock,  respectively.  The  maximum  centerline  Mach  number  is  2.36  at 
X  =  0.11  (see  Table  5.1). 


Table  5.1:  Aerodynamic  shock  properties  at  different  operating  conditions  with  x/d  = 
1.2  and  Po  —  1  atm _ 


X 

Mmax 

Re 

0.98 

0.14 

9.9  x  102 

0.50 

0.96 

5.8  x  103 

0.11 

2.36 

5.4  x  103 
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5.3.2  Particle  dynamics 


The  velocity  profile  of  the  particle  motion  is  used  for  the  calculation  of  mechanical 
forces  induced  oil  the  particle  (see  Section  5.3.2).  As  in  previous  chapters,  the  particle 
trajectory  is  solved  using  Newtons  second  law  with  the  gas  flow  field  solutions  from 
Section  5.3.1  [49,  51]: 


mp 


dvp 

dt 


2  C DPgAp(Vg  Vp)\vg  Vp\) 


(5.6a) 


—  =  Vp,  (5.6b) 

where  vp  and  vg  are  the  velocity  of  the  bacterial  particle  and  gas.  respectively,  mp 
is  the  mass  of  the  particle,  pg  is  the  density  of  the  gas,  Ap  is  the  projected  area  of 
the  bacterium  on  the  plane  perpendicular  to  the  flow  direction,  and  Co  is  the  drag 
coefficient.  The  thermophoretic,  gravitational,  and  Basset  history  term  are  neglected 
in  Eq.  5.6a.  The  drag  force  is  assumed  to  be  the  major  force  acting  on  the  bacterial 
particle  which  is  typical  for  impactor  systems.  A  high  Mach  number  drag  law  is  used 
to  calculate  the  drag  coefficient  (Co)  [9].  The  size  of  Bacillus  spores  is  in  the  range 
of  0.5  —  1.2  /zm  [10,  8].  The  bacterium  is  simulated  as  a  0.5  jam  sphere  with  a  density 
Pp  of  lg/mL  which  is  representative  of  values  reported  in  literature  [33,  59].  Values 
for  1  /am  particles  were  obtained  from  previous  studies  [49].  The  computational 
simulations  mimic  the  experimental  conditions  in  the  study  with  x/d  =  1.2  under 
shock  (x  =  0.11,0.50)  and  non-shock  (x  —  0.98)  conditions.  The  detailed  conditions 
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are  shown  in  Sections  5.3.2  and  5.3.2.  The  computational  results  predict  the  collection 
efficiency  of  the  impactor  system  at  different  operating  conditions  (see  Section  5.3.2) 
and  bacterial  envelope  instability  (see  Section  5.3.3). 

Particle  acceleration 

A  fifth-order  Runge-Kutta  method  is  used  to  numerically  integrate  Eq.  5.6a  and 
5.6b  using  the  gas  velocity  obtained  from  the  computational  model  described  in 
Section  5.3.1.  The  particle  velocity  is  calculated  every  5  x  10-5  —  1CT7  m  (for 
X  =  0.11  —  0.98)  along  the  centerline  of  the  impactor.  Fig.  5.2  shows  the  gas  and 
particle  velocity  profiles  along  the  centerline  of  the  impactor  in  the  region  of  the 
shock.  Under  shock  conditions  (x  =  0.11),  the  aerodynamic  shock  is  observed  as  a 
sharp  drop  in  gas  velocity  (see  Fig.  5.2-a).  Because  the  particle  does  not  adjust  to  the 
gats  stream  velocity,  it  experiences  a  relative  deceleration.  For  non-shock  conditions 
(X  —  0.98),  the  gas  velocity  decreases  at  slower  rate  compared  to  the  shock  allowing 
the  particle  to  adjust  its  velocity  (Fig.  5.2-c)  without  experiencing  a  relative  decelera¬ 
tion  as  strong  as  in  the  shock  case.  Fig.  5.2-b  shows  the  velocity  profiles  for  \  =  0.50 
where  the  particle  observes  an  intermediate  relative  deceleration. 

The  relative  deceleration  of  the  particle  is  determined  by  calculating  the  relative 
velocity  derivative  with  respect  to  time  (see  Eq.  3.8b)  [51  .  The  results  of  this  differen¬ 
tiation,  expressed  in  absolute  values,  are  shown  in  Fig.  5.3.  The  decelerations  provide 
a  measure  of  the  force  required  to  neutralize  bacterial  aerosols  (see  Section  5.3.3).  As 
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Figure  5.2:  Centerline  particle  and  gas  velocity  as  a  function  of  the  distance  from  the 
nozzle.  The  entrance  of  the  collection  tube  is  at  a  distance  6  x  10"4  m  from  the  nozzle 
exit.  Impactor  operating  conditions:  (a)  x  =  0.11,  (b)  \  =  0.5,  and  (c)  \  —  0*98. 

seen  in  Fig.  5.3-b  and  c,  the  absolute  value  of  acceleration  increases  at  5  x  10"4  rn  from 
the  nozzle  exit  after  the  main  deceleration  peak.  Since  Brownian  motion  and  gravita¬ 
tional  effects  were  neglected  in  Eq.  5.6a,  particles  artificially  stop  before  reaching  the 
collection  tube  bottom.  This  is  manifested  as  a  second  increase  in  the  absolute  value 
of  the  acceleration  at  the  entrance  of  the  deceleration  tube.  In  Chapter  4,  simulations 
with  1  /im  particles  showed  similar  increases  in  the  acceleration  before  the  inlet  of 
the  deceleration  tube  because  of  assumptions  in  the  computational  model  discussed 
in  Section  5.3.2  [49]. 
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(a)  (b) 


Figure  5.3:  Bacterial  particle  acceleration  (relative  to  gas  flow)  as  a  function  of  the 
distance  from  the  nozzle.  The  entrance  of  the  collection  tube  is  at  a  distance  6x10  4  m 
from  the  nozzle  exit.  Impactor  operating  conditions:  (a)  x  ~  0.11,  (b)  x  =  0-5.  and 
(c)  X  -  0.98. 

The  maximum  centerline  decelerations  ( amax )  reached  between  the  nozzle  exit  and 
deceleration  tube  entrance  for  the  three  different  operating  conditions  are  reported 
in  Table  5.2.  The  maximum  acceleration  and  weber  number  increase  with  decreasing 
X  (see  Chapters  3  and  4)  [49,  51].  At  x  —  0.11,  acceleration  increases  with  smaller 
particle  sizes,  which  is  a  trend  observed  for  x  —  0.14  at  x/d  =  1.0  and  1.4  in  Chapter  3 
[51].  At  x  ~  0.98  and  0.50,  the  trend  is  reversed  as  the  flow  transitions  to  subsonic 
velocities. 
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Particle  collection  efficiency 


In  Section  5.3.2  average  gas  velocities  were  used  to  calculate  the  particle  acceleration. 
However,  the  collection  efficiency  for  highly  turbulent  flows  and  low-density,  small- 
size  aerosols,  depends  on  particle  dispersion  due  to  gas  velocity  fluctuations.  The 
collection  efficiency  (fj)  is  defined  as  the  ratio  of  the  number  of  particles  that,  pass 
through  the  collection  tube  entrance  compared  to  the  number  of  particles  that  enter 
the  impactor  inlet.  The  bacterial  collection  efficiency  (77)  for  all  three  impactor  oper¬ 
ating  conditions  (y  =  0.11, 0.50,  0.98)  are  evaluated  to  determine  the  overall  bacterial 
aerosol  neutralization  rate  in  Section  5.4.  The  random  walk  model  is  employed  to 
study  the  effect  of  velocity  fluctuation  on  bacterial  aerosol  with  reasonable  computa¬ 
tional  load.  Stochastic  tracking  of  the  particles  is  used  to  account  for  the  turbulent 
dispersive  effects  in  the  gas  flow.  In  the  stochastic  tracking  model  the  gas  velocity 
is  represented  as  an  average  plus  a  fluctuating  term  when  calculating  the  particle 
velocity,  vV)  in  Eq.  5.6a.  The  fluctuating  velocity  is  kept  constant  for  time  intervals 
corresponding  to  the  lifetime  of  turbulent  eddies.  The  following  equations  describe 
the  stochastic  behavior: 


vg  =  v  +  v 


(5.7a) 


(5.7b) 


tL  =  0.15/c/e 


(5.7c) 


t(  =  -tL  log  r 


(5.7d) 
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where  vg  is  the  instantaneous  gas  velocity,  v  is  the  average  velocity,  and  v*  is  the 
velocity  fluctuation  term.  vf  is  a  function  of  (,  a  normally  distributed  random  number, 
and  va  is  the  local  root  mean  square  of  the  fluctuations,  v'2  is  calculated  based  on 
the  assumption  of  isotropy  and  is  given  by  va  —  2kj3 ,  with  k  being  the  turbulence 
kinetic  energy  in  the  k—e  model,  ti  is  the  lagrangian  time  interval,  e  is  the  turbulence 
dissipation  rate,  rc  is  the  lifetime  of  the  eddies  and  r  is  a  uniformly  distributed  random 
number  between  0  and  1  Fig.  5.4  shows  the  path  of  5  particles  released  from  the 


Figure  5.4:  Representative  trajectories  of  5  particles  released  from  the  inlet  of  the 
impactor  for  (a)  \  “  0.11  (b)  \  “  0.5  and  (c)  x  =  0.98.  The  collection  efficiencieb  as 
measured  by  the  number  of  particles  entering  the  collection  tube  are  0.37,  0.48  and 
0.26  for  x  —  0.11,  x  =  0.50  and  x  =  0.98,  respectively. 
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inlet  as  they  exit  the  nozzle  of  the  impactor  system  for  both  x  ~  0.11  (Fig.5.4-a), 
X  =  0.50  (Fig.  5.4-b)  and  x  =  0.98  (Fig.5.4-c).  In  calculating  77,  10,000  particles  are 
released  at  the  inlet  of  the  impactor.  The  particles  are  then  counted  as  they  pass  the 
virtual  surfaces  s  -  1  (nozzle  exit)  and  s  —  2  (deceleration  tube  entrance),  shown  in 
Fig.  5.5. 


Figure  5.5:  Computational  domain  of  the  impactor.  The  two  dashed  lines  indicate 
different  virtual  surfaces  (S-l  and  S-2)  where  particles  are  counted.  The  collection 
efficiency  is  calculated  by  dividing  the  number  of  particles  passing  S-2  by  that  of  S-l. 


Except  for  the  pressure  inlet  and  outlet,  all  the  boundaries  shown  in  Fig.  5.5  are 
assumed  to  be  reflection  boundaries.  The  collection  efficiency  is  defined  as  follows: 


(5.8) 


where  n,  is  the  number  of  particles  entering  the  impactor  and  n'c  is  the  number  of 
particles  that  enter  the  collection  tube.  n'c  is  calculated  under  the  assumption  of 
no  break-up.  77  is  used  to  determine  the  total  number  of  live  cells  that  would  be 
collected  if  the  shock  were  to  have  no  effect  (see  Section  5.4).  For  different  operating 
conditions,  the  collection  efficiency  (77)  for  both  1  fj,m  and  0.5  gm  size  particles  is 
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reported  in  Table  5.2.  The  values  for  1  pm  are  taken  from  Chapter  4  [49].  r/  increases 
with  diameter  as  the  inertia  of  the  particle  increases  and  prevents  it  from  following 
the  gas  streamlines  [21,  22  . 


Table  5.2:  Particle  dynamics  parameters  at  different  operating  conditions  with  x/d 
1.2  and  Pq  =  1  atm  for  dv  =  0.5  and  1  ^m 


X 

^ max 

0.5  fim 

(m/s2) 

1  /i  m 

We. 

0.5  fim 

max 

1  fim 

V 

0.5  pm 

1  fim 

0.98 

4.7  xIO^ 

1.7  x  1(T 

5.6  x  10"4 

1.70  x  10'2 

0.26 

0.35 

0.50 

1.5  x  108 

2.0  x  108 

5.4  x  10-1 

3.5  x  10° 

0.48 

0.72 

0.11 

5.9  x  109 

5.0  x  109 

1.7  x  10° 

3.1  x  10° 

0.34 

0.51 

5.3.3  Critical  accelerations  for  B.  atropheus 

Rayleigh-Taylor  instabilities  are  created  at  interface  between  two  fluids  of  different 
density  (air-bacterium)  when  they  are  at  relative  acceleration  with  respect  to  each 
other  [7,  28].  As  shown  in  Chapters  3  and  4,  The  acceleration  induces  waves  on  the 
interface  which  become  unstable  at  a  critical  acceleration  of: 


ac  =  47r2 


Pvdl 


(5.9) 


where  ac  is  the  critical  acceleration,  a  is  the  surface  tension  of  the  bacterium.  dp  is 
the  diameter  of  the  bacterium,  and  pp  is  the  density  of  the  bacterium.  The  Weber 
number  of  the  bacterium  should  also  exceed  the  critical  Weber  number  [28]: 


We  >  Wec  =  12(1  +  1.077 Oh1-6) 


(5.10a) 


100 


where 


and 


Wc=  PgK  ~%)2^P 
a 


(5.10b) 


Oh  = 


MP 

(PpdpCr)1/2 


(5.10c) 


where  We  is  the  Weber  number,  which  is  the  ratio  between  the  inertial  force  exerted 
on  the  bacterial  particle  and  the  particle  surface  tension  force.  Oh  is  the  Ohnesorge 
number,  which  is  the  ratio  between  the  viscous  forces  and  the  surface  tension  force. 
vPt  and  v9t  are  the  particle  and  gas  velocity  in  the  ^-direction  (i.e.,  direction  of  flow 
along  centerline),  respectively.  Therefore,  bacterial  membrane  break-up  is  inferred 
by  computing  the  relative  velocity  (vPx  —  v9i)  of  the  bacterium  and  the  maximum 
acceleration  achieved.  The  critical  acceleration  conditions  that  satisfy  the  conditions 
for  break-up  are  met  for  x  =  0.11;  however,  the  conditions  for  the  Weber  number 
are  not.  As  reported  in  Sections  5.4.2  and  4.4.3  [49],  cells  will  still  break-up  un¬ 
der  non-critical  Weber  number  conditions;  which  points  to  the  conservative  nature 
of  the  Weber  number  conditions.  In  this  study,  only  the  effect  of  Rayleigh-Taylor 
instabilities  is  considered. 


Table  5.3:  Cell  properties  for  E.  coli  and  5.  atropheus  [8,  61]  and  critical  shock 
properties  needed  to  induce  bacterial  break-up.  Values  are  reported  for  both  0.5  jum 
and  1  /v,m  for  spores.  (The  Gram  column  refers  to  the  Gram  stain) 


Cell 

Gram 

State 

dp  (/urn) 

a  (N/m) 

ac  (m/s2) 

Wec 

E.  colt 

-ve 

Vegetative 

1 

7.5  x  10"3 

3.0  x  108 

1.4  x  101 

B.  atropheus 

-fve 

Spore 

0.5 

9.6  x  10"2 

1.6  x  1010 

1.3  x  101 

B.  atropheus 

-f  ve 

Spore 

1 

9.6  x  lO"2 

3.9  x  109 

1.2  x  101 
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Bacterial  cells  with  a  lower  surface  tension  and  larger  diameters  require  lower 
accelerations  to  break-up.  It  is  important  to  note  that  the  surface  tension  values 
reported  in  Table  5.3  are  derived  from  the  turgor  pressure  [61]  and  contact  angle  [8]. 
However,  the  surface  tension  can  vary  with  deformation  generated  by  fluid  flow  over 
the  cell  membrane.  Such  a  variation  of  the  cell  surface  tension  under  fluid  flow  may 
limit  the  accuracy  of  the  critical  acceleration  predictions  of  Table  5.3.  In  such  a  case, 
more  accurate  models  of  the  variation  of  cell  surface  tension  under  stress  generated  by 
the  fluid  flow  (similar  to  the  works  of  [30,  31,  52,  60])  should  be  utilized  and  combined 
with  the  fluid  flow  and  particle  motion  models  developed  here  to  make  accurate  critical 
acceleration  predictions.  Previous  work  [31]  has  demonstrated  for  the  first  time  the 
use  of  inertial  forces  (expressed  in  terms  of  the  Weber  number  and  acceleration)  in 
achieving  cell  break-up  in  the  context  of  a  high-pressure  homogenization  system. 

5.4  Experimental  results  and  discussion 

5.4.1  Flow  cytometry  results 

Fig.  5.6  shows  a  representative  set  of  raw  data  obtained  using  the  protocols  described 
in  Section  5.2.3.  A  total  of  18  experiments  were  performed  with  6  different  starting 
spore  suspensions.  In  the  left  plots,  the  side  scatter  (SSC)  versus  the  forward  scatter 
(FSC)  (see  Fig.  5.6)  is  a  measure  of  the  relative  particle  sizes.  The  S  gate  indicates 
the  population  of  spores.  The  right  plots  in  Fig.  5.6  show  the  red  (FL2-H)  versus 
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green  (FL1-H)  fluorescence  emission  (see  Section  5.2.3)  of  the  spores  in  gate  S.  The 
negative  (Fig.  5.6-a,b)  and  positive  (Fig.  5.6-c,d)  control  define  the  live  and  dead  cell 
region  in  the  FL2-H  vs.  FL1-H  plots.  The  live  regions  are  indicated  by  the  gate  L  and 
the  dead  by  D.  The  fraction  of  cells  live  ( L/S )  is  determined  by  dividing  the  number 
of  spores  in  gate  L  by  the  total  number  of  spores  in  gate  S  ( S  =  L  +  D).  The  average 
value  of  the  fraction  of  live  cells  {L/S)  from  all  18  experiments  is  shown  in  Fig.  5.7 
and  Table  5.4.  L/S  is  not  the  true  measure  of  the  effect  of  the  aerodynamic  shock  on 
bacterial  aerosol  since  it  only  measures  intact  cells.  The  overall  neutralization  1  —  /; 
or  fraction  of  cells  remaining  intact  after  the  effect  of  the  impactor  fi  is  calculated 
through  a  mass  balance  on  the  experimental  system  and  is  discussed  in  Section  5.4.2. 

Using  a  two-tailed  t-test  with  unequal  variance,  the  fraction  of  live  cells  for  the 
negative  control  and  x  =  0-98  is  equal  with  a  p  >  95%  (see  Table  5.4).  However, 
the  mean  fraction  of  cells  live  are  unequal  for  x  —  0.98,  x  =  0.50,  and  x  =  0.11 
(L/5098  >  L/S050  >  L/S0  n)  with  a  p  >  95%.  In  Chapter  4  with  E.  coli,  L/V09g  > 
L/V q  jo  —  L/V 011  (w'here  V  —  L  D  represents  the  vegetative  cells  gate)  with  a 
p  >  95%  [49].  The  results  are  discussed  in  Section  5.4.2  in  comparison  to  the  fraction 
of  cells  that  remain  live  after  the  shock  (/;),  w-hich  accounts  for  the  total  mass  balance 
in  the  system. 
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Figure  5.6:  (left)  Side-scatter  versus  forward-scatter  plot  of  all  particles  in  the  flow 
cytometer.  The  S  gate  indicates  the  particle  sizes  which  represent  intact  spores.  PBS 
debris  indicates  the  debris  population,  (right)  Red  fluorescence  (FL2-H)  versus  green 
fluorescence  (FL1-H)  of  gated  cells  in  the  left  side  plots.  The  L  gate  represents  the 
live  cells  and  the  D  gate  represents  the  dead  cells.  (a,b)  spore  suspension  in  PBS  at 
OD$q o  =  0.05,  (c,d)  spore  suspension  after  exposure  to  heat,  (e,f)  cells  collected  after 
passing  impactor  operated  at  x  =  0.98,  (g,h)  x  =  0.50  and  (ij)  \  =  0.11. 
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Figure  5.7:  Fraction  of  live  cells  in  different  collected  samples  (L/S).  All  experiments 
start  with  a  measurement  of  the  fraction  of  live  cells  in  the  cell  culture  (live),  and 
the  negative  control  (heat).  The  cells  that  pass  through  the  impactor  system  and 
then  collected  in  the  collection  tube  at  x  =  0.98,  x  =  0.50,  and  x  =  0.11  with 
Po  =  1  atm.  Error  bars  show  the  95%  confidence  interval  using  a  normal  distribution 
for  18  measurements  per  sample. 
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Lysozyme  treatment  results 


The  effect  of  lysozyme  on  spores  (Fig.  5.8-c  and  d)  is  compared  to  the  effect  of 
lysozyme  on  vegetative  B.  atropheus  cells  (Fig.  5.8-g  and  h)  in  order  to  confirm  the 
presence  of  spores  in  the  starting  suspensions.  The  live  (L)  and  dead  (D)  gates  were 
set  by  running  a  positive  control  for  both  spore  and  vegetative  cells  (not  shown  in 
Fig.  5.8).  Ethanol  is  used  as  a  positive  control  for  vegetative  cells  as  in  Chapter  4  [49]. 
The  V  gate  indicates  the  vegetative  cell  population  versus  the  S  gate  for  the  spore 
population.  No  shift  in  the  spore  population  is  seen  after  treatment  with  lysozyme 
in  Fig.5.8-a  and  c.  For  vegetative  B.  atropheus  cells  the  forward  scatter  decreases 
(Fig.5.8-e  vs.  g)  after  lysozyme  treatment  indicating  smaller  cell  sizes  or  lysis. 

Furthermore,  the  vegetative  cells  treated  with  lysozyme  shift  to  the  D  gate.  Com¬ 
paring  both  the  scatter  and  fluorescence  plots  for  spores  and  vegetative  cells  it  is 
clear  that  the  starting  spore  suspensions  are  not  affected  by  lysozyme  treatment. 
The  fraction  of  live  spore  cells  ( L/S )  after  lysozyme  treatment  is  0.868  ±  0.253  (4 
measurements)  and  0.020  ±0.114  (2  measurements)  for  vegetative  cells  (L/V)  (see 
Fig.  5.9).  These  values  are  statistically  unequal  with  p  >  95%.  The  mean  value  of 
L/S  for  the  starting  spore  suspension  and  spores  treated  with  lysozyme  were  statis¬ 
tically  equal  (p  >  95%)  which  suggests  that  most  of  the  nebulized  cells  were  indeed 
spores. 
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Figure  5.8:  (left)  Side-scatter  versus  forward-scatter  plot  of  all  particles  in  the  flow 
cytometer  The  S  gate  indicates  the  particle  sizes  which  represent  intact  spores. 
The  V  gate  represents  the  intact  vegetative  cells.  PBS  debris  indicates  the  debris 
population,  (right)  Red  fluorescence  (FLH-2)  versus  green  fluorescence  (FLH-1)  of 
gated  cells  in  the  left  side  plots.  The  L  gate  represents  the  live  cells  and  the  D 
gate  represents  the  dead  cells,  (a.b)  spore  suspension  in  PBS,  (c,d)  spore  suspension 
with  lysozyme  treatment,  (e,f)  vegetative  suspension  in  PBS  and  (g,h)  vegetative 
suspension  with  lysozyme  treatment. 
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Figure  5.9:  Fraction  of  live  cells  (L/S)  after  lysozyme  treatment  for  both  spores 
and  vegetative  cells.  Error  bars  show  the  95%  confidence  interval  using  a  normal 
distribution  for  4  measurements  for  spores  and  2  measurements  for  vegetative  cells. 
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5.4.2  Effect  of  the  aerodynamic  shock  on  B.  atropheus  spores 


The  fraction  of  cells  that  remain  alive  after  the  effect  of  impactor  operation  (/*) 
provides  a  quantitative  measure  of  aerodynamic  shock  effect  when  compared  to  the 
operation  of  the  impactor  under  non-critical  acceleration  conditions.  Fluorescence 
flow  cytometry  measurements  (see  Section  5.2.3)  and  computational  particle  tracking 
(see  Section  5.3.2)  are  used  to  calculate  fi  using  the  following  expression  [49]: 


n c  _  nc 

K.  nsl<S*n 


(5.11) 


where  nc  is  the  number  of  live  cells  collected  after  the  impactor  operation  as  measured 
by  flow  cytometry  (see  Section  5.4.1)  and  is  equal  to  the  cells  in  gate  L  multiplied  by 
the  volume  of  PBS  remaining  in  the  collection  tube,  n'  is  the  number  of  cells  live  that 
would  be  collected  if  the  impactor  operation  had  no  effect.  n'c  can  not  be  determined 
from  direct  measurements  because  of  the  aerodynamic  shock  effect.  A  mass  balance 
accounting  for  the  losses,  shown  in  Fig.  5.1,  has  been  used  where  0  is  the  particle 
loss  before  the  impactor  (see  Section  5.2.1)  and  q  is  the  particle  collection  efficiency 
(see  Section  5.3.2).  n'c  becomes  nsi<*>q.  nsi  is  the  number  of  live  cells  in  suspension 
as  measured  by  flow  cytometry  (see  Section  5.4.1)  and  is  equal  to  the  cells  in  gate 
L  multiplied  by  1  mL  (volume  of  nebulized  suspension).  The  particle  losses  before 
the  impactor  are  calculated  in  Section  5.2.1  and  are  given  as  0  =  0.088  ±  0.029.  The 
collection  efficiency  (q)  for  both  0.5  and  1  /im  spores  in  the  impactor  system  is  given 
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in  Table  5.2.  As  seen  in  Fig.  5.10  and  Table  5.4,  fi  decreases  with  increasing  r)  (see 
Remark  13).  Both  fi  values  are  calculated  and  reported  in  Table  5.4  and  Fig.  5.10. 
A  detailed  derivation  of  Eq.  5.11  is  shown  in  our  previous  work  49]. 
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Figure  5.10:  Fraction  of  cells  that  remain  alive  after  the  effect  of  impactor  operation 
(/;)  and  after  exposure  to  heat.  Live  cells  in  these  samples  are  compared  to  live 
cells  in  the  negative  control  by  accounting  for  losses  before  the  impactor  and  for  the 
collection  efficiency  of  the  impactor  for  (a)  dp  =  0.5  pm  and  (b)  dp  =  1  pm.  Error  bars 
show  the  95%  confidence  interval  using  a  normal  distribution  for  18  measurements 
per  sample. 


The  critical  acceleration 
3.9  x  109  m/s2,  respectively 


(ac)  for  0.5  and  1  pm  spores  are  1.6  x  1010  m/s2  and 
(see  Table  5.3).  At  x  =  0.98,  the  maximum  computed 
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Table  5.4;  Measured  fraction  of  live  cells  ( L/S )  and  the  fraction  of  spores  surviving 
the  shock  ( fi )  for  both  dp  =  0.5  and  1  jim.  The  fi  for  the  -fve  control  is  not  size 
dependent  but  is  shown  under  the  0.5  /im  column  for  convenience. _ 


Condition 

L/S 

//(0.5/um) 

-ve  control 

0.971  ±0.012 

- 

- 

X  =  0.98 

0.983  ±0.011 

0.715  ±0.235 

0.531  ±0.176 

X  =  0.50 

0.895  ±  0.036 

0.033  ±0.011 

0.022  ±0.008 

X  =  0.11 

0.765  ±  0.082 

0.030  ±0.010 

0.020  ±  0.007 

±ve  control  (heat) 

0.014  ±0.006 

0.008  ±0.003 

- 

acceleration  in  the  experimental  impactor  is  1.7  x  106  m/s2  (see  Table  5.2).  The  com¬ 
puted  accelerations  arc  three  orders  of  magnitude  lower  than  the  critical  accelerations 
required  for  break-up.  The  computational  results  therefore  predict  that  break-up  of 
spores  is  not  achieved  under  these  conditions.  Experimental  values  of  fi  for  x  =  0.98 
closely  match  the  predictions  with  0.715  ±  0.235  fraction  of  the  cells  surviving  the 
impactor  operation.  Additionally,  for  x  =  0.11,  where  the  aerodynamic  shock  is 
strongest,  spores  reach  maximum  accelerations  of  5.9  x  10°  m/s2  and  5.0  x  109  in/s2 
for  0.5  and  1  /rm  spores,  respectively  (see  Table  5.2).  The  critical  acceleration  for 
0.5  fi m  spores  is  63%  greater  (less  than  an  order  of  magnitude)  than  the  computed 
accelerations  in  the  impactor.  The  computations  provide  insight  into  the  flow  struc¬ 
tures  in  the  impactor  system  and  are  not  exact  representations  of  the  real  fluid  flow. 
Therefore,  a  63%  difference  in  value  is  not  significant  and  is  reflected  in  the  experi¬ 
mental  fi  values  for  x  =  0.11  shown  in  Table  5.4  and  Fig.  5.10.  The  same  argument 
cannot  be  made  for  x  —  0.50  where  the  computed  amax  is  an  order  of  magnitude 
lower  than  ac,  yet  /;  shows  experimental  break-up  of  the  cells.  As  mentioned  in  Sec- 
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tion  5.3.3,  the  values  for  surface  tensions  are  calculated  either  from  turgor  pressure  or 
contact  angle  measurements  and  might  not  reflect  the  actual  surface  tensions  under 
shear  conditions.  The  surface  tension  values  of  spores  under  shear  can  change  the 
calculated  critical  accelerations  (see  Section  5.3.3).  However,  obtaining  these  values  is 
outside  the  scope  of  this  work.  Therefore,  the  computational  along  with  experimental 
results  provide  a  framework  for  assessing  the  effect  of  the  aerodynamic  shock  on  the 
spores  rather  than  an  exact  estimate.  For  both  0.5  and  1  pm  cells,  fi  at  x  —  0.11  is 
less  than  fi  at  x  =  0.50  with  p  >  95%.  For  the  non-shock  case  the  fi  is  significantly 
higher.  The  same  trends  were  observed  in  Chapter  4  with  E.  coli  [49]. 

As  discussed  in  Section  5.4.1,  B.  atropheus  spores  that  survive  break-up  in  the 
impactor  system  and  are  collected  in  the  deceleration  tube  are  more  likely  to  be 
alive  than  E.  coli  vegetative  cells  (L/5011  >  L/VQbQ).  The  data  indicates  that  the 
spore  coat  is  more  resistant  to  membrane  damage  if  the  particle  remains  intact  (see 
Remark  12).  The  damage  increases  as  the  x  decreases  while  in  vegetative  experiments 
the  damage  is  already  at  a  maximum  at  x  =  0.5  [49].  In  addition,  fi  is  equal  for  both 
X  =  0.11  and  x  =  0.5.  Although,  the  same  number  of  cells  are  breaking  up  the 
remaining  cells  are  weaker  as  the  acceleration  is  higher.  Operating  the  impactor  at 
lower  x  (i-eM  higher  accelerations)  can  give  more  insight  into  the  mechanism  of  break¬ 
up  for  both  vegetative  cells  and  spores.  We  expect  fi  values  to  be  lower  as  cells  that 
currently  survive  the  shock  get  weaker  and  break-up. 

Remark  12  The  accelerations  calculated  in  Section  5.3.2  are  centerline  a ccelera- 
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tions.  In  reality  the  bacterial  particles  experience  a  distribution  of  accelerations  with 
the  maximum  at  the  centerline.  For  x  =  0.5  the  accelerations  away  from  the  centerline 
can  be  at  values  where  the  spore  coat  of  intact  cells  is  not  affected. 

Remark  13  In  order  to  calculate  the  exact  value  for  fi ,  the  average  collection  effi¬ 
ciency  has  to  be  calculated  from  the  particle  size  distribution  of  the  spore  suspension. 
More  measurements  are  required  to  obtain  an  accurate  size  distribution.  For  exam¬ 
ple,  an  aerosol  particle  size  analyzer  could  be  used.  However,  the  focus  of  this  study 
is  to  show  the  range  of  possible  fi  values,  instead  of  a  very  accurate  value.  This 
range  is  sufficient  to  determine  the  effect  of  the  aerodynamic  shock  and  comparing 
the  experimental  results  to  computational  results. 

5.5  Conclusions 

The  aerodynamic  shock  was  experimentally  shown  to  cause  break-up  of  B.  atropheus 
spores  at  x  =  0.11  and  x  =  0.50.  The  study  provides  a  framework  that  yields  the 
parameters  needed  to  guide  designs  of  systems  that  can  neutralize  spore  aerosols 
at  their  release.  The  computational  model  predicts  the  break-up  of  the  spores  at 
X  =  0.11  where  accelerations  in  the  aerodynamic  shock  reach  values  higher  than 
the  ac  (>  3.9  x  10*0  m/s2).  The  experimental  results  confirm  the  predictions  with 
fi  =  0.030  ±  0.010.  The  model  also  predicts  the  survival  of  spores  at  accelerations  of 
106  rii/s2  by  operating  the  impactor  at  x  =  0.98.  This  prediction  was  also  confirmed 
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experimentally  with  ft  —  0.715  ±0.235  [50].  The  critical  accelerations  were  computed 
using  static  calculations  of  surface  tension.  When  surface  tension  values  under  shear 
conditions  become  available  for  B.  atropheus  spores,  the  predictions  can  be  adjusted 
using  the  framework  provided  in  this  work  [49,  51 
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Chapter  6 


Conclusions  of  the  dissertation 


It  has  been  shown,  for  the  first  time,  that  aerodynamic  shocks  can  neutralize  vegeta¬ 
tive  and  spore  forms  of  bacterial  aerosol.  The  study  provides  a  theoretical  framework 
to  predict  the  decelerations  experienced  by  a  bacterial  particle  as  it  passes  through  an 
aerodynamic  shock.  The  study  considers  the  break-up  of  a  particle  as  the  mechanism 
of  neutralization  which  is  more  stringent  than  loss  of  viability  or  infectivity  of  a  cell. 
The  model  predictions  can  be  enhanced  by  obtaining  more  accurate  surface  tension 
values  for  cells  under  fluid  flow. 

The  computational  framework  developed  in  Chapters  2  and  3  was  validated  ex¬ 
perimentally  using  two  test  bacteria.  In  Chapter  4,  vegetative  Gram-negative  E.  coh 
cells  were  found  to  break-up  under  conditions  producing  accelerations  higher  than 
the  critical.  In  Chapter  5,  Gram-positive  B.  atropheus  spores  were  also  found  to 
break-up  under  conditions  predicted  by  the  computational  model.  However,  spores 
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had  a  higher  percentage  of  live  cells  in  the  population  of  cells  that  did  not  break-up. 


These  results  can  be  used  in  two  ways.  First,  the  shock  strength  of  an  explosive 
device  can  be  designed  to  produce  the  decelerations  needed  for  bacterial  break-up  in 
the  context  of  defense  applications.  These  situations  can  include  the  design  of  explo¬ 
sives  that  minimize  collateral  effects  of  bacterial  aerosols  in  cases  where  a  biological 
stockpile  is  being  destroyed.  Second,  nozzles  can  be  designed  to  avoid  the  break-up  of 
cells  during  sampling  where  the  collection  of  intact  live  bacteria  is  needed  to  perform 
subsequent  analysis. 
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Appendix  B 


Fluent  and  Matlab  Codes 


The  appendix  contains  sample  FLUENT  and  Matlab  codes  used  to  solve  the  equations 
in  Appendix  A. 

************************  Sample  FLUENT  code  ********************************* 


START  OF  JOURNAL  FILE 


/file/ start-transcript  "C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases 

\P0_0992_Pl_0114\transcript l_build_model . trnn 

yes 


STEP  1:  Build  Model 


read  mesh 


/file/read- case  "C: \Fluent . IncXWorkXPatrickXMeshXnozzle.xdl^Xxdl^.e.vlO . rush" 
yes 


Scale 


/grid/scale 

0.001 

0.001 


;  change  unit  of  pressure 
/define/units 


VIII 


pressure 

atm 


define  solver 


/def ine/models/solver/pressure- based 
yes 


;  axis syme trie 
» 

/def ine/models/axi symmetric 
yes 


;  use  Green-Gauss  Node-Based  gradients 

j 

/solve/ set/gradient-scheme 
yes 


;  enable  energy  equation 

/def ine/models/ energy 

yes 

no 

no 

no 

yes 


turbulent  model 


/def ine/models/viscous  ke-standard 
yes 


Materials 


/def ine/materials/change- create  air 

air 

yes 

ideal-gas 

no 

no 

yes 

Sutherland 

three-coeff icient -method 

1 . 716e-05 

273.11 


IX 


110.56 

no 

no 

no 

no 

no 

no 


;  operating  conditions 

/define/ operating- conditions  operating-pressure 
0 

/def ine/ operating- conditions  gravity 
no 

/def me/operating- conditions  reference-pressure-location 

0 

0 


Boundary  conditons 


/define/boundary-conditions/pressure-inlet  pressure_inlet . 1 
no 

0.992 

no 

0.553 

no 

300 

no 

yes 

no 

no 

yes 

1 

1 


;  Boundary  conditions:  outlet 

/define/boundary-conditions/pressure-outlet  pressure.outlet . 2 
no 

0.114 

no 

300 

no 

yes 

yes 


X 


no 

1 

no 

1 

no 


Define  centerline 


/ surface/line-surface 
centerline 
0.062 
0 

0.0671 

0 


Define  Pitot 


/ surface/point-surface 

pitot 

0.0671 

0 


monitor->residual 


/ solve/monitors/ surface/ clear-monitors 

/solve/monitors/residual/check- convergence 

yes 

yes 

yes 

yes 

yes 

yes 

/ solve /monitors /residual/ convergence -criteria 
0.0001 
0.0001 
0.0001 
le-6 
0.0001 
0.0001 

/ solve/monitors/residual/plot 
yes 


Set  solvers 


/solve/set/p-v-coupling  1 


XI 


/solve/set/under-relaxat ion/body-force  1 
/solve /set /under-relaxation/density  1 
/solve/set/under-relaxation/epsilon  0 . 8 
/solve/set/under-relaxation/k  0.8 
/solve/set/under- relaxation/mom  0 . 7 
/solve/ set/under-relaxation/pressure  0 . 3 
/solve /set /under-relaxation/temperature  1 
/solve/set/under- relaxation/turb-viscosity  1 


First-order  discretization 


/solve/set /discretization-scheme/density  0 
/ solve/set/discretization-scheme/epsilon  0 
/solve/set/discretization-scheme/k  0 
/sol ve/set /discretization-scheme /mom  0 
/solve/set/discretizat ion-scheme /pressure  10 
/solve /set /discretization- scheme /temperature  0 
/file/write-case  MC : \Fluent . Inc\Work\Patrick 
\SimQ04_Two_cases\P0_0992_Pl_0114\case0 . cas" 
yes  ;  OK  to  overwrite 
/file/ st op- transcript 

/f ile/start-transcript  HC : \Fluent . Inc\Work\Patrick 

\Sim004_Two_cases\P0_0992_Pl_0114\transcript2_ini . trnM 
yes 


STEP  2:  fmg  initialization 


initialize 


/solve/ initialize/ compute-default/pressure-inlet  pressure_inlet . 1 
/solve/ initialize/ initialize -flow 
/ solve/ initialize/ set-f mg-initialization 
5 

0.001 

100 

0.001 

100 

0.001 

100 

0.001 

500 

0.001 

500 

0.75 


XII 


yes 

solve/initialize/fmg-initialization 

yes 


;  Save  for  debug 

wed  "C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases\P0_0992_Pl_0114\ini . cas 
yes 

/file/stop-transcript 

/file/start-transcript  "C : \Fluent . Inc\Work\Patrick 

\Sim004_Two_cases\P0_0992_Pl_0114\transcnpt3_F0 . trn" 
yes 

;  STEP  3:  Solve  flow  equation  only 


Disable  energy  and  turbulent  equation 


/solve/set/equations/ke  no 
/solve/set/ equations/temperature  no 


iterate  2000  times 


/solve/iterate  2000 


save 


wed  "C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases\P0_0992_Pl_01 14\F0 . cas" 

yes  ;  OK  to  overwrite7 

/file/stop-transcript 

/file/ start-transcript  "C : \Fluent . Inc\Work\Patrick 

\Sim004_Two_cases\P0_0992_Pl_0114\transcript4_CP_SKE. trn" 
yes 


STEP  4:  Solve  coupled  equations 


Enable  energy  and  turbulent  equation 


/solve/set/equations/ke  yes 
/solve/set/ equations/temperature  yes 


Solve 


/solve/iterate  2000 


XIII 


Save 


wed  "C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_P1_01 14\F0_CP_SKE . cas M 
yes  ;  OK  to  overwrite? 

/f ile/ stop-transcript 

!  RD  "C: \Fluent . Inc\Work\Patrick\Sim004_Two_cases\P0_0992_Pl_0114\SD" 
1  MD  MC: \Fluent . Inc\Work\Patrick\Sim004 JTwo_cases\P0_0992_Pl_01 14\SD" 
/f ile/st art- transcript  "C: \Fluent . Inc\Work\Patrick 

\Sim004_Two_cases\P0_0992_Pl_01 14\SD\transcript5_SD. trn" 
yes 

;  STEP  5:  Second-order  discretization 


Second-order  discretization 


/ sol ve/set /discretization-scheme/density  1 
/solve/set/discretization-scheme/k  1 
/solve/set/discretization-scheme/epsilon  1 
/solve/set/discretization-scheme/temperature  1 
/solve/set/di scretizat ion-scheme /mom  1 
/solve/set/discretizat ion- scheme /pres sure  12 


Change  under-relaxation  factor 


/solve/set/under-relaxat ion/body-force  1 
/solve/set /under-relaxation/density  1 
/solve/set/under-relaxation/epsilon  0 . 5 
/solve/set/under-relaxation/k  0.5 
/solve/set/under-relaxation/mom  0 . 5 
/solve/ set/under-relaxation/pressure  0 . 2 
/solve/ set/under-relaxatlon/temperature  1 
/solve/set/under-relaxation/turb-viscosity  1 


Solve 


/solve/iterate  2000 


Save 


wed  "C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_P1_01 14\SD\F0_CP_SD . cas M 
yes  ;  OK  to  overwrite'’ 


XIV 


/file/stop-transcript 

!  RD  "C : \Fluent . Inc\Work\Patrick\Sim004_Tvo_cases\P0_0992_Pl_0114\ADl " 
!  MD  "C: \Fluent . Inc\Work\Patrick\Sim004_Two_cases\P0_0992_Pl_0114\ADl" 
/file/start-transcript  "C : \Fluent . Inc\Work\Patrick 

\Sim004_Two_cases\P0_0992_Pl_0114\ADl\transcript6_ADl . trn" 
yes 


STEP  6:  ADAPTION 


;  Adaption 

/ad apt /ad apt -to -gradients 
pressure 
space-gradient 
1 

no 

0.7 

0 

no 

0 

0 

0 

yes 


iterate 


/solve/iterate  2000 


Save 


wed  "C : \Fluent . Inc\Work\Patrick\Sim004 _Two_cases 
\P0_0992_P1_0114\AD1\F0_CP_SD_AD1 . casM 
yes 

/file/stop-transcript 

!  RD  "C:\Fluent.Inc\Work\Patrick\Sim004jryo_cases\P0_0992_Pl_0114\AD2" 
'  MD  MC : \Fluent . Inc\Work\Patrick\Sira004_Two_cases\P0_0992_Pl_0114\AD2H 
/file/ start-transcript  "C : \Fluent . Inc\Work\Patrick 

\Sim004_Two_cases\P0_0992_Pl_0114\AD2\transcript7_AD2 . trn" 
yes 

;  STEP  7:  ADAPTION  2 

red  "C : \Fluent . Inc\Work\Patrick\Sim004 JTwo_cases 
\P0_0992_P1_0114\AD1\F0_CP_SD_AD1 . cas " 


XV 


Adaption 


/ad apt /adapt -to -gradients 
pressure 
space-gradient 
1 

no 

0.7 

0 

no 

0 

0 

0 

yes 


iterate 


/solve/iterate  2000 


Save 


wed  "C: \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_P1_01 14\AD2\F0_CP_SD_AD2 . cas" 
yes 

/f ile/ st op- transcript 

!  RD  "C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases\P0_0992_Pl_0114\AD3 
!  MD  "C: \Fluent . Inc\Work\Patrick\Sim004_Two_cases\P0_0992_Pl_0114\AD3 
/file/start-transcript  "C: \Fluent . Inc\Work\Patrick 

\Sim004_Two_cases\P0_0992_Pl_0114\AD3\transcript8_AD3 . trn" 
yes 

;  STEP  8:  ADAPTION  3 


red  "C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_P1_0114\AD2\F0_CP_SD_AD2 . casM 


Adaption 


/adapt /adapt -to-gradients 
pressure 
space-gradient 
1 

no 

0.7 


XVI 


0 

no 

0 

0 

0 

yes 


iterate 


/solve/iterate  2000 


Save 


wed  "C : \Fluent , Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_P1_01 14\AD3\F0_CP_SD_AD3 . cas  M 
yes 

/file/ stop-transcript 


Postprocessing 


Second  order  discretization 


Massflow  rate  across  inlet  and  outlet 


red  "C : \Fluent . Inc\Work\Patrick\Sim004_Two_ cases 
\P0_0992_P1_0114\SD\F0_CP_SD . cas" 
/report/fluxes /mass -flow 
no 

pressure.inlet .  1 
pressure.outlet . 2 
0 
yes 

"C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_Pl_0114\SD\massf low.SD . f lp" 
no 
yes 


Maximum  mach  number  along  center  line 


/report /surf ace- integrals /f acet-max 
centerline 
0 

mach- number 


XVII 


yes 

"C: \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_Pl_0114\SD\mach_SD . srp" 
no 
yes 


;  Pitot  pressure 

/report/ surf ace-integrals/f acet-avg 
pitot 

0 

pressure 

yes 

"C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_Pl_0114\SD\Pp_SD . srpM 

no 

yes 


Adapt  once 


Massflow  rate  across  inlet  and  outlet 


red  "C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_P1__0114\AD1\F0_CP_SD_AD1 . cas" 
yes 

/report /f  luxes/mass- f  low 
no 

pressure_inlet . 1 
pressure_outlet . 2 

0 

yes 

"C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992 _P1_01 14\ADl\massf low.ADl . f lp” 
no 
yes 


Maximum  mach  number  along  center  line 


/ report/ surface-integrals/facet-max 
centerline 
0 

mach-number 

yes 

"C : \Fluent . Inc\Work\Patr ick\Sim004_Tvo_cases 


XVIII 


\P0_0992_Pl_0114\ADl\mach_ADl . srp" 
no 
yes 


Pitot  pressure 


/report/ surf ace- integrals/ facet -avg 
pitot 
0 

pressure 

yes 

"C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_Pl_0114\ADl\Pp_ADl . srp" 
no 
yes 


Adapt  twice 


Massflow  rate  across  inlet  and  outlet 


red  "C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_P1_01 14\AD2\F0_CP_SD_AD2 . casH 
yes 

/report /f luxe s/mass -flow 
no 

pressure_inlet . 1 
pressure_outlet . 2 
0 

yes 

MC:\Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_P1_01 14\AD2\massf low_AD2 . f lp" 
no 
yes 


Maximum  mach  number  along  center  line 


/report/ surf ace -integrals /facet -max 
centerline 
0 

mach-number 

yes 

,#C :  \Fluent .  Inc\Work\Patrick\Sim004 _Two_cases 
\P0_0992_Pl_0114\AD2\mach_AD2.srp" 
no 


XIX 


yes 


;  Pitot  pressure 

/report/ surf ace-integrals/f acet-avg 
pitot 
0 

pressure 

yes 

”C : \Fluent . Inc\Work\Patr ick\Sim004_Two_cases 
\P0_0992_P1_01 14\AD2\Pp_AD2 .  srp" 
no 
yes 


Adapt  the  third  time 


Massflow  rate  across  inlet  and  outlet 


red  "C: \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_P1_01 14\AD3\F0_CP_SD_AD3 . cas” 
yes 

/report /f luxe s/mass-f low 
no 

pressure_inlet . 1 
pressure_outlet . 2 
0 
yes 

"C : \Fluent . Inc\Work\Patrick\Sim004_Two_cases 
\P0_0992_P1_01 14\AD3\massf low_AD3 . f lp" 
no 
yes 


;  Maximum  mach  number  along  center  line 

/report/ surf ace-integrals/f acet-max 
centerline 
0 

mach-number 

yes 

"C: \Fluent . Inc\Work\Patr ick\Sim004_Two_cases 
\P0_0992_P1_01 14\AD3\mach_AD3 . srp” 
no 
yes 


XX 


;  Pitot  pressure 


/report /surface- integrals /facet -avg 
pitot 
0 

pressure 

yes 

"C : \Fluent . Inc\Work\Patr ick\Sim004_Two_cases 
\P0_0992_Pl_0114\AD3\Pp_AD3 .  srp" 
no 
yes 

********************  Sample  Matlab  code  ***************************** 
function  [acc , acc_max]  =  Acceleration(data_g, data.p , interp.method, gap) 

%  [acc ,acc_max]  =  Acceleration(data_g,data_p , interp_method , gap) 

7,  Calculates  maximum  acceleration 

% 

7.  INPUTS: 

y,  data_g:  returned  from  ReadGasPhase () 

*/,  data_p:  returned  from  ReadParticle () 

*/,  mterp_method :  same  as  interplO,  e.g.,  1  linear  *  ,  9  cubic  *  . 

'/,  gap:  a  two-element  vector  indicating  the  region  in  which  to  search  for 
'/,  find  maximum  acceleration.  Unit  is  meter. 

y. 

y. 

time  =  data_p{l}; 
x_p  =  data_p{2}; 
u_p  -  data_p{5}; 
x_g  =  data_g{l} ; 
u_g  =  data_g{2} ; 

u_g_new  =  interpl (x_g,u_g ,x_p , mterp_method) ; 
delta_u  =  (u_g_new  -  u_p) ; 

•/.  Acceleration  calculation  method  #1 
acc  =  zeros (length(u_g_new) -1 , 1 ) ; 
for  i  =  1 : (length(u_g_new)-l) 

acc(i)  =  (delta_u(i+l) -delta_u(i) )/(time (i+l)-time (i) ) ; 

end 

acc_abs  =  abs(acc); 

'/,  acc_max  is  the  maximum  acceleration  between  the  gap 

index  =  (x_p(l : (length(x_p) -1) )  >  gap(l))  k  (x_p( 1 : (length(x_p) -1) )  <  gap(2)); 
acc.max  =  max (acc_abs( index) ) ; 
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Appendix  C 


Distribution  List 


1.  DTIC-OCP,  8725  John  J.  Kingman  Road  Suite  0944,  Fort  Belvoir,  VA  22060- 
6218,  two  copies 

2.  DTRA,  ATTN:  NTRM,  6200  Meade  Road,  Ft.  Belvoir,  VA  22060-5264,  one 
copy 

3.  DTRA,  ATTN:  BDCT,  6200  Meade  Road,  Ft.  Belvoir,  VA  22060-5264,  one 
copy 

4.  DTRA,  ATTN:  BMCI,  6200  Meade  Road,  Ft.  Belvoir,  VA  22060-5264,  one 
copy 
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